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"We are at the very beginning of time for the human race. It is not unreasonable that 
we grapple with problems. But there are tens of thousands of years in the future. Our 
responsibility is to do what we can, learn what we can, improve the solutions, and pass 
them on." 

Richard Feynman 
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Abstract 



In this thesis I investigate aspects of two problems. In the first part of this thesis, I will 
investigate how an effective field theory can be constructed. One of the most fundamental 
questions in physics is how new degrees of freedom emerge from a fundamental theory. In 
the Hamiltonian framework this can be rephrased as finding the correct representation for 
the Hamiltonian matrix. The similarity (not essentially unitary) renormalization group 
provides us with an intuitive framework, where a transition from a perturbative region 
to a non-perturbative one can be realised and physical properties can be computed in a 
unified way. In this context, we have shown that the well-known coupled-cluster many- 
body theory techniques can be incorporated in the Wilsonian renormalization group to 
provide a very powerful framework for construction of effective Hamiltonian field theo- 
ries. Eventhough the formulation is intrinsically non-perturbative, we have shown that a 
loop-expansion can be implemented. 

The second part of my thesis is rather phenomenologically orientated. In this part, 
I will employ an effective field-theoretical model as can be constructed by means of the 
techniques of the first part of my thesis, a quark-confining non-local Nambu-Jona-Lasinio 
model and study the nucleon and diquarks in this model. For certain parameters the 
model exhibits quark confinement, in the form of a propagator without real poles. After 
truncation of the two-body channels to the scalar and axial- vector diquarks, a relativistic 
Faddeev equation for nucleon bound states is solved in the covariant diquark-quark pic- 
ture. The dependence of the nucleon mass on diquark masses is studied in detail. We 
find parameters that lead to a simultaneous reasonable description of pions and nucleons. 
Both the diquarks contribute attractively to the nucleon mass. Axial-vector diquark corre- 
lations are seen to be important, especially in the confining phase of the model. We study 
the possible implications of quark confinement for the description of the diquarks and the 
nucleon. In particular, we find that it leads to a more compact nucleon. 
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General remarks 



This thesis is organized as follows: in the first three chapters, we concentrate on renor- 
malization group methods in Hamiltonian framework. In chapter 1, we introduce the 
coupled-cluster theory. In chapter 2, we show how renormalization group can be em- 
ployed in the context of the coupled-cluster theory. In order to highlight the merits and 
the shortcomings of our approach over previous ones, in Sec. 2.2, we review different RG 
methods in Hamiltonian framework. Different aspects of our approach is introduced in 
sections 2.4-2.8. In sections 2.9 and 2.10, as illustrative examples, we apply our formu- 
lation on the theory and an extended Lee model. In chapter 3, we pursue a different 
approach for the renormalization of the many-body problem. We show that a combination 
of the coupled-cluster theory and the Feshbach projection techniques leads to a renormal- 
ized generalized Brueckner theory. 

In the second part of this thesis, we investigate the baryon structure in a chiral quan- 
tum chromodynamics model based on the relativistic Faddeev approach. In sections 4. 1 
and 4.2, we introduce the most important properties of QCD which are needed for the 
modelling of hadrons. In sections 4.3 and 4.4, we show how an effective low-energy field 
theory can be constructed from the underlying QCD theory. In chapter 5, we introduce 
alternative field theoretical approaches for describing baryons, such as Skyrme models, 
bag models and diquark-quark models in the context of the relativistic Faddeev approach. 
Finally, in chapter 6, we study baryons based on the diquark-quark picture in a quark- 
confining non-local NJL model. The non-local NJL model is introduced in Sec. 6.2. In 
Sec. 6.3, we discuss the pionic sector of the model. In Sec. 6.4 the diquark problem is 
solved and discussed. In Sec. 6.5 the three-body problem of baryons is investigated. The 
numerical technique involved in solving the effective Bethe-Salpeter equation is given 
and the results for three-body sector are presented. 
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Chapter 1 



Basic structure of the coupled-cluster 
formalism 

1.1 Introduction 

In order to understand fully the properties of quantum many-body systems, various meth- 
ods have been developed which aim to go beyond perturbation theory. One of the simplest 
approaches has been the so-called configuration-interaction method which diagonalises 
the Hamiltonian in a finite subspace of the full many-body Hilbert space. An exten- 
sion of this method has been introduced via various versions of coupled-cluster meth- 
ods [1, 2, 3, 4]. The coupled cluster method (CCM) in its simplest form originated in 
nuclear physics around forty years ago in the work of Coester and Kiimmel [1]. The 
configuration-interaction method (CIM), and various version of coupled-cluster methods: 
normal coupled cluster method (NCCM)[1, 2] and the extended coupled cluster method 
(ECCM)[3, 4] form a hierarchy of many-body formulations for describing quantum sys- 
tems of interacting particles or fields [5]. They are denoted generically as independent- 
cluster (IC) parametrizations, in the sense that they incorporate the many-body correla- 
tions via sets of amplitudes that describe the various correlated clusters within the inter- 
acting system as mutually independent entities. The intrinsic non-perturbative nature of 
the methods is considered to be one of their advantages which make them almost univer- 
sally applicable in many-body physics. The IC methods differ from each other in the way 
they incorporate the locality and separability properties; in a diagrammatic language they 
differ in their linking properties. Each of the IC methods has been shown [4] to provide an 
exact mapping of the original quantum mechanical problem to a corresponding classical 
mechanics in terms of a set of multiconfigurational canonical field amplitudes. The merit 
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of IC has been outlined in Ref. [5] and literature cited therein. 

1.2 Formalism 

In this section we concentrate on the NCCM and the ECCM from a formal viewpoint. 

Exponential structures arise frequently in physics for similar underlying fundamental 
reasons. For example, in the Ursell-Mayer theory in statistical mechanics, in the Gold- 
stone linked-cluster theorem [6] and the Gell-Mann and Low theorem [7]. The complex- 
ity of the vacuum (the ground state) of an arbitrary many-body system in the NCCM 
parametrization [1, 2, 5] is expressed by an infinite set of correlation amplitudes {sj, sj} 
which have to be determined by the dynamics, 

\^) = K{t)e''\^l>o), 

Here K{t) is a time-dependent scale factor. The coefficients sj and s> are time depen- 
dent. The intermediate normalization condition (V'lV') — 1 is explicit for all times t. 
We restrict ourselves to the non-degenerate system, so that the exact states of the sys- 
tem may sensibly be refereed to some suitably chosen single reference state denoted as 
\^o)- The state \iJjq) is a ground state, e.g. a special (Hartree) bare vacuum in quan- 
tum field theory (QFT). The function can be chosen rather generally, but is tied to 
the choice of generalized creation operators {C]}; the state jt/'o) is annihilated by {Cj} 
V/ 7^ (where Co = 1, the identity operator) and is a cyclic vector in the sense that 
the algebra of all possible operators in the many-body Hilbert space H is spanned by the 
two Abelian subalgebras of creation and annihilation operators defined with respect to 
it. In this way we can define proper complete orthonormal sets of mutually commuting 
configuration creations operators {C|} and their Hermitian adjoint counterparts {C/}, 
defined in terms of a complete set of many -body configuration {/}. These are, in turn, 
defined by a set-index /, which labels the cluster configuration created by Cj with re- 
spect to the reference state li^o). Therefore, {/} defines a subsystem or cluster within the 
full system of a given configuration and the actual choice of these clusters depends upon 
the particular system under consideration. We assume that the creation and annihilation 
subalgebras and the state |^o) are cyclic, so that all ket states in the Hilbert space H can 
be constructed from linear combinations of states {C'||V'o)}; and for the bra states with 
respect to states {{ipolCi}. It is well-known in the many-body application that the above 
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parametrization Eq. (1.1), guarantees automatically proper size-extensivity (see section 
1.3) and conformity with the Goldstone linked-cluster theorem [6] to all levels of trunca- 
tion. In contrast, the configuration-interaction method is size extensive and linked only 
in the full, infinite-dimensional space [5]. The NCCM parametrization of bra- and ket- 
states Eq. ( 1 . 1 ) , in its asymmetrical (independent) form, does not manifestly preserve their 
Hermitian congugacy, hence we have here a biorthogonal formulation of the many-body 
problem. However, this is the most reasonable parametrization if one is to preserve the 
canonical form of the equations of motion with respect to phase space {s/, §/} and the 
Hellmann-Feynman theorem ^ [4, 5]. Nevertheless, non-hermiticity is negligible if the 
reference state and its complement are not strongly correlated [40]. We may hope that 
S and S are small (in a somewhat ill-defined non-perturbative sense), in other words, we 
may require that some of the coherence has already been obtained by optimizing the ref- 
erence state. (This can be done, for example, by a Hartree-Bogolubov transformation). 
Then the remaining correlations can be added via the CCM^. Therefore, defining a good 
reference state can in principle control the accuracy of CCM. 

In the CIM, one defines the ket and bra states as follows, 

l^) = ^l^o); {^\ = {i^o\F, 

^ = E/^^1; F = E//C,, (1.2) 

I I 

where the normalization condition = 1 can not be imposed trivially. Although the 

CIM has a simpler parametrization than the CCM, it does not satisfy the size-extensivity 
after truncation and contains unlinked pieces emerging from the products of noninteract- 
ing subclusters. A formal relation between the CIM and coupled cluster theory will be 
demonstrated in section 1.3. 

While all ground-state expectation values {%jj\A\il)) = A{sj, sj) and amplitudes {s/} 
are linked in NCCM, the amplitudes {s7} contain unlinked terms. This is resolved in the 
ECCM [3, 5] where we reparametrize the Hilbert space such that all basic amplitudes are 



'According to Hellmann-Feynmaii theorem, if we perturb the Hamiltonian, H ^ H' = H + XA (where 
A is infinitesimally small quantity) then the ground state energy changes as Eq ^ Eo + XdEo /dX + 0(A^ ) 
with dEo/dX = {ipldH/dXl^jj}. 

^The application of this procedure to two-dimensional theory has been shown to give a stable result 
outside the critical region [41]. 
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linked, 

1^) = K{t)e'\i;o); = J^i^o\e^^-', 

E\^o) = Qeh-^S\^o); g = 1 - \^o){ipo\, 

E = ^a;C]; E=^a,C,. (1.3) 

The inverse relationships between the ECCM and the NCCM counterparts are given by 
ai = (V'olC/^^lVo), SI = (V'olC/e^SlV'o). (1-4) 

The ECCM amplitudes {a/, o-/} are canonically conjugate (which comes from time de- 
pendent variation)''. It must be clear that this parametrization for the ECCM is not unique, 
however it is complete and sufficient to specify the ECCM phase space. It has been 
pointed out earlier that this choice of parametrization is the most convenient one [3]. The 
ECCM is believed to be the unique formulation of quantum many-body with full locality 
and separability at all levels of approximation. The individual amplitudes {sj, §/} (or 
{o'l, CT/}) in both parametrizations are determined independently by solving an infinite 
set of non-linear equations which emerge from the dynamics of the quantum system. In 
practice, one needs to truncate both sets of coefficients. A consistent truncation scheme 
is the so called SUB(n) scheme, where the n-body partition of the operator {sj, sj} (or 
{aj, CT/}) is truncated so that the general set-index {/} contains up to n-tuple excitation 
(e.g., of single particle for bosonic systems with a reference state). Determination of the 



^The equation of motion for the ECCM amplitudes can be obtained from a variational principle [4] by 
requiring the action-like functional 

A = j dt{i)\id/dt- H{t)\ip), (1.5) 
to be stationary against small variations of amplitudes. After some straightforward algebra, one can obtain 



A = j dt -i'Y^ ^i(Ti - {H) 



1^0 

The stationary conditions lead us to the pair of equations of motion, 



(1.6) 



S{H) 5{H) 
i(Ti = tai = (1.7) 

0(7/ oaj 

Therefore, {aj, aj} are obviously canonically conjugate to each other in the usual terminology of classical 
Hamiltonian mechanics. We will elaborate more on this property of the coupled-cluster theory in the context 
of the renormalization group in the chapter 2.7. 
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amplitudes corresponds to summing infinite sets of diagrams which in perturbation lan- 
guage take into account arbitrary high-order contributions in the coupling constant, there- 
fore the NCCM and ECCM are not an expansion in this coupling constant. This property 
demands a precise truncation scheme for hierarchies without losing the renormalizability 
of the theory. We will consider this problem in the next chapter. 



In the following, we confine our consideration to a real Klein-Gordon field, as an 
example. The configuration operators are specified by {Cj — > ak^...akj, C\ a\^...a\^} 
with subalgebra 

K, 4'] = ^kk'^ K' ^k'] = 0. (1.8) 

The corresponding S and S operators are 

n=l qi,..qn 
S ^ l + ^Sn, Y ^Sn{qi,..qn)aqi-aq„, (1.9) 

n. 

n=i qi,--qn 

and iV'o) is the Fock vacuum. The individual amplitudes §„} which describe excita- 
tions of n Fock particles have to be fixed by the dynamics of the quantum system. Using 
Fock states in QFT has often been ambiguous due to problems connected with Haag's 
theorem [9] (Haag's theorem says that there can be no interaction picture - that we can- 
not use the Fock space of noninteracting particles as a Hilbert space - in the sense that 
we would identify Hilbert spaces via field polynomials acting on a vacuum at a certain 
time). It is well-known that the algebraic structure, does not, in general, fix the Hilbert 
space representation and therefore dynamical considerations are required. The IC formal- 
ism invokes dynamics rigorously. The dynamical principle to fix the physical vacuum is 
Poincare invariance, 

H\^) = PIV') = LIV') = K\^) = 0, (1.10) 

where H, P, L and K are generators of the Poincare group (Hamiltonian, momentum, an- 
gular momentum and boost operators, respectively). The excited states are no longer in- 
variant under these symmetry operations and the spectrum can be obtained in an extended 
version of IC [10]. By putting the Ansatz Eqs. (1.1,1.9) into the conditions Eq. (1.10) one 
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n 

{qi,...,qn\e-^Pe^\ipo) = {qi, ...,qn\{P + [P, S])\iPo) = [ g^] g„) = 0, 

i=l 

(1.11) 

n 

{^o\Se~'^Pe^\qi,...,qn) = {^o\S{P + [P, S])\qi, ...,qn) = [ ^ ■-, ^n) 

i=l 

+ E ^"^-^^^^-^-^-^ [E = 0, (1.12) 

m 1=1 

(gi, gf„|e~'^Le^|V'o> = ^nMoc)] r,, , gn(gi, ^n) = 0. (1.13) 

In the same fashion one can impose the condition {'4>Q\Se~^ Le^\q-i, ...,qr^ — 0, which 
leads to the following condition, having made used of the equation (1.13): 

" d 

E ^ijK?a)j^^^Sn(gi, qn) = (1.14) 

Similarly for Hamiltonian and boost operators we have 

(gi, gje-^i/e^l^o) = (^ol^e-^^^e^ki, Qn) = 0, (1.15) 
{qi, g„|e-^Xe^|Vo) = (MSe-^ Ke^lQi, Qn) = 0. (1.16) 

The equations (1.11,1.12,1.13,1.14) lead us to 

n n 

Sn{qi,--;qn) = S^^qi\sn{{qi.qj}), 5„(gi, g„) = S^'^q^^Sn{{qi.qj}), 

1=1 1=1 

(1.17) 

which means that {s/, 5/} depend on scalar quantities only and momentum is preserved. 
The same result can be obtained for the ECCM parametrization. To complete the deter- 
mination of phase space, we use the energy hierarchy Eq. (1.15), where i7 is a normal 
ordered Hamiltonian and the vacuum energy vanishes. Conceptually, it is evident that 
Eq. (1.15) suffices to fix all amplitudes without invoking Eq. (1.16), but explicit verifica- 
tion to all orders seems to be impossible. 



The fully linked feature of the e ^He^ term in Eq. (1.15) can be made explicit by 
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denoting it as {He^}c, which can be written as a set of nested commutators, 

e-^He^ = {He^}c ^H+[H,S] + ^[[H, S],S] + ... . (1.18) 

This procedure is still rigorous. Lorentz symmetry, stability and causality are examples 
of features normally expected to hold in physical quantum field theories. In renormalized 
QFT stability and causality are closely intertwined with Lorentz invariance. For example, 
stability includes the need for energy positivity of Fock states of ordinary momenta, while 
causality is implemented microscopically by the requirement that observables commute 
at spacelike separation [11], so-called microcausality. In addition, both are expected to 
hold in all inertial frames. A stable and causal theory without Lorentz symmetry could in 
principle still be acceptable [12]. In the framework of many-body theory, medium contri- 
butions always affect the local and global properties of hadrons especially at high density. 
In dense matter, the Pauli principle and cluster properties can affect causality since they 
restrict the permissible process in a scattering reaction. To understand these effects we 
need firstly to consider if our formalism itself can in principle preserve the causal struc- 
ture of given physical system. Let us introduce a new set which are connected to 
previous {a^, a[,} defined in Eq. (1.8) via a generalized Bogolubov transformation, 

bl = Akial + Bkiai + Di, (1.19) 

This is the most general linear transformation, which preserves commutator relations (bk 
is an annihilation operator and b^ is the Hermitian conjugate of bk), 

[bk, bl,] = Skk', [bk, bk'] = 0, (1.20) 

provided that AA^ — BB^ — 1. This obviously preserves the commutator relations for 
boson fields and preserves microcausality explicitly. The bare "a vacuum" defined by 
afe|V'o)a — 0, is replaced by a bare "6 vacuum" which satisfies 6^1 "00)6 — 0. Using 
Thouless theorem [14], 1^);, can be written as 

|V^);, = iV-l/Vi+^^|V^)„, (1.21) 

with 

•^i = E '^^(^' ^' S2 = J2 l^Slk'iA, B, D)aial, (1.22) 

k kk' 

where and S'^ are known functions of the matrices A and B and the vector D [14]. 
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It is obvious that is a low-order approximation to the CCM wave function Eq. (1.1) 
or (1.3). The above-mentioned parametrization of vacuum wave function can be general- 
ized to the CCM wave function, however it may require a nonlinear transformation from 
which it can be constructed. The inspiration for this transformation can be taken from the 
extension IC formulation for excited states [10] 

oo 

bl = e'F{k), F{k) = 4 + ^ k)al...al^_^, (1.23) 

n=3 

where the correlation operator S is known from Eq. (1.9) and F is a new amplitude 
which includes momentum conservation. This new amplitudes has to be determined and 
it changes the momentum of the bosons before creation. It is not hard to derive equations 
for the energy spectrum which lead us to A^-body effective Hamiltonians [15] and yields 
the folded diagrams of degenerate many-body perturbation theory [16]. This nonlinear 
transformation manifestly invalidates the commutator relation Eq. (1.20) and accordingly 
the microcausal commutation relations of the boson field. It should be noted that causality 
of the underlying theory can not be fully determined at this level and one needs to take 
into account the dynamics of the underlying quantum system. Causality in the context of 
IC formulation can be ensured by requiring 

{i^mX), (t>{Y)M =0, {X- Yf < 0, (1.24) 

where Y denotes space-time coordinates (y, yo) and 4>{Y) is Klein-Gordon field operator 
which is expressed in terms of Fock space operators by 

<t>{Y) = Y.mkak + il{Y)al (1.25) 

k 

where ik{Y) form a complete orthonormal set of states. Eq. (1.24) can be evaluated 
explicitly by using the following identity [4, 5] 

/Tir. Till ,\ r , „i OA OB OA OB 

im,Bm = lA,B]=j:—--——, (1.26) 
where {xj, xj} are the canonical coordinate and momenta of the NCCM or the ECCM 
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parametrization. By making use of Eqs. (1.9,1.18) and (1.25) one can find 

miVm = 0(F) = ^ ^ efc(>^)^n-i(gi, qn-i)xniqi, qn-1, k) + il{Y)x,{k), 

n=l k 

(1.27) 

where {xn — > s„, (7„, x„ ^ s„, On} and we define — 1. By exploiting Eq. (1.26) one 
can show 

WXhW)] = Y.\.^k{X)ii{Y) - il{X)i,,{Y)] = A{X -Y)- A{Y - X). (1.28) 

k,k' 

The function A{X — Y) is the Pauli- Jordan function [11] defined for field operators ex- 
panded in plane- wave basis. When [X — Y)"^ < 0, we can perform a Lorentz transfor- 
mation on the second term, taking [X — Y) — > —{X — Y). The two terms are therefore 
equal and cancel to give zero, hence, Eq. (1.24) is satisfied. It should be noted that for 
a general quantum system by considering just the dynamics of the underlying system, 
determining the amplitudes {xi, xi} and verifying Eq. (1.24) at every level of truncation 
one can ensure causality. Obviously this might introduce a lower limit of truncation in a 
consistent SUB(n) scheme which contains causality. 

In relativistic quantum mechanics there is another distinct type of causality, the fact 
that there is a well-posed initial value problem, so-called Cauchy causality. In the local 
field theory the Poincare invariance implies the existence of a unitary representation of the 
Poincare group that acts on the Hilbert space. In other words, the transformed final state 
is uniquely determined by time evolving the transformed initial state. Therefore, Cauchy 
causality is a consequence of Poincare invariance and is independent of any consideration 
concerning microcausality. 

One of the standard approaches in the nuclear many-body theory is to introduce an ef- 
fective interaction and reduce the full many-body problem to a problem in a small model- 
space. As an example for the meson-nucleon system below threshold, it was shown that 
one has either hermitian effective operators with all desirable transformation properties 
under the Poincare group or non-hermitian ones obtained by the CCM with desirable sim- 
plicity [13]. Thus effective operators introduced by the CCM can not manifestly preserve 
Cauchy causality, at any level of truncation. However, it economically disentangles and 
reduces the complexities of the many-body problem. It is notable that in some cases, the 
violation of relativistic invariance due to truncation is consistent with errors introduced 
via the approximation [14]. In practice one should not expect that Poincare invariance 
holds for an approximation, and forcing the approximated wave function to obey Lorentz 
invariance might lead to inconsistent results. Having said that, we will show that in the 
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context of coupled-cluster renormalization group framework, one can in principle define 
a truncation scheme where Poincare invariance is preserved. 

1.3 Linked-cluster theorem and Wightman functionals 

In this section we proceed from a formal view to show the relation between the CIM and 
the coupled cluster method. We will show that the CCM is a natural reparametrization 
of the CIM which incorporates the full size-extensivity. For this purpose we exploit the 
reconstruction theorem [17], well known in axiomatic Quantum field theory [18]: 

Denote the state vector of the vacuum by (in the language of the CCM \n) is the 
full interacting ground state |^)). The physical vacuum expectation value of products of 
local fields 

w'\Xi...Xn) = m{xi)..4{Xnm, (1.29) 

are tempered distributions'^ over M''". These w"- are called Wightman distributions. If 
the hierarchy of distributions = 0, ...) is known then the Hilbert space can be con- 
structed. This is the so-called reconstruction theorem [17]. 

Now consider a configuration (a;i, Xn) consisting of several clusters. A cluster here 
is a subset of points so as all the points in one cluster have a large space-like separation 
from all the points in any other clusters. One may expect that in infinite separation be- 
tween the clusters, one has 

w''{xi,...,Xn) ~ Ylw'''^{yr,l,-yr,nr), (1-30) 

r 

the index r and rir denote the cluster and number of points in the r-th cluster, respectively. 
?/r,nfc indicates a subset of the Xk belonging to this cluster. Thus we require that in the 
vacuum state the correlation of quantities relating to different regions decreases to zero as 
the space-like separation of the regions increases to infinity. Therefore it makes sense to 
introduce another hierarchy of functions w^, 

W^{x,,...,Xr,) ^J^U'^T^yr^k)- (1-31) 

p r 

Here p denotes a partition of the set of points Xi into subsets labeled by the index r. 
The sum is over all possible partitions. The objects wJ^ are called truncated functions or 

^Tempered distributions generalize the bounded (or slow-growing) locally integrable functions; all dis- 
tributions with compact support and aU square-integrable functions can be viewed as tempered distributions. 
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correlated functions^. Within the hierarchy {w"} the truncated functions can be computed 
recursively: 

w^ixi, X2) = w'^{xi, X2) - {xi)w^ {X2) (1-32) 



The asymptotic property Eq. (1.30) of the is converted to the simpler property 

«;?(xi,...x„)^0, (1.33) 

if any space like separation Xi — xj goes to infinity. 

Let us now consider a hierarchy of functions P'^{ki, ...kn) (for simplicity totally sym- 
metric under permutation of their arguments). We define a generating functional V{f} 
for an arbitrary function f{k): 

^^/> ^T.^J ^"(^1' ■■■kn)f(k,)...f(kn). (1.34) 

The P" can be found from 

Now one can simply use the definition Eq. (1.31) to obtain a relation between the hi- 
erarchy {P"} and the hierarchy of the truncated functions {Pt} in terms of respective 
generating function, namely 

V{f} = e^^^f\ (1.36) 

This relation is called the Linked cluster theorem [18]. Having in mind the relations 
Eqs. (1.34) and(1.36), one observes that the CCM parametrization introduced inEqs. (1.1) 
and (1.3) is a natural reparametrization of the CIM defined in Eq. (1.2) which incorporates 
the size-extensivity [5], at finite levels of truncation. This leads to a cluster decomposition 
and linked decomposition property of CCM at any level of truncation. Thus for extensive 
variables as the energy, the linked terms lead to contributions which obey the proper linear 
scaling in the particle number. The CIM contains unlinked diagrams for the ground-state 
energy expectation value, and thereby suffers from the so-called size-extensivity prob- 



^For free fields all trancated functions with n ^ 2 vanish, thus it suffices to know the two-point func- 
tions. 



12CHAPTER 1. BASIC STRUCTURE OF THE COUPLED-CLUSTER FORMALISM 



lems. 

The cluster decomposition condition Eq. (1.30) is fundamental to a quantum field 
theory. It may break down partly when, in a statistical-mechanical sense, the theory is 
in a mixed phase [19]. This implies that there is more than one possible vacuum state, 
and therefore the cluster decomposition should be restored in principle if one builds the 
Hilbert space on one of the vacua. 

Notice that in QCD, the cluster decomposition for colour singlet objects can break 
down due to the confinement. In axiomatic local quantum field theory with an indefinite 
metric space u (e.g., Minkowski space) the following theorem holds; for the vacuum 
expectation values of two (smeared local) operators A and B with spacelike distance R 
from each other [20, 21]: 

\{Q\A{x)B{0)\Q) - {Q\A{x)\Q){Q\B{0)\Q)\ 

^ r const X i?-3/2-i-2Arg-MH ^^lerQ is a mass gap M, 
~ \ const X if there is no mass gap, 

where M is the mass gap. (Herein, we assume = 0, however, a positive N is possi- 
ble for the indefinite inner product structure in u.) In order to avoid the decomposition 
property for product of unobservable operators A and B which together with Kugo-Ojima 
criterion^ [20, 22] for the confinement is equivalent to failure of the decomposition prop- 
erty for coloured clusters, there can not be a mass gap in the indefinite space u. This 
would thus eliminate the possibility of scattering a physical state into color singlet states 
consisting of widely separated colored clusters (the "behind-the-moon" problem). How- 
ever, this has no information for the physical spectrum of the mass operator in Ti. which 
indeed does have a mass gap. In another words, there is a mass gap (but not in the full 
indefinite space u) in the semi-definite physical subspace induced by the BRST charge 
operator Qb, t'phy = KerQs, where states within are annihilated hy Qb, 

t'phys = {\ip) e V : Qslip) = 0} = KerQfi, 
n{QB,v) = KerQs/Imgs, (1.38) 

where states in ImQ^ = (KerQs)^ are called BRST-coboundaries in the terminology of 
de Rham cohomology and do not contribute in z/phys [23]. It has been shown [24] that there 
is a connection between Kugo-Ojima criterion and an infrared enhancement of the ghost 
propagator. In Landau gauge, the gluon-ghost vertex function offers a convenient pos- 



*The Kugo-Ojima confinement scenario describes a mechanism by which the physical state space con- 
tains only colourless states, the colored states are not BRS-singlets and therefore do not appear in ^-matrix 
elements, since they are confined [20, 22]. 
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sibility to define a non-perturbative running coupling. The infrared fixed point obtained 
from this running coupling determines the two-point color-octet interactions and leads 
to the existence of unphysical massless states which are necessary to escape the cluster 
decomposition of colored clusters [20]. Notice that a dynamical mechanism, responsible 
for breaking of cluster decomposition for colored objects, is yet to be discovered, since 
our knowledge about color confinement is still preliminary. 
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Chapter 2 

Hamiltonian renormalization Groups 



2.1 Introduction 

The main goal of traditional renomialization theory is to determine when and how the 
cancellation of divergences originating from the locality of quantum field theory may 
occur. This is essential if one wants to have meaningful quantitative results. However, it is 
by no means obvious how quantum fluctuations associated with short distance scales can 
be incorporated and controlled through the choice of only a few parameters, typically the 
bare masses and coupling constants, or by the counterterms in renormalized perturbation 
theory. 

The development of Wilson's renormalization group (RG) formalism [25] allowed 
physicists to produce a logically consistent picture of renormalization in which perturba- 
tion theory at arbitrary high energy scale can be matched with the perturbation expansion 
at another scale, without invoking the details of intermediate scales. In the Wilsonian 
approach, all of the parameters of a renormalizable field theory can be thought of as 
scale-dependent objects and their flows are governed by the so-called RG differential 
equations. 

Motivated by Wilson's picture, effective field theory (EFT) approaches have been in- 
troduced to replace complicated fundamental theories with simpler theories based only 
on the relevant degrees of freedom at the physical scale of interest [26]. The basis of 
the EFT concept is the recognition of the importance of different typical energy scales in 
nature, where each scale has its own characteristic degrees of freedom. In strong inter- 
actions the transition from the fundamental to the effective level is induced by a phase 
transition that takes place around Aqcd — IGeV via the spontaneous breaking of chiral 
symmetry, which generates pseudoscalar Goldstone bosons. This coincides, of course, 
with the emergence of nuclei and nuclear mather, as opposed to the quark-gluon plasma 
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and quark matter expected to occur at high temperature and high density. Therefore, at 
low energies (E < Aqcd), the relevant degrees of freedom are not quarks and gluons, but 
pseudoscalar mesons and other hadrons. The resulting description is a chiral EFT, which 
has much in common with traditional potential models. In particular, there might be an 
intermediate regime where a non-relativistic model is inadequate but where relatively few 
hadronic degrees of freedom can be used to faithfully describe both nuclear structure and 
response. If this is the case, one should be able to describe strongly interacting hadronic 
systems in this effective model. 

The power of Hamiltonian methods is well known from the study of non-relativistic 
many-body systems and from strongly-interacting few particle systems, even though a La- 
grangian approach is usually chosen for relativistic theories. One might prefer to obtain 
the effective interaction using the covariant Lagrangian formalism and then consider the 
ground state and collectively excited states in the non-covariant Hamiltonian formalism 
by exploiting many -body techniques. For the description of physical states and in partic- 
ular bound states, the Hamiltonian formalism is preferable over the Lagrangian one. This 
is due to the fact that such problems are not naturally defined covariantly. For a bound 
state the interaction time scale is infinite, and thus a time-independent approach is better 
suited. However, the Lagrangian can not generally be converted to a Hamiltonian if the 
effective Lagrangian contains higher-order time-derivatives, since no Legendre transfor- 
mation exits for such a case. Therefore one may wish to obtain the effective interaction 
in a unified self-consistent way within the Hamiltonian formalism. The renormalization 
group transformation is used to derive the physical Hamiltonian that can describe exper- 
iment. To implement a RG calculation, one should define a space of Hamiltonians and 
find a certain RG transformation which maps this space into itself. Then one should study 
the topology of the Hamiltonian space, by searching for the fixed points and studying the 
trajectories of the Hamiltonian with respect to these fixed points. At the fixed points we 
have a scale-invariant quantum field theory. Near the fixed point, the irrelevant operators 
in the original Hamiltonian have small coefficients, and we are only left with the relevant 
and marginal operators. The coefficients of these operators correspond to the parameters 
of the renormalized field theory. In this way, in a renormalizable field theory, we disre- 
gard information about the evolution of irrelevant terms. The fact that irrelevant terms 
can be dropped near the fixed point does not necessarily imply that they are unimportant 
at low-energy scales of experimental interest. This depends on how sensitive the physical 
observables are to physics near the scale of the cutoff. Therefore, it is of interest to de- 
velop a RG method in which the irrelevant variables are treated on an equal footing with 
the relevant and marginal variables (the schemes presented here have this advantage). 

Hamiltonian methods for strongly-interacting systems are intrinsically non-perturbative 
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and usually contain a Tamm-Dancoff type approximation, in the sense that one expands 
the bound state in states containing a small number of particles. This truncation of the 
Fock space gives rise to a new class of non-perturbative divergences, since the truncation 
does not allow us to take into account all diagrams for any given order in perturbation the- 
ory. Therefore renormalization issues have to be considered carefully. Two very different 
remedies for this issue are the use of light-front Tamm-Dancoff field theory (LFFT) [27] 
and the application of the coupled cluster method (CCM) [5, 28]. However, both meth- 
ods are too complicated to attack the issue in a self -consistent way. In the last decade 
extensive attempts have been made to give a workable prescription for renormalization 
within the Hamiltonian formalism [29, 30, 31, 32]. Commonly unitary transformations 
are used to decouple the high- and low-energy modes aiming at the partial diagonalization 
of the Hamiltonian. One of the most elegant approaches in this context is the similarity 
renormalization group (SRG) proposed by Glazek and Wilson [29] (and by Wegner [30] 
independently). The SRG [29, 30] is designed to be free of small energy denominators 
and to eliminate interactions which change the energies of the unperturbed states by a 
large amount. However, there are several problems with this approach: it is hard to in- 
corporate loop expansions within the method, the SRG can not systematically remove 
interactions which change the number of particles (i.e, when the Hamiltonian is not diag- 
onal in particle number space), and most importantly, the computations are complex and 
there is no efficient non-perturbative calculation scheme. 

In this chapter we introduce a new method [32, 33] for obtaining the low-energy ef- 
fective operators in the framework of a CCM approach. The transformation constructed 
avoids the small denominators that plague old-fashioned perturbation theory. Neither 
perturbation theory nor unitarity of the transformation are essential for this method. The 
method is non-perturbative, since there is no expansion in the coupling constant; nonethe- 
less, the CCM can be conceived as a topological expansion in the number of correlated 
excitations. We show that introducing a double similarity transformation using linked- 
cluster amplitudes will simplify the partial diagonalization underlying renormalization 
in Hamiltonian approaches. However, a price must be paid: due to the truncation the 
similarity transformations are not unitary, and accordingly the hermiticity of the resul- 
tant effective Hamiltonian is not manifest. This is related to the fact that we have a 
biorthogonal representation of the given many-body problem. There is a long tradition 
of such approaches. The first we are aware of are Dyson-type bosonization schemes [34]. 
[Here one chooses to map the generators of a Lie algebra, such that the raising generators 
have a particularly simple representation.] The space of states is mapped onto a larger 
space where the physically realizable states are obtained by constrained dynamics. This 
is closely related to CCM formalism, where the extended phase space is a complex man- 
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ifold, the physical subspace constraint function is of second class and the physical shell 
itself is a Kahler manifold [35]. The second is the Suzuki-Lee method in the nuclear 
many -body (NMT) problem [36, 37], which reduces the full many-body problem to a 
problem in a small configuration space and introduces a related effective interaction. The 
effective interaction is naturally understood as the result of certain transformations which 
decouple a model space from its complement. As is well known in the theory of effective 
interactions, unitarity of the transformation used for decoupling or diagonalization is not 
necessary. Actually, the advantage of a non-unitarity approach is that it can give a very 
simple description for both diagonalization and ground state. This has been discussed by 
many authors [38] and, although it might lead to a non-hermitian effective Hamiltonian, 
it has been shown that hermiticity can be recovered [35, 39]. Notice that defining a good 
model space can in principle control the accuracy of CCM [40, 41]. 

To solve the relativistic bound state problem one needs to systematically and simul- 
taneously decouple 1) the high-energy from low-energy modes and 2) the many- from 
the few-particle states. We emphasize in this chapter that CCM can in principle be an 
adequate method to attack both these requirements. Our hope is to fully utilize Wilso- 
nian Exact renormalization group [25] within the CCM formalism. Here the high energy 
modes will be integrated out leading to a modified low-energy Hamiltonian in an effective 
many-body space. Notice that our formulation does not depend on the form of dynamics 
and can be used for any quantization scheme, e.g., equal time or light-cone. 



2.2 Traditional approaches and their problems 

The Tamm-Dancoff approximation [42] was developed in the 1950's to describe a rela- 
tivistic bound state in terms of a small number of particles. It was soon revealed that the 
Tamm-Dancoff truncation gives rise to a new class of non-perturbative divergences, since 
the truncation does not allow us to take into account all diagrams at a given order in per- 
turbation theory. On the other hand, any naive renormalization violates Poincare symme- 
try and the cluster decomposition property (the cluster decomposition means that if two 
subsystems at very large space-like separation cease to interact then the wave function 
becomes multiplicatively separable). One of the simplest example of the Tamm-Dancoff 
approximation is the constituent picture of QCD where one describes a QCD bound state 
within a truncated Fock-space, 



IV') = 01 k?) + <l>2\qq9) + (pslqqqq) + ■ 
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we use a shorthand notation for the Fock-space where g is a quark, q an antiquark, and 
g stands for a gluon. Now the bound-state problem is solved via the Schrodinger equa- 
tion H'^^^\ip) = E\ip). It is well known that in the complicated equal-time vacuum 
bound states contain an infinite number of particles that are part of the physical vacuum 
on which hadrons are built. On the other hand, interactions in a field theory couple states 
with arbitrarily large difference in both free energy and numbers of particles. Thus any 
Fock-space expansion can hardly be justified without being supplemented with a pre- 
scription for decoupling of the high-energy from the low-energy modes and the many- 
body from few-particle states. In the context of the Hamiltonian formulation this problem 
can be expresed by asking how the Hamiltonian matrix is diagonalized in particle- and 
momentum-space. 

In his earliest work Wilson[43] exploited a Bloch type transformation [44] to reduce 
the Hamiltonian matrix by lowering a cutoff which was initially imposed on the individual 
states. Later Wilson abandoned this formulation in favour of a Lagrangian one. The most 
important reason was that the Bloch transformation is ill-defined and produces unphysical 
divergences. These divergences emerge from denominators which contain a small energy 
difference between states retained and states removed by the transformation, and appear 
across the boundary line at A . 

Two remedies for this issue are the use of light front coordinates [27] and application 
of the CCM [28]. In the light-front Tamm-Dancoff field theory (LFFT) the quantization 
plane is chosen to coincide with the light front, therefore the divergences that plagued the 
original theories seem to disappear [45] since here vacuum remains trivial. Furthermore, 
not having to include interactions in boost operators allows a renormalizable truncation 
scheme [46]. One of the most important difficulties in LFFT is the complicated structure 
of the renormalization process [47]. In principle, ad-hoc counterterms can not be pre- 
vented if one is to preserve the underlying symmetry. In the standard form of CCM, on 
the other hand, the amplitudes obey a system of coupled non-linear equations which con- 
tain some ill-defined terms because of ultraviolet divergences. It has been shown [48] that 
the ill-defined amplitudes, which are also called critical topologies, can be systematically 
removed, by exploiting the linked-cluster property of the ground state. This can be done 
by introducing a mapping which transfers them into a finite representations without mak- 
ing any approximation such as a coupling expansion. Thus far this resummation method 
has been restricted to superrenormalizable theories due to its complexity. 

Recently Wilson and Glazek [29] and independently Wegner [30] have re-investigated 
this issue and introduced a new scheme, the similarity renormalization group (SRG). The 
SRG resembles the original Wilsonian renormalization group formulation [43], since a 
transformation that explicitly runs the cutoff is developed. However, here one runs a 
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Figure 2. 1 : Two ways to run a cutoff on free energy. In top a cutoff on the magnitude 
of the energy is lowered from A to A, leading to a small matrix (this scheme corresponds 
to the original Wilson RG). In down we show similarity renormalization group scheme 
where a cutoff is imposed on the off-diagonal part of the matrix and diagonalization is led 
to a narrow matrix with a band width A. CTa denotes the counterterms. 
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cutoff on energy differences rather than on individual states. In the SRG framework one 
has to calculate a narrow matrix instead of a small matrix, and the cutoff can be conceived 
as a band width (see Fig. [2.1]). Therefore by construction the perturbation expansion for 
transformed Hamiltonians contains no small-energy denominators. Here we review the 
general formulation of the SRG. 



2.2.1 Glazek- Wilson formulation 

The detailed description of Glazek- Wilson RG method can be found in Ref. [29]. Here we 
only concentrate on the key elements of their method. We introduce a unitary transforma- 
tion aiming at partially diagonalizing the Hamiltonian so that no couplings between states 
with energy differences larger than A are present. The unitary transformation defines a set 
of Hamiltonians Hx which interpolate between the initial Hamiltonian (A = A ) and the 
effective Hamiltonian H^. This A plays the role of a flow parameter. We will assume that 
Hx is dominated by its diagonal part which we will denote Hqx with eigenvalues and 
eigenstates 

{f\Hox\i)^Eix{f\i). (2.1) 

Therefore the Hamiltonian at scale A can be written Hx — Hqx + Hjx, where Hjx is 
non-diagonal part. Notice that Hqx is not necessarily the bare free Hamiltonian which is 
independent of A. In order to introduce the infinitesimal unitary transformations which 
induce the infinitesimal changes in Hx when A changes by an infinitesimal amount, we 
need to define various zones of the operators We introduce an auxiliary function Xijx = 
{i\xx\j) of the states with labels i and j for a given A. If we denote the eigenvalue of Hqx 
with Eix, X is defined as 

Eix + Ejx + A 

The modulus of the function Xijx is close to 1 when one of the energies is much larger than 
the other and also much larger than the cutoff A and x approaches when the energies are 
similar or small in comparison to the cutoff. One then introduces smooth projectors 

Uijx = {i\ux\j) = u{xijx), (2.3) 

and 

rijx = 1 - Uijx = r{xijx). (2.4) 

By means of u and r one can separate each matrix into two parts M — D{M) -\- R{M), 
where we define 

D{M)ij = UijxMij, (2.5) 
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and corresponding 

R{M)ij = VijxMij. (2.6) 

The functions u and r are needed in order to ensure smoothness and differentiabiUty, and 
consequently a continuous transition between different parts of the Hamiltonian matrix. 
They are the key elements in making the diagonalization free of small-energy denomina- 
tors. These functions are implemented as a "form factor" in every vertex of the interaction. 
We construct an infinitesimal unitary transformation eliminating the part of Hamiltonian 
which has only non-zero elements far away from the diagonal, thus they can not produce 
small-energy denominators. This continuous unitary transformation satisfies, 

^ = [T,,H,]. (2.7) 
The generator Tx is anti-Hermitian and is chosen in such a way that 

= D{Qx) = Qxux, (2.8) 

where Q is arbitrary in the far-off-diagonal region. In terms of Q, the matrix elements of 
Eq. (2.7) satisfy 

-^^QijA + Uijx = Tijx{Ejx - Eix) + \Tx, Hjxlij. (2.9) 

Notice that Qx can be arbitrary in the far-off diagonal region where ux = ^ = if it is 
finite and its derivative is finite. In the above equation we have two unknowns, and 
Tx. As additional input we use the fact that, for a given A, the Hamiltonian Hx can be 
additionally unitary transformed without violating the relation D[Qx) = Hx- We regroup 
the above equation in the form, 

Uijx — Tijx{Ejx — Eix) — [Tx, Hix]ij -^^Qijx — Gijx- (2.10) 

The unknowns on the left hand-side are determined by first solving Eq. (2.10) neglecting 
the commutator on the right-hand side; the result is then substituted into the right-hand 
side, and solved iteratively until convergence is reached. The D and R parts of the opera- 
tor G can be defined, 

Uijx ^^^^ = D{Gx)ij, 
Tijx{Ejx — Eix) = —R{Gx)ij- (2.11) 
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By evaluating the matrix elements of both sides of the above equations in different zones 
of the operators, one obtains differential equations for matrix elements of Qx. Therefore, 
from Eqs. (2.10, 2.11) one can immediately obtain the generator T and the Hamiltonian 
flow equation in terms of the matrix elements {f\Hx\i) = Hxfi and {f\Tx\i) = Txfi, 

TijX = '-E-x (j^^'^^^^'^ d\^^'^^ ' ^^'^^^ 

— Uijx[J-x,Rix\ij + rijx — — Rijx- (2.13) 

It is obvious from the above equations that no small-energy denominators Eix — Ejx arises 
out-side of the band (off-diagonal region), and within the band zone, we have T — 0. The 
equations (2.12,2.13) can be solved iteratively. 



2.2.2 Wegner Formulation 

Wegner's formulation of the SRG is defined in a very elegant way aiming at diagonaliza- 
tion of the Hamiltonian in a block-diagonal form with the number of particles conserved 
in each block. Again, a unitary transformation is used with flow parameter s that range 
from to oo, 

^^[T{s),H{s)]. (2.14) 

We separate the Hamiltonian in a diagonal part H^, and the remainder H^. We use the fact 
that tr H"^ is invariant under the unitary transformation, therefore we have 

tr + tr = tr = const. (2.15) 

This means that tr falls monotonically if tr Hj increases. One can use Eq. (2.14) to 
obtain 

~~ds^ ^ 2^ ^ ^ ^ "^'^ ^ (^ij-^ji ~ Hij'^ji) = 2 ^ TjiHji{Hjj - Ha). 

i i j ij 

(2.16) 

In order to ensure that falls monotonically, one can simply choose the generator 

as Tji = Hji{Hjj - Hii) or 

T{s)^[Ha{s),Hr{s)]. (2.17) 
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One can show by making use of this definition that. 



dHij{s) 



ds 



j 



d_ 
ds 



(2.18) 



k ij 



Because Ylii^j Hfj falls monotonously and is restricted from below. When s — > oo the 
derivative vanishes and we have — > 0, at this limit the procedure of block-diagonalization 
is completed. At this point, the unitary transformation Eqs (2.14,2.17) is completely de- 
fined. The only freedom is in the choice of separation of the Hamiltonian into a "diagonal" 
and a "rest" part. Of course this depends on the given physical problem. As a illustration 
of the method, we show how perturbation theory can be applied in this formalism. For a 
given values of s we have. 



where the superscript denotes the order in the bare coupling constant. The part is the 
free Hamiltonian. The index r denotes the rest of the Hamiltonian. Note that generally 
the diagonal part in the flow equation is the full particle number conserving part of the 
effective Hamiltonian. The choice of only H^'' as the diagonal part leads to the simplest 
band-diagonal structure where the particle number is conserved. We use the basis of the 
eigenfunctions of the free Hamiltonian H^^\i) — Ei\i) to obtain the matrix elements of 
Eqs (2.14,2.17), 



T{s) = ms), Hr{s)] = [H^f\ //(^)] + [H^f\ HP] + ... . 



(2.19) 



— J'(l) _|_ 2^(2) _|_ ^ 



(2.20) 



dHii 



-{E, - E.fH^ + [T('\Hilj] -{E,- E.fH^. + (2.21) 
{Ei - Ej)Hiij + {E, - Ej)Hl^j + (2.22) 



ds 



where the energy differences are given by 



ni n2 

Ei — Ej = ^ Eik — ^ Ejk, 

k=l k=l 



(2.23) 
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and Eik and Ejk are the energies of the creation and annihilation particles, respectively. 
To leading order in perturbation theory, one finds 



ds 

Uij{s) = e-(^*-^^)'^ (2.24) 

The flow-parameter s has dimension 1/ (energy)^ and is related to the similarity width 
A (ultraviolet cutoff) by s = This implies that matrix elements of the inter- 

action which change the number of particles are strongly suppressed, since we have 
\Ei — Ej\ > X. The similarity generator in Wegner's formulation corresponds to the 
choice of a gaussian similarity function with uniform width. In the next leading order, 
one has to deal separately with the diagonal and rest parts. In analogy to the Glazek- 
Wilson method we introduce H^"^^ — u(s)Hr'^\s), and the solution reads, 

hS]{s) = H^^fjis = 0) + ds'u{s')[T('\H^\,j{s'), 

Jo 

H^M = = 0) + f ds'[T^'\H^%,As'). (2.25) 

It is obvious that for the non-diagonal term a smooth form factor appears in a natural way 
to suppress the off-diagonal interaction. In other words, the particle number changing in- 
teractions are eliminated while a new terms are produced. The procedure can be extended 
to arbitrarily high orders. The counterterms can be determined order-by-order using the 
idea of coupling coherence, namely that under similarity transformation Hamiltonian re- 
mains form invariant (see next section). 



2.3 Coupling coherence condition 

One of the most severe problems for the traditional light-front RG is that an infinite num- 
ber of relevant and marginal operators are required [49]. This is due to the fact that 
light-front cutoff violates the underlying symmetry, e.g., Lorentz invariance and gauge 
symmetries. Since these are continuous symmetries, their violation in principle leads to 
infinite number of symmetry violating counterterms (with new couplings), in order to 
maintain the symmetry of the effective Hamiltonian. In terms of the effective field theory 
approach, some sort of fine tuning is required to fix the strength of the new couplings 
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so that the underlying symmetry is restored. One should note that in order to reduce the 
number of momentum degrees of freedom, one must introduce a real cutoff, such as a 
momentum cutoff or a lattice cutoff. However, dealing with divergences do not require 
necessarily a decrease in degrees of freedom. In fact, one may even increase the degrees 
of freedom, e.g., the Pauli-Villars or the dimensional regularization methods [50]. The 
main idea behind the coupling coherence renormalization condition [51, 52] is that the 
Hamiltonian is form-invariant on the RG trajectory. This condition isolates and repairs 
the hidden symmetries [51, 52]. 

H{A)=n{fx), (2.26) 

In order words, rewriting the Hamiltonian in different degrees of freedom does not change 
the operator itself. One may think of Hdi) as QCD written in terms of constituent quarks 
and gluons and H{A) as the same QCD Hamiltonian written in terms of canonical quarks 
and gluons, associated with partons and current quarks. The SRG and the coupled-cluster 
RG with coupling coherence allows one to construct effective theories with the same 
number of couplings as the underlying fundamental theory, even when the cutoff violate 
symmetries of the theory. This does not preclude the emergence of the new couplings, 
however, they depend on the original coupling and will vanish if the fundamental cou- 
plings are turned off. This boundary condition together with the RG equations determines 
their dependence on the fundamental coupling. 



As an illustrative example [52], we consider following interaction 

V{<f>)-^<Pt+^<Pl + ^<Pl4, (2.27) 

where 0i and 02 are scalar fields. We want to investigate that under what conditions the 
couplings are independent of each other. Consider the Gell-Mann-Low equations [53] up 
to one-loop in perturbation theory, ignoring the masses, 

dt 3 



7TeA^ + eA3(Ai + A2), (2.28) 



where t — log(A//x) and ^ = /i/(167r^). Assume that there is only one independent 
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coupling A = Ai, and A2, A3 are functions of A. Now one can simplify Eq. (2.28), 

V 12 7 dX ^ 12 ^' 

(SA^ + ^A^)^ = ^A^ + AA3 + A2A3. (2.29) 

In the leading order, the above equations have two distinct solutions, one is when A2 = A 
and A3 = 2A. In this case we have 

n0) = ^(0? + 02)'- (2.30) 

Therefore we recover the 0(2) symmetric theory. The other solution is A2 = A and 
A3 = 6A which leads to two decoupled scalar fields, 

^(0) = ^ ((</>! + <l>2)' + (01 - h)') , (2.31) 

One can conclude that A2 and A3 do not run independently with the cutoff if and only 
if there is a symmetry which connects their strength to Ai. The condition that a limited 
number of couplings run with cutoff independently reveals the symmetries broken by the 
regulator and repairs them. More interesting, it may be used as well to uncover symme- 
tries that are broken by the vacuum. This may reconcile the trivial vacuum in light-front 
field theory and vacuum symmetry breaking problem. 



2.4 General formulation of the similarity renormaliza- 
tion group 

In this section we pave the way for an introduction of the coupled-cluster RG, and consider 
the similarity renormalization group in a more general framework without requiring the 
unitarity. The discussion in this section is partially based on the work of Suzuki and 
Okamoto [37]. Let us consider a system described by a Hamiltonian H{A) which has, 
at the very beginning, a large cut-off A. We assume that the renormalized Hamiltonian 
H^^{A) up to scale A can be written as the sum of the canonical Hamiltonian and a 
"counterterm" He (A.), 

H^^{A) = H{A) + Hc{A) . (2.32) 
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Our aim is to construct the renormalized Hamiltonian by obtaining this counterterm. Now 
imagine that we restrict the Hamiltonian to a lower energy scale (fi), where we want to 
find an effective Hamiltonian H^^{ii) which has the same energy spectrum as the original 
Hamiltonian in the smaller space. The cut-off /i can be conceived as a flow parameter. 
The value of ^ = A corresponds to the initial bare regulated Hamiltonian. Formally, 
we wish to transform the Hamiltonian to a new basis, where the medium-energy modes 
fjL < k < A, decouple from the low-energy ones, while the low-energy spectrum re- 
mains unchanged. We split the Hilbert space by means of flow-parameter n into two 
subspaces, the intermediate-energy space Q containing modes with fx < k < A and a 
low-energy space P with k < fi. Our renormalization approach is based on decoupling of 
the complement space Q from the model space P. Thereby the decoupling transformation 
generates a new effective interaction 6H{ii, A) containing the effects of physics between 
the scales A and fi. One can then determine the counterterm by requiring coupling co- 
herence [31, 51, 52], namely that the transformed Hamiltonian has the same form as the 
original one but with A replaced by // everywhere. (This is in contrast to the popular 
Effective Field Theory approach, where one includes all permissible couplings of a given 
order and fixes them by requiring observable computed be both cutoff-independent and 
Lorentz co variant.) 

We define two projector operators, also called P and Q, which project a state onto the 
model space and its complement, satisfy P^ = P,Q'^ = Q, PQ = and P + Q = 1. We 
introduce an isometry operator G which maps states in the P- onto the Q- space. 



The operator G is the basic ingredient in a family of "integrating-out operators", and 
passes information about the correlations of the high energy modes to the low-energy 
space. The operator G obeys G = QGP, GQ = 0, PG = and G" = for n ^ 
2. The counterintuitive choice that G maps from model to complement space is due to 
the definition Eq. (2.36) below (c.f. the relation between the active and passive view of 
rotations). In order to give a general form of the effective low-energy Hamiltonian, we 
define another operator X{n, ji, A), 



\q) = G\p) 



(|g) e g, \p) e P) . 



(2.33) 



X{n, At, A) = (1 + G){1 + G^G + GG^f . 



(2.34) 



(n is a real number.) The inverse of X{n^ A) can be obtained explicitly. 



X-^{n, A) = (1 + G^G + GGt)-"(l - G) . 



(2.35) 
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The special case n = is equivalent to the transformation introduced in ref. [54] to relate 
the hermitian and non-hermitian effective operators in the energy-independent Suzuki- 
Lee approach. We now consider the transformation of H{A) defined as 

H{n, /X, A) = X-\n, /x, A)H{A)X{n, /x, A) , (2.36) 

where we have 

H{A) H{n, 11, A) = H{iJi) + 5H{iJi, A) . (2.37) 
One can prove that if H{n, ji, A) satisfies the desirable decoupling property, 

QH{n,ii,A)P ^0, (2.38) 

or more explicitly, by substituting the definition of X{n,n,A) and X~^{n, n, A) from 
Eqs. (2.34)-(2.35), 

QH{A)P + QH{A)QG - GPH{A)P - GPH{A)QG = , (2.39) 

that H^^{ij) = /x) = PH{n, ji, A)P is an effective Hamiltonian for the low energy 
degrees of freedom. In other words, it should have the same low-energy eigenvalues as 
the original Hamiltonian. The proof is as follows: 

Consider an eigenvalue equation in the P space for a state \(t>{k)) e P, 

PH{n, /X, A)P\(j){k)) = EkPX-\n, /x, A)X(n, /x, A)P|0(A;)) . (2.40) 

By multiplying both sides by X{n, /x, A) and making use of the decoupling property Eq. 
(2.38), we obtain 

H{A)X{n, /X, A)P\<P{k)) = EkX{n, /x, A)P\<P{k)) . (2.41) 

This equation means that Ek in Eq. (2.40) agrees with one of the eigenvalue of H{A) and 
X{n, ji, A)P\(f){k)) is the corresponding eigenstate. Now we demand that 

H'\fx) = H{fj,) + Hc{fi). (2.42) 

This requirement uniquely determines the counterterm He. 

In the same way we can also obtain the Q-space effective Hamiltonian, from the defi- 
nition of H{n, ji, A). It can be seen that if G satisfies the requirement in Eq. (2.39), then 
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we have additional decoupling condition 

PH{n,i^,A)Q = 0. (2.43) 

Although the above condition is not independent of the decoupling condition previously 
introduced in Eq. (2.38) in the exact form, but this is not the case after involvement of 
some approximation (truncation). We will argue later that both of the decoupling condi- 
tions Eqs. (2.38) and (2.43) are necessary in order to have a sector-independent renormal- 
ization scheme. The word "sector" here means the given truncated Fock space. Let us now 
clarify the meaning of this concept. To maintain the generality of the previous discussion, 
we use here the well known Bloch-Feshbach formalism [44, 55, 56]. The Bloch-Feshbach 
method exploits projection operators in the Hilbert space in order to determine effective 
operators in some restricted model space. This technique seems to be more universal than 
Wilson's renormalization formulated in a Lagrangian framework. This is due to the fact 
that in the Bloch-Feshbach formalism, other irrelevant degrees of freedom (such as high 
angular momentum, spin degrees of freedom, number of particles, etc.) can be systemat- 
ically eliminated in the same fashion. 



Assume that the full space Schrodinger equation is H\iIj) = E^\%1)) and for simplicity 
has been normalized to one. The similarity transformed Schrodinger equation now 
reads H\'il)x-^) = E^\il)x-^), where we defined ItZ-'x-i) = ""^1^) and H h & similarity 
transformed Hamiltonian. This equation is now separated into two coupled equations for 
P- and Q-space. 

{E^ - PHP)P\'il;x-^) = PHQliPx-^), (2.44) 
{E^ - QHQ)Q\tlJx-^) = QHP\^x-^). (2.45) 

We may formally solve Eq. (2.45) as 

and substitute this equation into Eq. (2.44) to obtain a formally exact uncoupled equation 
in P- space, 

H^''P\^x-^)=E^P\ijx-^), (2.47) 

where we have 

H^^ = PHP + PHQ ^-^^QHP. (2.48) 

E^ — QHQ 
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The effective Hamiltonian H^^ constructed in this fashion is explicitly energy dependent. 
This equation resembles the one for Brueckner's reaction matrix (or "G"-matrix) equation 
in nuclear many-body theory (NMT). In the same way for arbitrary operator O (after a 
potential similarity transformation), we construct the effective operator. Let us define a 
similarity transformed operator O, 

^ = + ^)^(^ + Q)l^a:-l)(^a.-l|, (2.49) 

where and |^) are eigen functions of Hamiltonian and we have — 
Yl = 1 ^^d P + Q = 1. We now write (^^-ilQ and Q|^/'x-i) in terms of 

their solution in the P-space obtained in Eq. (2.46). One can immediately show 

(0,-i|pog|v^,-i) - {cj>,-^\P0Q J^"^ p\^,-.), 
((/),-i|gop|v^,-i) = (0,-i|p ^^^ gopiv^.-i), 
(0,-i|gog|v^,-i) = (0.-i|p ^ goQ ^^^ p|v^.-i). (2.50) 

By plugging the above equations into Eq. (2.49), one can obtain the effective operator in 
the P-space 

^=EEi'^--^)('^--^i^<^'''^i^--^)(^--^i = (2.51) 

where we have 

O^ff ^ POP + PHQ ^-^^QOP + POQ ^-^^QHP 

— QHQ — QHQ 

+ PHQ ^-^^QOQ ^-^^QHP. 

E^ — QHQ E^ — QHQ 

(2.52) 



Notice that Eq. (2.52) can be converted into the form of Eq. (2.48) with the effective 
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Schrodinger equation Eq. (2.47) when O ^ H^. 

The E'-dependence in Eqs. (2.48) and (2.52) emerges from the fact that the effective 
interaction in the reduced space is not assumed to be decoupled from the excluded space. 
However, by using the decoupling conditions introduced in Eqs. (2.38) and (2.43), we 
observe that energy dependence can be removed, and the effective operators become 

= PHP ^n{n,n) , 
O^ff = POP^O{n,^). (2.55) 

The decoupling property makes the operators in one sector independent of the other sec- 
tor. The effects of the excluded sector is taken into account by imposing the decoupling 
conditions. This is closely related to the folded diagram method in NMT for removing 
energy-dependence [57]. (It is well-known in NMT that ^'-dependence in the G-matrix 
emerges from non-folded diagrams which can be systematically eliminated using the ef- 
fective interaction approach). The above argument was given without assuming an explicit 
form for X and thus the decoupling conditions are more fundamental than the prescription 
used to derive these conditions. We will show later that one can choose a transformation 
X, together with the model space and its complement, which avoids the occurrence of 
"small energy denominators". We now show that Lorentz covariance in a given sector 
does not hinge on a special form of similarity transformation. Assume the existence of 
ten Poincare generators Lj satisfying 

[L,,Lj]=Y,4l,, (2.56) 
where are the structure coefficients. One can show that if the operators Li satisfy the 



'To this end, we organize Eq. (2.52) when O ^ H, 

1 1 
^E^-QHQ ^ E^-QHQ E^ - QHQ^^' ^ E^ - QHQ^ 



QcB ^ POP + PHq( + ^^^ + h^^QHQ h^^)QHP, 

\E^- OHO E,,. - OHO E^ - OHO K,. - OHO J 



= POP + PHQ ( ^-^^ + =—) QHP. (2.53) 

^\E^-QHQ {E^-QHQ){E^-QHQ)J^ 

One now can use Eq. (2.46) to rewrite the last part of the above equation into the form of 

where we used Q -I- -P = 1 andQ^ = Q. Having made used of Eqs. (2.51,2.53,2.54), one can immediately 
obtain Eq. (2.47) with the effective Hamiltonian in the form of Eq. (2.48). 
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decoupling conditions QLiP — and PLiQ — Oit follows that 

[Lf,Lf] = J24Lf. (2.57) 

This leads to a relativistic description even after simultaneously integrating out the high- 
frequency modes and reducing the number of particles. Indeed we conjecture that requir- 
ing the decoupling conditions makes the effective Hamiltonian free of Lorentz-noninvariant 
operators for a given truncated sector regardless of the regularization scheme. However, 
one may still be faced with an effective Hamiltonian which violates gauge invariance (for 
e.g., when sharp cutoff is employed). 

Note that the solution to Eq. (2.39) is independent of the number n. One can make use 
of Eq. (2.39) and its complex conjugate to show that for any real number n, the following 
relation for the effective low-energy Hamiltonian 

n{n,fi)^n\-n-l,n). (2.58) 

The case n = — 1/2 is special since the effective Hamiltonian is hermitian, 

n{-l/2,ii) = {P + G''Gy/^H{A){P + G){P + G^G)-^/'^ . (2.59) 

Hermiticity can be verified from the relation [58] 

e^P = (l + G')(P + G'^G')-^/^ (2.60) 

where, 

T = arctan(G' - G^) = ^ ^ — L- [GiG'^GY - h.c.) . (2.61) 

n=0 ^ 

Since the operator T is anti-hermitian, is a unitary operator. From the above expression 
Eq. (2.59) can be written in the explicitly hermitian form 

n (^-^, i)j - Pe-^H{K)e'P . (2.62) 

As was already emphasized, renormalization based on unitary transformations is more 
complicated and non-economical. Thus we will explore a non-unitary approach. An 
interesting non-hermitian effective low-energy Hamiltonian can be obtained for n = 0, 



n{0,f^,A)^PH{A){P + QG) 



(2.63) 
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This forai resembles the Bloch and Horowitz type of effective Hamiltonian as used in 
NMT [55], and was the one used by Wilson in his original work on quantum field the- 
ory [44](see section II). In the context of the CCM, this form leads to the folded diagram 
expansion well known in many-body theory [16]. It is of interest that various effective 
low-energy Hamiltonians can be constructed according to Eq. (2.55) by the use of the 
mapping operator G which all obey the decoupling property Eq. (2.38) and Eq. (2.43). 
Neither perturbation theory nor hermiticity is essential for this large class of effective 
Hamiltonians. 



2.5 The coupled-cluster renormalization group 

The CCM approach is, of course, just one of the ways of describing the relevant spec- 
trum by means of non-unitary transformations. According to our prescription the model 
space is P : {|-^) <8) |0, b)h, L < //}, where |0, b)h is a bare high energy vacuum (the 
ground state of high-momentum of free-Hamiltonian) which is annihilated by all the high 
frequency annihilation operators {C/} (for a given quantization scheme, e.g., equal time 
or light-cone) , the set of indices {/} therefore defines a subsystem, or cluster, within 
the full system of a given configuration. The actual choice depends upon the particular 
system under consideration. In general, the operators {Cj} will be products (or sums of 
products) of suitable single-particle operators. We assume that the annihilation and its 
corresponding creation {C]} subalgebras and the state |0, 6)/^ are cyclic, so that the linear 
combination of state {Cj|0, b)h} and { ^(6, 0|C/} span the decoupled Hilbert space of 
the high-momentum modes, where ji < H < K. It is also convenient, but not 

necessary, to impose the orthogonality condition, (0|C/Cj|0) = J), where 8{I, J) is 
a Kronecker delta symbol. The complement space is Q : {\L) ® {\H) — \Q,h)h)}- Our 
main goal is to decouple the P-space from the Q-space. This gives sense to the partial di- 
agonalization of the high-energy part of the Hamiltonian. The states in full Hilbert space 
are constructed by adding correlated clusters of high-energy modes onto the P-space, or 
equivalently integrating out the high-energy modes from the Hamiltonian, 



I/) 
(/I 



X{^x, A)|p) = e^e-^'|0, 6)^ (g) |L) = e^|0, 6)^ (g) |L) , 
(L| ^(6,0|X-i(//, A) = (L| ft(0|e^'e-^, 



(2.64) 
(2.65) 
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Figure 2.2: Shows the Wilsonian shells and the low-energy phase space {s™, sf'} which 
is induced by integrating out the high-energy modes {C/, C]}, please see the text for 
details. 

where the operators X{ijl, A) and A) have been expanded in terms of independent 

coupled cluster excitations /, 

I 

Y^sTCi. (2.66) 
I 

Here the primed sum means that at least one fast particle is created or destroyed (/ ^ 
0), and momentum conservation in P © Q is included in s/ and s'j. Sm{S'^ are not 
generally commutable in the low-energy Fock space, whereas they are by construction 
commutable in the high-energy Fock space. One of the crucial difference between our 
approach and the traditional CCM [3] is that, here sf and s'p are only c-number in the 
high-energy Fock space, but they are operators in the low-energy Fock space (unless 
when ^ — > 0). We have in the Wilsonian high-energy shell [5*, C\\ — [S' , C/] = 



m=0 



m=0 
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(see Fig. 2.2). Therefore, one can guarantee the proper size-extensivity and consequently 
conformity with the linked-cluster theorem at any level of approximation in the Wilsonian 
high-energy shell. Nevertheless one can still apply the standard CCM to the induced 
effective low-energy Hamiltonian to extend the proper size-extensivity to the entire Fock 
space. 

It is immediately clear that states in the interacting Hilbert space are normalized, 
(/I/) = h{b,0\0,b)h = 1. We have two types of parameters in this procedure: One 
is the coupling constant of the theory (A), and the other is the ratio of cutoffs (/x/A). The 
explicit power counting makes the degree of divergence of each order smaller than the 
previous one. According to our logic, Eq. (2.63) can be written as 

n{fi) = h{b,0\X-\fx,A)H{A)X{fx,A)\0,b)h, (2.67) 

with X{ij,, A) and X~^{ij,, A) defined in Eq. (2.64) and Eq. (2.65). We require that effec- 
tive Hamiltonian H{ijl) obtained in this way remains form invariant or coherent [31, 51, 
52]. This requirement satisfies on an infinitely long renormalization group trajectory and 
thus does constitutes a renormalized Hamiltonian. Thereby one can readily identify the 
counter terms produced from expansion of Eq. (2.67). 

The individual amplitudes for a given m, {sf', s'^} = {^/) s^}^, have to be fixed by 
the dynamics of quantum system. This is a complicated set of requirements. However, 
we require less than that. Suppose that after a similar transformation of Hamiltonian, H, 
we obtain an effective Hamiltonian of the form 

H = if (low) + i/free(high) + ClVijCj , (2.68) 

where Vjj is an arbitrary operator in the low frequency space. The I and J indices should 
be chosen such that the last term in Eq. (2.68) contains at least one creation- operator and 
one annihilation- operator of high frequency. By using Rayleigh-Schrodinger perturbation 
theory, it can be shown that the free high-energy vacuum state of iJfree(high) is annihi- 
lated by Eq. (2.68) and remains without correction at any order of perturbation theory. 
Having said that, we will now consider how to find the individual amplitudes {s/, s'j}m 
that transfer the Hamiltonian into the form Eq. (2.68). We split the Hamiltonian in five 
parts: 

H^H, + H^'^ihigh) + Vc{Cj) + Va{Ci) + Vb , (2.69) 

where Hi contains only the low frequency modes with k < fi, H2 is the free Hamiltonian 
for all modes with n < k < A, Vc contains low frequency operators and products of 
the high frequency creation operators C| and Va is the hermitian conjugate of Vc- The 
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remaining terms are contained in Vb, these terms contain at least one annihilation and 
creation operators of the high energy modes. Our goal is to eliminate Vc and Va since Vb 
annihilates the vacuum. The ket-state coefficients {s/}m are worked out via the ket-state 
Schrodinger equation H{A)\f) = E\f) written in the form 

(0 1 Cie-^He^ 1 0) = , V7 ^ . (2.70) 

The bra-state coefficients {sj, s'j}m are obtained by making use of the Schrodinger equa- 
tion defined for the bra-state, {f\H{A) = {f\E. First we project both sides on C] 1 0), then 
we eliminate E by making use of the ket-state equation projection with the state (0|e'^'c] 
to yield the equations 

(0|e^'e-^[//,Cj]eV^'|0) = 0, V7 ^ . (2.71) 

Alternatively one can in a unified way apply e^e~^' Cl\0) on the Schrodinger equation for 
the bra-state and obtain 

(0|e^'e-^//e^e-^'cJ|0) = , V7 ^ . (2.72) 

Equation (2.70) and Eqs. (2.71) or (2.72) provide two sets of formally exact, microscopic, 
operatorial coupled non-linear equations for the ket and bra. One can solve the coupled 
equations in Eq. (2.70) to work out {sjjm and then use them as an input in Eqs. (2.71) 
or (2.72). 

It is important to notice that Eqs. (2.70) and (2.71) can be also derived by requir- 
ing that the effective low-energy Hamiltonian defined in Eq. (2.67), be stationary (i.e. 
5'H{iJi) — 0) with respect to all variations in each of the independent functional {s/, s'j}^. 
One can easily verify that the requirements 5'H{ij)/5sf — and SHdj^/Ss'p — yield 
Eqs. (2.71) and (2.70). The combination of Eqs. (2.70) and (2.71) does not manifestly 
satisfy the decoupling property as set out in Eqs. (2.38) and (2.43). On the other hand 
Eqs. (2.70) and (2.72) satisfy these conditions. Equations (2.70) and (2.72) imply that all 
interactions including creation and annihilation of fast particles ("J") are eliminated from 
the transformed Hamiltonian H{fi) in Eq. (2.67). In other words, these are decoupling 
conditions leading to the elimination of Vc and Va from Eq. (2.69), which is, in essence, 
a block-diagonalization. Therefore it makes sense for our purpose to use Eqs. (2.70) 
and (2.72) for obtaining the unknown coefficients. We postpone the discussion of the 
connection between decoupling and variational equations in next section. 

So far everything has been introduced rigorously without invoking any approxima- 
tion. In practice one needs to truncate both sets of coefficients {sj, s'j}jn at a given order 
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of m. A consistent truncation scheme is the so-called SUB(A/^, A^) scheme, where the 
n-body partition of the operator {S*, S"} is truncated so that one sets the higher partition 
with I > J\f to zero up to a given accuracy m = Ai. Notice that, Eqs. (2.71) and (2.72) 
provide two equivalent sets of equations in the exact form, however after the truncation 
they can in principle be different. Eqs. (2.70) and (2.72) are compatible with the decou- 
pling property at any level of the truncation, whereas Eqs. (2.70) combined with (2.7 1) are 
fully consistent with HFT at any level of truncation. Thus the low-energy effective form 
of an arbitrary operator can be computed according to Eq. (2.55) in the same truncation 
scheme used for the renormalization of the Hamiltonian. In particular, we will show that 
only in the lowest order (m = 0), equations (2.71) and (2.72) are equivalent, independent 
of the physical system and the truncation scheme. 

Although our method is non-perturbative, perturbation theory can be recovered from 
it. In this way, its simple structure for loop expansion will be obvious and we will observe 
that at lower order hermiticity is preserved. Now we illustrate how this is realizable in our 
approach. Assume that Vc and Va are of order A, we will diagonalize the Hamiltonian, at 
leading order in A up to the desired accuracy in /x/A. We use the commutator-expansion 
Eq. (1.18) to organize Eq. (2.70) perturbatively in order of m, aiming at elimination of the 
high momenta degree of freedom up to the first order in the coupling constant, thus yields 

m = 0: {0\Ci{Vc+[H2,So]m = 0, 

m = 1 : {0\Ci{[H^, So] + [H2, S,] + [Va, ^1] + [Vc, Si])\0) - , 

m = n : (0|C,([//i, Sn-i] + [H^, Sn] + [Va, Sn] + [Vc, ^„])|0) = , (2.73) 

where / 7^ 0. Notice that 5*0 is chosen to cancel Vc in the effective Hamiltonian, hence 
it is at least of order of A, consequently it generates a new term [Hi , So] which is of 
higher order in ///A and can be canceled out on the next orders by iSi. The logic for 
obtaining the equations above is based on the fact that Sn should be smaller than Sn-i (for 
sake of convergence) and that the equations should be consistent with each other. Since 
H2,Va,Vc ~ A and Hi ji, from Eq. (2.73) we have the desired relation Sn ~ x'^n-i- 
Notice that if ~ A then one needs to keep all coupled equations in Eq. (2.73) up to 
m = n. In other words, there is no perturbative expansion in ratio of cutoffs. This may 
introduce small-energy denominator i. e, 1/(A — /i), since one may obtain the solution 
of these equations, e. g.. So in form of a geometrical series which can be resummed and 
leads to a small-energy denominator ^. Having said that as far as the renormalization is 



■^Thanks to Prof. F. Wegner, for bringing this point to my attention. 
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concern we are always interested on condition that A >> j2. Furthermore, for ~ A one 
should resort to the non-perturbative decoupling equations which are by construction free 
of small-energy denominator. The same procedure can be applied for Eq. (2.72) which 
leads to the introduction of a new series of equations in order of m, 

m = 0: {0\iVA-[H2,S',])C\\0)^0, 

m = 1 : (0|([//i, S'o] + S[] + [Vc, S[] + [Va, S[] - [Va, S,])C\\0) = , 

m = n: (0|([i/i, S'^_,] + [H^, S'J + [Vc, S'J + [Va, S'J - [Va, Sn])Cj\0) = 0(2.74) 

Alternatively, we can use Eq. (2.71) to yield the equations 

m = : {0\{[Va, C]] - [[H^, C]],S',] ) |0) = , 

m = l: {0\{[[Va,C\],S,] - [[H2,C\],S[] - [[Va,CI], S[])\0) ^ , 



m — n 



{0\{[VA,Cj],S^] - [[H„Cj],S',] - [[Va,CI],SX\0)^0. (2.75) 



It is obvious that at order m = 0, Eqs. (2.14) and (2.75) are the same and S'q = Sq, which 
indicates that the similarity transformation at this level remains unitary. It should be 
noted that diagonalization at first order in the coupling constant introduces a low-energy 
effective Hamiltonian in Eq. (2.67) which is valid up to the order A^. In the same way, 
diagonalization at second order in A modifies the Hamiltonian at order A^ and leads gen- 
erally to a non-unitarity transformation. In this way one can proceed to diagonalize the 
Hamiltonian at a given order in A with desired accuracy in ///A . Finally, the renormal- 
ization process is completed by introducing the correct Z{A) factors which redefine the 
divergences emerging from Eq. (2.67). 



2.6 The decoupling conditions versus the variational prin- 
ciple 

In section 2.3 we introduced the effective low-energy phase space {sj, s'jjm induced by 
integrating out the high-energy modes. We argued later that the induced low-energy phase 
space can be obtained either from Eqs. (2.70, 2.72) which are manifestly compatible with 
the decoupling conditions or from Eqs. (2.70, 2.71) compatible with the variational prin- 
ciple and the Hellmann-Feynman theorem. Here we will show that in fact this two for- 
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mulations are related through an isomorphic transformation (or in a strict sense by a sym- 
plectomorphism). In what follows, for simplicity, we absorb the factor )™ in Eq. (2.66) 
into the definition of Sm and S'„^, hence we assume SmiS'^) — (f )™, therefore we have 
desirable convergence relation S„i+i{Sl^_^_i) ~ ^SmiS'^). In the spirit of Schrodinger 
representation in quantum field theory, one may now define a generalized many-body ket 
and bra wave functional at a given scale 

1^^^) = eS(M,A)g-5'(M,A)|o) ^ e^('^'^)lO), (^^1 = (0|e^'('^'^)e-^('^'^\ (2.76) 

It is clear that by construction, we have {ip^\ij^) = 1 at any level of approximation. Notice 
that the bra and ket are parametrised independently, since they are not hermitian-adjoint 
of each other. Therefore, we have a biorthogonal representation of the many-body system. 
We will illuminate below the main underlying reasons behind this parametrization. With 
this definition the effective low-energy Hamiltonian Eq. (2.67) is rewritten as, 

n{fi) = (^^1^14)- (2.77) 

Here, the bra and ket wave functional are built by adding independent clusters of high- 
energy correlation (where < A; < A) in the vacuum of the free high-energy Hamiltonian, 
or by integrating out the high-energy modes from the Hamiltonian. The unknown low- 
energy operators {s™, s'J"'} are obtained by solving the Schrodinger equation in the high- 
energy shell (iJ, < k < A) and are decoupling equations, 

Di = QHP = {0\Cie-'^He^\0) = , V/ (2.78) 
D2 = P^g = (0|e^'e-%eV^'C||0) = 0, V/ ^ . (2.79) 

The decoupling equations Eqs. (2.78,2.79) imply that all interactions containing creation 
and annihilation of "/" high-momentum particles are eliminated from the transformed 
Hamiltonian while it generates new low-momentum interactions through {sY, s'f}. These 
equations show the changes of the generalized wave functional and accordingly the effec- 
tive Hamiltonian with the flow-parameter Alternatively one can obtain {s^*, s'p} by 
requiring that the effective low-energy be stationary SHdi) — with respect to all varia- 
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tions in each of the independent functional {sf-, s'J^}, 

^2M. ^ ^ {Q\e^'e-^[H, C\]e^e-^'\Q) = , V7 ^ (2.80) 

br/(n.'\ 

^ Di = (2.81) 



where Di is introduced in Eq. (2.78). Notice that we have already shown that Eq. (2.80) 
is as well derivable from dynamics . The decoupling conditions Eqs. (2.78, 2.79) seems 
generally to be in conflict with variational equations (2.80, 2.81). Here we will show 
that in fact this two formulations are related via an isomorphic transformation (or in a 
strict sense by a symplectomorphism). Let us introduce a new set of variables {cr'j' , cr'j'} 
within the low-energy phase space, by transforming the induced low-energy phase space 
operators set {s'f-, s'p} into a new set {a'^^a^}. 



(2.82) 



'I — ) 

where the low-energy operator ujji is defined as 

u;ji = {0\Cie^'cl\0). (2.83) 
In the same way, one may conversely write s in terms of a, 

STi^:^)= Y^^7^Ji, (2.84) 

J 

where we have used the following orthogonality property, 

£'j/ = (0|C,e-^'cJ|0), 

y^^LbjKU^KI = y^^^JKl^KI = S{J,I). (2.85) 
K K 

We use the canonical transformation Eq. (2.82) to rewrite the decoupling equations Di 
and D2 in terms of derivative with respect to new variables, one can immediately show, 

= (0|Ge^'e-%(A)e*e-'|0) = |M . |^ + E'*r..||. (2.86) 

I I J J 
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In the same fashion, having made use of the closure identity / in the high-energy Fock 
space Q, and Eq. (2.85), one can rewrite the other decoupling equation D2 in the form of 

D2 = (0|e^'e-%(A)e^(/)e-^^'cl|0), 

= ?Y(/.)(0|e-^'C||0)+ 5^'(S,^{(0|e^'e-^[i/(A),Cj]e^|0), 

J 

+ (0|e^'C]e-^'(/)e^'e-%(A)eV^'|0)}, 

- + 1^ ^^^^ + 2^ ^^J-^kLki, (2.87) 

^ J J,K 

where in final step we have made use of the orthogonality relation Eq. (2.85) and the 
relation in Eq. (2.86). The operator Lu is defined as 

Uj = J2{0\e''cl,e-''cl\0){0\CKe-''Cj\0). (2.88) 

K 

The operator L/j is symmetric and doubly linked. It is obvious from Eqs. (2.86, 2.87) that 
the requirement of STi.{fi)/dd-J' = and di-i^ji) / Sa^ = lead directly to the decoupling 
conditions Di = and D2 = 0. However the reverse is not always correct. Therefore, 
the decoupling conditions are weaker than the variational equations. A set of similar 
canonical variables as in Eq. (2.82) (where the variables are c-numbers), was introduced 
by Arponen, Bishop and Pajanne [3] in the context of the traditional CCM and turned out 
to be quite practical. 

We have already proved (in section 2.4) that the effective low-energy operators Eq. (2.67) 
or Eq. (2.77) equipped with decoupling conditions Eqs. (2.78,2.79) have the same low- 
energy eigenvalues as the original Hamiltonian. The decoupling conditions are thus suffi- 
cient requirements to ensure partial diagonalization of the Hamiltonian in the particle and 
momentum space. 

Having obtained the effective bra and ket Eq. (2.76) at a given scale ji, one can com- 
pute the effective low-energy of an arbitrary operator A, 

A{sirs'i) = {hM'>P,^)^ (2-89) 

One may now pose the question if the approximation (truncation) used to obtain the in- 
duced bra- and ket-state, and accordingly the effective Hamiltonian, is sufficient to obtain 
the effective low-energy of a given operator by Eq. (2.89). On the other hand, one might 
be curious about the necessity of the double similarity transformation and as well the in- 
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troduction of two independent sets of variables {s^, s'™}, since one could have defined a 
single transformation H e^^He^ where S is defined in Eq. (2.66) (without introduc- 
ing a new set of variables {s'p}, and invoking the definition of the bar state). In this case, 
the variables set {s™} would be determined by Eq. (2.78) alone. Let us pursue this idea 
and apply the Hellmann-Feynman theorem^ to such a parametrization. In our framework, 
the Hellmann-Feynman theorem states that if we perturb the Hamiltonian H ^ H' — 
H + J A, where J is an infinitesimally small (a source term) and A is an arbitrary oper- 
ator, such that the effective Hamiltonian changes as — > 7Y' = + JdH/dJ + 0( J^) 
then we have, 

dn/dJ = i = {ilj^\dH/dJ\iP^), (2.90) 

where we define an effective low-energy operator A — {ip^\A\ipij). We can find the 
effective low-energy operator A by using the Hellmann-Feynman theorem, 

rl - - . . - - i^s™ 

A = —{0\e-^{H + JA)e^\0) = (0|e-W|0) + (0|e-^i/e^C]|0)-f , (2.91) 
dJ 6 J 

where we used [S, C\] = 0. If we now calculate 5s™ from the first decoupling equation 
(2.78) which involves only S, having retained only 0{J, 5S) terms, we find 

(0|C/e-(^+^^)(^+'^^)e^^'^^^^|0> = {0\Cie-^JAe^\0) + {0\Cie-^[H,C\]e^\0)5s'y = 0. 

(2.92) 

Therefore, one can show that, 

Cjio)^ ^ Q(7Y - Qe-^He^Q)-^Qe-'^Ae^\0), (2.93) 

where the operator Q = 1 — |0)(0| = ^ Cj\0){0\Cj is introduced. Now one can 

J 

make use of Eq. (2.93) to show that the right-hand side of Eq. (2.91) can be recast in the 
following form, 

i= (0|e^'e-^^e^|0), (2.94) 
where we introduced the notation e^' by, 

(0|e^' = (0| + {0\e-^He^Q{n - Qe'^ He^ Q)'^ Q. (2.95) 

Interestingly, Eq. (2.95) satisfies the second decoupling condition Eq. (2.79). Thereby 



^The Hellmaim-Feymnaii theorem is originally introduced for the ground state, however its proof is 
more general and can be applied here. 
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Eq. (2.94) becomes exactly equivalent to Eq. (2.89). Therefore the use of the the Hellmann- 
Feynman theorem and assuming the coupled-cluster parametrization leads naturally to the 
definition of the bra-state in Eq. (2.76) and emergence of a new set of variable S'. In 
other words, no other bra-state parametrization (including hermitian- adjoint of ket-state) 
is compatible with the Hellmann-Feynman theorem. 

It is well known that the traditional multiplicative renormalization is not sufficient for 
the renormalization of more than one composite operator inserted into the renormalized 
Green functions. To avoid ad hoc subtractions, one can introduce the composite operators 
into the Lagrangian with a space-time dependent source (coupling). It has been shown 
that the renormalization of the source produces requires counter terms to render all Green 
functions containing the insertion of composite operators renormalized [59]. These new 
counterterms do not affect the renormalization of the original theory. 

The main advantage of the compatibility of our parametrization with the Hellmann- 
Feynman theorem underlies that here the effective low-energy of an arbitrary operator, 
can be obtained in the same truncation scheme SUB(A/', M) used for the Hamiltonian 
matrix ^. Moreover, the above-mentioned technique used in Lagrangian formalism, can be 
employed in our framework by means of the Hellmann-Feynman theorem. Therefore, the 
renormalization of an arbitrary operator can be calculated in a unified way. This is indeed, 
one of advantages of the non-unitary parametrization of the similarity renormalization 
group. 



2.7 The symplectic structure 

Despite extensive progress in development of various RG techniques, little is still known 
about the geometrical interpretation of the RG [60] . Here, we introduce the geometrical 
structure emerging from our approach. We define a low-energy action-like functional A, 
having integrated out fast modes and making use of the reparametrization Eq. (2.82), 

= / dt{i^ ai&i -n{iJ,,a,a)), 

di(-i J]'^jaj-7Y(//,(7,a)), (2.96) 



I 



'^As Thouless pointed out, the Hellmann-Feynman theorem immediately imphes that an expectation 
value {A) of an arbitrary operator is computed diagrammatically from the same set of Goldstone diagrams 
as for the energy (H), where all interaction is replaced by the operator A. 
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where in the final step we employed integration by parts. The operator TC{fi, a, a) is the 
effective low-energy Hamiltonian defined in Eq. (2.67) and ai{d-i) are defined, 

M M 
m=0 m=0 

The stationary of A with respect to the complete set of variables {ai,aj} for a given 
truncation SUB(j\^, M) yields 



(2.98) 



These equations can be obtained for the set {s/, s'j} as well, without invoking the canon- 
ical transformation Eq. (2.82). However, if one wants to obtain the above equations in 
a straightforward manner from time-dependent Schrodinger equation, then the canonical 
transformation Eq. (2.82) is necessary. 

We note that although our transformation is not unitary and the parametrization of S 
and S' have been introduced in a very asymmetric fashion, their fundamental dynamics 
in the low-energy phase space follows the canonical equation of motion. Therefore the 
induced low-energy operators ai and a/ are canonically conjugate in the terminology of 
the classical Hamiltonian mechanics. This can be made even more suggestive by defining 
a new set of variables, the generalized field 0/ and their canonically conjugate generalized 
momentum densities tt/. 

In terms of the new operators the equation of motion Eq. (2.98) can be rewritten into the 
form 

^Jl = ii„HW} = -'-m, (2.100) 

at S(f)i 

where a generalized Poisson bracket {A, B} for two arbitrary operators is defined as 
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For the nonzero mutual Poisson brackets of canonical coordinates and momentums we 
have {0/, TTj} = 6{I, J). We can also look at the behaviour of the low-energy effective 
operator for a product of operators. One can write the product of operators after a double 
similarity transformation as a product of transformed operators, 

{■^^\AB\iff,) = {0\e^'e-^ABe^e-^'\0) = {0\AB\0), (2.102) 

= AB + {0\e''e-'Ae^e-''Cj\0){0\Cje^'e-'Be^e-''\0), 
I 

where we have used the closure identity in Q-space. The effective operator A and B are 
defined by Eq. (2.89). Now we express the last part of Eq. (2. 102) in terms of a derivatives 
of the effective operators A{B) with respect to low-energy phase space {ai, aj}. In fact, 
we have already shown such relations in Eqs. (2.86, 2.87). The same relationships are 
valid for any arbitrary operator for new variables {a/, a/}. Therefore with a replacement 
of if — > A, S and H ^ A,Bin Eqs. (2.86, 2.87) we have the desired equations. The same 
procedure can be carried out for evaluating of (t/S^ISAI^^). After some straightforward 
algebra, one finds that, 

4^\AB - BA\,p^) = (4l[^,i^]|V^M) = HA-B}, (2.103) 

where the Poisson bracket is as defined in Eq. (2. 101). 

Therefore our parametrization is clearly suggestive that the effective low-energy phase 
space induced by integrating out the high-energy particles can be described in terms 
of highly non-linear classical dynamics with the canonical coordinates {^i,t^i), with- 
out losing any quantum-field-theoretical information. The canonical equation of motions 
Eq. (2.100) are valid at any level of approximation SUB(AA, M) on the renormalization 
group trajectory. In this way, we are led to a manifestation of the correspondence prin- 
ciple in a more generalized form, that is the induced low-energy Hamiltonian obtained 
by integrating short-distance modes are governed by a hierarchy of non-linear classical 
mechanical equations for quasi-local fields (j)i and tt/. One may hope to recover the full 
classical mechanics (with c-numbers variables as coordinates) at the other extreme of the 
RG trajectory. 

The general form for the action completely determines the symplectic structure of 
our low-energy phase space (0/,7r/). The phase space is the cotangent bundle^ of the 
configuration space C,T — T*{C) (the coordinates 0/ label the points of configuration 



^In differential geometry, the cotangent bundle is the union of all cotangent spaces of a manifold [23]. 
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space C) [23]. The Poisson bracket Eq. (2.101) induces a symplectic structure 

CO ^^cuabdx'' Adx^, (2.104) 

where the coordinates on the symplectic manifold F are denoted by e {4>i\ ^iid '^a6 is 
the inverse of the symplectic matrix 

= {x\ x^} ^ {A, B} = (2.105) 

Moreover, a symplectic form defines an isomorphism between the tangent and cotangent 
spaces of F. One may associate a vector field Xj to every function / e C°°(F) by 

i{X^)u = -d/, (2.106) 

where %{X) and denote the interior product and the symplectic gradient of /, respec- 
tively. X/ is the so-called the Hamiltonian vector field of /, and it generates a flow on F 
which leaves uj invariant, since the Lie derivative of uj along X/ is zero. In this way, one 
can rewrite the Poisson bracket in the form 

{/,^} = i{Xf)i{X,)u = i{X^)dg = u{X^,X,) e C°°(r), (2.107) 

which shows the change of g along Xf. We require that uj to be closed {duj = 0)^, which 
implies the Jacobi identities {duj){Xf, Xg, X^) = 0. Therefore, the existence of Poisson 
bracket in our definition of the coordinates of the phase-space introduces a symplectic 
manifold for the phase-space. 



2.8 The constrained induced phase space 

In this section we closely follow Ref. [35] (one should pay some extra care here, since the 
phase space here is operatorial rather than c-number). Let us assume that the renormalized 
form of an arbitrary operator A, Eq. (2.89) can be obtained by, 

^^{0|eVWe-^'|0).i^jgM0>. ^^^^^^ 

symplectic manifold (F, ui) is a smooth real A'^-dimensional manifold without boundary, equipped 
with a closed non-degenerate two-form w, i. e., duj = where d is the exterior differential [23]. 
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Generally the left-hand side does not agree with right-hand side (it is by no means clear 
that this relation will be held after a truncation). In order to ensure the unitarity of the 
similarity transformation we require, 

(Oe^ = ^ - — , = ■ (2.109) 

^ ' (Ole^+e^lO) ' ' (0|e^^e^|0) 

Thereby, the ket-state and bra-state defined in Eq. (2.76) become hermitian- adjoint of 
each another. We assume the induced low-energy phase space to be a complex mani- 
fold {{sj, s'j), {s*j, s'j)}. The hermiticity conditions are introduced by constraint functions 

Xiifi) and x*M 

M = ((0|e^^e^|0))-^(0|e^^C]e^|0) - {0\S'Cj\0), 

= ((0|e^'e^|0))-^(0|e^'C,e^|0) - (0|C,^'*|0). (2.110) 

Therefore the physical submanifold shell is defined through: 

M 



Xj^O^ PX-\i,, A)C\X{i,, A)P = J2 ^ s'j, 

m=0 
M 

X*j^O^ PX-\pi,A)CiX{pi,A)P = ^s';- = (2.111) 



m=0 



where X(ijl, A) is defined in Eq. (2.64, 2.65). This implies that in the physical subman- 
ifold where we have exact hermiticity, and there is an isomorphic mapping between a 
cluster of high-energy creation C\ (annihilation Cj) operators and a low-energy operators 
s'j{s'j). This isomorphism is invariant under the renormalization group transformation. 
We introduce a Poisson bracket for the complex representation of the phase space, 

' ^~'Y\Ssi5s'j Ss'jSsi^ Ss'j* 5s} SsjSs'j*)' ^' ^ 

The non-zero commutators of the canonical coordinates follow the canonical symplectic 
structure {s/, s'j} = {s}, s'}} = S{I, J). The nature of the constraint can be revealed by 
considering the commutators between the constraints functional %/ and x}. After some 
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tedious but straightforward algebra, one obtains, 

{XhXj} = {x*i,X*j} = 0, 

{Xi,Xj} = 2 ((0|e^'e^C,Cye-^'|0) - ^j^}) . 

This impUes that the constraints are of second class and do not correspond to any gauge 
symmetry degrees of freedom. However, these superfluous degrees of freedom can be 
eliminated by Dirac bracket technique [61]. Of course, since the constrained manifold 
N has dimension less than the full manifold M, one may define a puUback map /* : 

T*{M) T*{N) to obtain the induced symplectic structure of the constraint surface ci>°, 
a)° = f*uj. In analogy to ordinary CCM [35], we define a symplectic two-form in the full 
manifold, 

uj = Y^{dsi A ds'j + ds'; A ds}). (2.113) 
I 

The induced symplectic two-form on the physical shell can be found by substituting the 
value of ds'j and ds'j in terms of dsj and ds} by using on-shell condition dxi — dx} — 0, 
hence we find 

= 2 w^'^dai A da*j, (2.114) 
i,j 

where w^'^ = consequently is closed and a) is positive matrix. Therefore the in- 
duced physical low-energy phase space is a Kahler manifold ^[23]. This defines a positive 
hermitian metric in the physical shell ^. 

In the following next two sections, we apply our RG formalism to compute the effec- 
tive Hamiltonian for cf)'^ and extended Lee theory up to two- and one-loop order, respec- 
tively. 



2.9 Example I: theory 

In this section, we obtain the effective Hamiltonian of 0^ theory up to two-loop order in 
equal-time quantization. In following we will quote from Ref. [33]. The bare 0^ theory 



Kahler manifold is a Hermitian manifold (M,g) whose Kahler form CI is closed: dfl = 0. The metric 
g is called the Kahler metric of M [23]. 

^Notice that it is well known that the geometrical quantization can be applied on Kahler manifold since 
it has a natural polarization. 
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Hamiltonian is [50] 

j d?x Q7r2(x) + - + m^)(i){x) + . (2.115) 

According to our logic the ultraviolet-finite Hamiltonian is obtained by introducing coun- 
terterms, which depend on the UV cutoff A and some arbitrary renormalization scale. This 
redefines the parameters of the theory and defines the effective low-energy Hamiltonian. 
The renormalized Hamiltonian has the form 

H = j d^x {^t:\x) + ^ v^0(x) ( - + Z^m^) v^(/)(x) + Z,Zlgct>\x) + ...^ 

(2.116) 

Each of the Z-factors has an expansion of the form. 

Z = 1 + /i(A)A + /2(A)A2 + ..., (2.117) 

where A is a generic coupling constant of theory and has been defined at a given renor- 
malization scale M. The functions will be obtained order-by-order, by summing up 
contributions of the fast modes between ji and A, in the sense that Z{K) — > Z{iJi) and 
/„(A) — > fn{lA- This means that the low-energy correlation functions are invariants of 
the renormalization group flow. One can therefore assume that the Z's are initially 1 and 
choose the corresponding /'s from the condition that the cut-off dependence be cancelled 
out after computing the effective Hamiltonian in the desired loop order. Even though the 
newly generated interactions are sensitive to the regularization scheme (as is well known 
[62], a sharp cutoff may lead to new non-local interaction terms), nevertheless one can 
ignore these if they are finite and do not produce any divergence as the cutoff A ap- 
proaches to infinity. We now split field operators into high- and low-momentum modes; 
(f){x) — (pLix) + (t>H{x), where ^^(a:) denotes modes of low-frequency with momentum 
k < n and (pnix) denotes modes of high-frequency with momentum constrained to a 
shell /i < A; < A. The field ^^(a;) can be conceived as a background to which the ^//(x)- 
modes are coupled. Therefore, in the standard diagrammatic language, integrating out the 
high-frequency modes (pnix) implies that only high-frequency modes appear in internal 
lines. The field (Ph{x) is represented in Fock space as 

M^)- E -^{ake''^ + 4e-''n, (2.118) 

IJ.<k<A ^ 
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where Uk = Vk^ + rn^ and the operators and a\ satisfy the standard boson commu- 
tation rules. From now on all summations are implicitly over the high-frequency modes 
fi < k < A. The Hamiltonian in terms of high- and low-frequency modes can be written 
as, after normal ordering with respect to high-frequency modes, 

H ^Hi + H2 + Vb + Vc + Va, (2.119) 



where we define. 



^i(j>+q-k)x ^i{k-p)x 



2uJr 



Va = 

^—i(k+p+q)x ^ ^-i{k+p+q+r)x 

The high-energy configurations in the Fock space are specified by {C/ — > ni=i ^fei} 
and {Cj — > ni=i ^L)}- Up to two-loop expansion, our renormalization scheme requires 
to keep Sis') at least to order n = 4, which allows us to eliminate the pure terms Vc 
and Va at a lower level of expansion. The <S'(S") operators consistent with a SUB{A, m) 
truncation scheme are. 



= J Yl {^^ + Si alal + ajajaj + S^ alalalal^ , 
^'m ^ J Y (^^m + S'^ akap + S'^ ahttpttq + S^ akapttgUr^ . (2.121) 
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We split the diagonalization of the Hamiltonian matrix in an upper and lower triangle 
part, by using the double similarity transformation. One may notice that the "most non- 
diagonal" terms in the Hamiltonian are Vc and Va (in the light-front Hamiltonian such 
terms do not exist because modes with longitudinal momentum identically zero are not 
allowed). The potential Vb is already partially diagonalized and does not change the 
vacuum of the high-energy states. Therefore, here we employ a minimal scheme, aiming 
at removal of Va and Vc only. 



We restrict ourselves to the elimination of the high-energy degrees of freedom up to 
the first order in the coupling constant g and second order in the ratio of cutoffs /x/A. 
Therefore, our truncation scheme is called SUB {A, 2). For m — one finds. 



^0^ 



r-3 



-9—, 



-9 



v^ 



Si = -9 



St = -9 



\/2 



V^ 



■,(2.122) 



where the V^~^ are defined in Eq. (2. 120). Here, one has S'^ — Sl. At this stage the results 
for the one-loop renormalization can be computed. We evaluate the effective Hamiltonian 
by substituting <S'(<S") from Eqs. (2.121) and (2.122) into Eq. (2.67). In order to achieve 
renormalization, one should identify the potentially divergent terms ( when A oo) in 
the expansion of H^^{ii). Such a process generally can be done by inventing a power- 
counting rule, using the property Sn ~ ^>S'„_i. Here we take Uk ^ \k\ for 3> m and 
replace by J (|^- The standard tadpole one-loop mass renormalization arises from 
Vb due to normal-ordering. We add this divergent term to Hi and renormalize the bare 
mass 



^^l-loop 



— 



{0\Vb\0) 
27r2 



6g 



2ujk 47r2 



(A^ 



(fx(f)'^{x), 



(A^ - 1.' 



(2.123) 



In this order the contribution of the terms [Vc, S], [Va, S'] and [Hi, S{S')] are zero, after 
projection on to the high-energy vacuum. The only divergent contributions come from 
[yf^^\ 'S'q^^^] due to a double and third contraction of the high-frequency fields respec- 
tively. There are two other divergent terms, {[Vq^^\ S'^^^^j, however they are harmless 
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and are cancelled out by the divergence of [[H2, Sq], Sq^^^]. One thus obtains 

18^ / <P\x)(j)\y) ,ik+p)(x-v) 



(27r)9 } 



^ I it}\ .1 . \ tt)\ If \ , , , , , , 

^(k+v+DK^-y) _ (2.124) 



In general evaluation of integrals like Eq. (2.124) may produce non-localities. This is due 
to the fact that the total momentum in integrands of Eq. (2.124), namely ri = p + q and 
r2 = A; + p + g are in the low-momentum space. To evaluate such integrations, one can 
firstly reduce the potential divergent integrals by a change of variable, for example for the 
first integrand we use p, ? — > p, ri, and then expand the integrand in ri/p. Therefore, after 
expansion and evaluating the momentum integrals, one may be faced with non-analytic 
terms in the low-momentum space. However here these are non-divergent and will thus 
be ignored. We find 



^^i-ioop ^ -gl„(^^^|d3:,0^(:,)-g(21n2-l)A2y(i3x<i 



+ TZ-^^^i-] / ci^a;(V0(a;))2 + finite terms. (2.125) 



One can immediately deduce the renormalization factors Zg and from above expres- 
sion 



Zg = l + ;^ln(-), (2.126) 



l-^lnfA 1 (2.127) 



The unknown coefficients in expression 5*1 is computed by making use of Eq. (2. 122) and 
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solving coupled equations (2.70), therefore one may yield, 

si = T^^^^ fl-i27ri(a;)0aa;)) ^ f [y^, ^j] + V ^-^15^^ 



S\ = ^^Vc~"^^l\ (2-128) 

In the above expression summation over dummy momentum indices is assumed. One can 

find S'y in the same manner by exploiting Eq. (2.72) and using Eq. (2. 128) as an input, 
which leads to 

S'^ = {S'^y + = 1, 4, (2.129) 

with the notations, 

3 



2 



^/4a _ ^ rV[VT,5f-"h. (2.130) 

(a;fe + + + u;^) ^ 

The only divergent contribution up to order g'^ arises from, 

5H^-{Q\[H^,S^IS'M. (2.131) 
After the evaluation of the leading divergent part, we find that 

5// = -#^lnf-'l / d'x^ixl (2.132) 



167r4 \^ 
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which contributes to the two-loop wave-function renormalization Z^. By comparing 
Eqs. (2.127) and (2.132), one may conclude that Z^, = Z^^, as it should be. To finish 
the renormalization up to two-loop order, one should also take into account the contribu- 
tion at order g^. The divergent terms at this level originate from 



[{Va + 1/2Vc + Vb),So],S', |0). 



(2.133) 



After a straightforward but lengthy computation one can obtain the leading divergent 
parts, 

this term should be added to Eq. (2. 124), therefore one can immediately deduce the correct 
total renormalization factor Zg up to two-loop order. 



7 ^ g' 





(2.135) 



One can now immediately obtain the well-known [50] two-loop /3-function and anoma- 
lous dimension by making use of Eqs. (2.127, 2.135). 
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dg 
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IdXogZ^ 
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9 
2^ 
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167r4 



9 



(2.136) 
(2.137) 



It is important to point out that the diagonalization at first order in the coupling constant 
defines a correct low-energy effective Hamiltonian which is valid up to order (/'\ Having 
said that, from Eq. (2.129) one can observe that the non-hermiticity of the S operator 
appears at order g"^ and in a lower order of /i/A. As we have shown, non-hermiticity 
is negligible up to two-loop order (asymmetric terms appear in irrelevant contributions 
(which are non-divergent and vanish as A goes to infinity). We conjecture that, for the 
present model, non-hermitian terms only appear in irrelevant contributions, whatever the 
order of truncation. 



2.10 Example II: Extended Lee Model 

As another illustrative example, we will now apply coupled-cluster RG to determine the 
effective Hamiltonian for an extended Lee model (ELM) up to the one-loop order. 
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We define four kinds of particles, the ^-particle and A/^-particle as two different 
fermions and the 9 and ^ as a scalar boson and anti-boson respectively. Here a{k), a'^{k) 
and b{k), b'^(k) are the annihilation and creation operators which satisfy boson commuta- 
tor rules. The V{p), and N{p), N'^{p) define the fermion sector and obey the usual 
anticommutator rules. The bare ELM Hamiltonian then reads 

H — Ho + Hi, 

Ho = j d'puv{p)V\p)V{p) + j d'pujN{p)N\p)N{p) 
+ j d^k ue{k)a\k)a{k) + J d^k ujg{kp{k)b{k), 

+ X,(2w)-V' j j^^^Nt(j,)V(j,-k)l,(k)+h.c.. (2.138) 



The kinetic energy generically is defined uoik) = ^/k"^ + rriQ where the indice O can be 
either {V,N,9,9). The interaction term in Hi describes the processes; 

V^N + 9, (2.139) 
N^V + 9. (2.140) 

The crossing symmetry become manifest if we take Ai = A2 and equal masses for boson 
and anti-boson. For sake of generality we will ignore crossing symmetry at the moment. 
The Lee model [63] can be recovered if we decouple the anti-boson 9, A2 0. In the Lee 
model the virtual process Eq. (2.140) is not included and thus the A^-particle state remains 
unrenormalized and the model become exactly solvable. 

It is believed that the Lee model is asymptotically free for space-time dimension D 
less than four [64]. With on-shell renormalization one can show that the Lee model for 
D > A (odd D) is ultraviolet stable and not asymptotically free [65]. It is well known that 
such a model in four dimension exhibit a ghost state as the cutoff is removed. The Hamil- 
tonian Eq. (2.138) exhibits two symmetries; it is straightforward to verify that following 
operators commute with H 

B = J d^pV\p)V{p) + J d^pN\p)N{p), 

Q = J d^pN^{p)N{p) + J d^kb^{p)b{p) - J d^ka\p)a{p). (2.141) 
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Clearly i? is a baryon number operator and Q is a charge operator. We assign the charges 
1,0,-1 and 1 to the N, V, 9 and 9, respectively. The sectors of the ELM are labeled by the 
eigenvalue (6, q) of the operators [B. Q). According to our formulation the ultraviolet- 
finite Hamiltonian is obtained by introducing Z-factors, which depend on the UV cutoff A 
and some arbitrary renormalization scale M in such way that effective Hamiltonian does 
not depend on A. The bare Hamiltonian can be rewritten 

H ^ j d'p Z^ZMyU;vip)VHp)Vip) + J d^pZ%ZM^u;MNHp)N{p) 
+ J d^k ZlZMeOJe{k)a\k)a{k) ^ J d^k ZlZM,oJ-e{k)h\k)h{k) 

+ h.c. . (2.142) 
We split the original Hamiltonian in the form of Eq. (2.69); 

Hi = HiT), 
Jo 



2 ^ Hoif ), 

Jim 



'0 Jfi 
f/i /"A 



Va = Vl:, 



'0 J 
r-A /-A 



+ (2(2Xt))V^ ^-^''''"^'^-*^^"''+''"- 



d^pd^k 



58 



CHAPTER 2. HAMRTONIAN RENORMALIZATION GROUPS 



Here p' and k' stand for low momenta (p', k' < jj). If the arguments of an operator are 
all low momenta {p' or k'), this indicates low momentum operators. The arguments in 
H{Jq) and Ho{J^) means that all the momentum integrations involved in Eq. (2.138) 
are running between < < /x for the former and fi < p < A for the latter, respec- 
tively. The configuration space of the high momentum operators are specified by {C/ — > 

yni^n2Qn3^n4^(-t ^ c] ^ ( V^)"! (]Vt)«2 (at)«3 (^t)n4 1 ^^jh m + Us + + n4 = /. 

Aiming at a one-loop expansion the corresponding S and S' operators which preserve the 
symmetry property Eq. (2.141), can be chosen as 

= j d'p'd^kSl^{p')V^{p' -k)h\k) + j d^p'd'kSl,{p')N'^ip' -k)a\k), 

Si = S^N{p') + Sl'V{p'-k')h{k') + SX:V{p' + k')a\k'), 

Si - SlV{p') + Sl^^N{p'~k')a{k')+Sl'N{p' + k')h\k'). (2.144) 

We have ignored the / = 1 configuration, since there are no tadpole type diagrams. ( The 
truncation of 5"/ in configuration space should be consistent with our loop expansion.) 
Here we confine our attention to the elimination of the high-momentum degrees of free- 
dom up to the first order in coupling constant and second order in jj,/ K. The unknown 
coefficients in Eq. (2.144) can be obtained by making use of Eq. (2.73), 
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{2{2T:fuje{k)YI\ujN{p' - k) + uje{k)) ' 
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{2{27ryujs{k)y/^{ujv{p' -k)+ ug{k)) ' 
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Xiuvip') 
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A1A2 


2{2T:Y{uj-e{k')uJe{k)Y'\ujN{p' - A;) + a;,(A;))2 
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2{2'nY{^e{k')ue{kY^/\ujN{p' - k) + uJe{kY^ 
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(2.145) 
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It is easy to observe that Eq. (2.74) will be satisfied if we require S'^ = Sl^, since, up 
to the first in A the similarity transformation introduced in Eq. (2.67) remains unitary. 
Equally one could use Eq. (2.75) to obtain S', it is obtained that 5*0 = Sq and S[ = 0, we 
will show that the renormalization feature of our model up to this order will remain un- 
changed, however the effective low-energy Hamiltonian will be different. As was already 
pointed out, this is because Eq. (2.75) requires a different truncation scheme. The effec- 
tive Hamiltonian is now produced by plugging the S and S' defined in Eq. (2.144) into 
Eq. (2.67). With naive power-counting one can identify the potentially divergent terms. 
At the lower order of expansion, the divergent term is (0|[l/4, -S'oJlO), the divergence in 
this term arises from a double contraction of high-energy fields. 

At this step the contributions of the terms [Vb, So{Sq)] and [Hi, So{S'q)] are zero, 
after projection on to the high-frequency vacuum. There is one other divergent term, 
{0\[Vc, S'o\\0), but this is harmless and will be cancelled out by {0\[[H2, So\,S'q] |0). We 
thus obtain 

From this expression one can immediately deduce the renormalization factors Zj^v and 
ZruN' we take (vq — for ^ mo, therefore 

Zmv — 

There is no mass renormalization for 9 and 9 and accordingly there are no vacuum polar- 
ization type diagrams. Thus 9 and 9 remain unrenormalized, Z^ — Z^ng — Zg — Z^^ — 1. 
The other contribution of at one-loop which are not zero after projecting on to vac- 
uum come from 

m{\) = -(0| [[Hi, Si], S'^] + (0| [[Hi, Si],S'i] |0). (2.148) 

The divergent contribution emerges from the first terms, the leading divergence of this 
expression is logarithmic which means that we can neglect the difference between p and 
k — p' (for the divergent contribution only). After evaluating a momentum integral we 



1 + 4(A-A^). 



87r2 

A 

87r2 



(A-/.). 



(2.147) 
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finally get, 

SH = 



A? 
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Af + Aj 
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2 In 
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(2.149) 



From this expression we deduce the renormalization factor Zx^ , Zx^ , Zy and Zat: 



r.1 A2 A 

Z?. = 1 + n-, 
^ 167r2 /i ' 



Zx^ ^Zx,^l + 



Xl + Xh A 
— ^In— . 



(2.150) 



It is obvious from the equations above that one can define renormalized coupling constants 
in terms of the bare couplings and wave function renormalization 

Xi^X^jZyZj,, i^l,2. (2.151) 

This definition corresponds for Ai = A2 with the renormalization introduced to compute 
the T-matrix for the N — 9 interaction [63, 66]. The one-loop /3-function and anomalous 
dimension 7 are 
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(A? + A^), 



1,2. 



(2.152) 



Since the fixed points of the theory are the zero solutions of the /3-function, one imme- 
diately identifies the trivial solution Ai = A2 = (we ignore the nonphysical imaginary 
solution). It is now obvious from Eq. (2.152) that — In — ^ only at trivial so- 
lutions of the /3-functions. This result is in correspondence with property that for real 
field theories the 7-function is not zero unless at trivial fixed point of the theory [67]. 
Now to investigate the behaviour of the theory at high momentum, we must compute the 
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momentum-dependent effective coupling constant Xi{k) and A2(A;) by 

= Ai(^)|,=i = Af 



(2.153) 



where Af and Af are dimensionless physical renormalized coupling constants defined at 
the renormalization scale k = 1. The coupled equations (2.153) can be solved by going 
to polar coordinates r'^{k) = Xl{k) + Xl{k) and e{k) = tan"^ 

V - 
XAk) — — = sin 9, 

yi- (167r2)-if2inA; 

X2(k) = ' cos 6, (2.154) 

^ ^ yi-(167r2)-if2 1nA; 

where f and ^ denote the value at the renormalization scale. The behaviour of the ELM 
in the deep-Euclidean region is obtained by allowing k oo. From Eq. (2.154) one 
observes that Ai(/c) and X2{k) in this region are imaginary. This means that the effective 
Hamiltonian is non-hermitian and the theory generates ghost states when the cut-off is 
removed. The ghost state appears as a pole in V and A^-propagators. Since a theory is 
said to exhibit asymptotic freedom if (i) ^|a(oo) < (ultraviolet stability at the fixed 
point A(oo) and (ii) A(oo) = hmfc^oo X{k) — 0, Eqs. (2.154) indicate that the ELM can 
not exhibit asymptotic freedom at D — A. 



2.11 Conclusion 

In this chapter we have reviewed the merits and shortcomings of several approaches for 
the construction of the effective field theories in the Hamiltonian framework. 

We have outlined a new strategy to derive effective renormalized operators. The for- 
mulation is not restricted to any quantization scheme (e.g., equal time or light cone). The 
effective low-frequency operator is obtained by the condition that it should exhibit de- 
coupling between the low- and high-frequency degrees of freedom. All other irrelevant 
degrees of freedom like many-body states can be systematically eliminated in the same 
way. We have shown that the similarity transformation approach to renormalization can 
be systematically classified. The non-hermitian formulation gives a very simple descrip- 
tion of decoupling, leading to a partial diagonalization of the high-energy part. 

The techniques proposed are known from the coupled cluster many-body theory. We 
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fully utilized Wilsonian Exact renonnalization group within the CCM formalism. Our 
approach invoke neither perturbation nor unitarity transformation. It can be conceived as 
a topological expansion in number of correlated excitation of the high-energy modes. We 
showed that our formalism can be solved perturbatively. In this way, it was revealed that 
diagonalization at first order in coupling constant defines a correct low-energy effective 
Hamiltonian which is valid up to the order A^. 

We showed that non-unitarity representation inherent in our formulation is in favour 
of economic computation and does not produce any non-hermiticity in the relevant terms. 
One can show that the non-hermiticity of the effective Hamiltonian is controllable and 
might appear in higher order which is beyond our approximation or in irrelevant terms 
which can be ignored in renonnalization group sense. We argued that our formulation 
is free of any small-energy denominator plaguing old-fashion perturbation theory. We 
showed that the non-hermiticity of the coupled-cluster parametrization leads to the com- 
patibility of the formulation with the Hellmann-Feynman theorem and it also induces a 
symplectic structure. More importantly, it provides a simple framework for the renonnal- 
ization of an arbitrary operators. Notice that all these features are connected with each 
other and one can not give up any of them without spoiling the others. One may conceive 
that the non-hermiticity adds a auxiliary (non-physical) sector to the physical phase space, 
thereby, it makes the whole phase space geometrically meaningful and moreover it gives 
enough room to keep the formulation to be conformed with the cluster decomposition 
property and Poincare invariance regardless of a regularization method. There is a long 
tradition behind such approaches, of course with different motivation, e.g., in the BRST 
formulation, the phase space is enlarged by anti-commuting canonical coordinates, an- 
other example is the bosonization of spin algebraic or fermionic system, where one maps 
the original Hilbert space of the system into a boson Hilbert space which turns out 
to be larger than original Hilbert space, in the sense that physically realizable states in 
the original space map into a subspace of H^. Interestingly, in this approach as well, the 
boson Hamiltonian can be either Hermitian or non-Hermitian [34]. 

Notice that our RG method is non-perturbative although we have already shown that 
perturbation expansion in coupling can be easily implemented [32, 33]. We successfully 
applied our RG formalism to compute perturbatively the effective Hamiltonian for cf)'^ and 
extended Lee theory up to two- and one-loop order, respectively. We have employed a 
sharp cutoff, however this idealization should be removed since generally it may lead 
to pathologies in renonnalization, since it induces non-locality and moreover potentially 
violates the gauge symmetry. 

One of the key features which has not yet been exploited is the non-perturbative aspect 
of the method; it may well be able to obtain effective degrees of freedom that are very dif- 
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ferent from the ones occurs at the high-energy scale. This is a promising avenue for future 
work. Another interesting question is that the connection between our non-perturbative 
truncation scheme SU B (J\f, Ai) and other non-perturbative scheme e. g., the large Nc 
truncation. A systematic scheme which relates the large Nc limit with an approximate RG 
equation remains yet to be discovered. 
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Chapter 3 

Renormalization problem in many-body 
system 

3.1 Introduction 

In this chapter we shall concentrate on the renormalization problem in many-body theory 
of nuclear matter. As we already pointed out, renormalization of many-body system in 
a truncated (in number of particles) space is problematic, the so-called Tamm-Dancoff 
problem. The question is how can one renormalize the many-body problem equations 
obtained by non-perturbative approaches (such as the CCM, Brueckner's reaction matrix 
(or G-matrix) theory [68] and etc.,) in a truncated Fock space. 

In the last chapter we showed that a renormalized effective interaction in small num- 
ber of particles can be obtained by imposing certain decoupling conditions between the 
model- and excluded-spaces. In this sense, Feshbach formalism is in contrast with the 
effective interaction theory since it is not derivable from such decoupling conditions. 
Notice that the energy-dependence of the Feshbach formalism (and any Green func- 
tion type formulations, e.g., Schwinger-Dyson resummation, Faddeev approach) emerges 
from the fact that the effects of the excluded Hilbert space is taken into account by a 
"quasi-potential", while in the effective interaction approach the latter is taken into ac- 
count by imposing a certain decoupling conditions. Therefore a given truncated Hilbert 
space becomes independent of the remaining sectors and accordingly it can in principle 
be described by an energy-independent prescription. 

In order to clarify the differences between an energy-dependent and an energy-independent 
formulation, here we investigate how can one resolve the renormalization problem by 
fully utilizing the Feshbach projection operator technique [56] in the framework of the 
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CCM. Therefore, we pursue an inverse of the EI approach. With a field theoretical con- 
sideration, we show that the coupled-cluster formalism by means of Feshbach projection 
technique leads to a renormalized generalized Brueckner (^'-dependent) theory. 

The Feshbach projection technique was introduced to treat nuclear reactions with 
many channels present. It was originally formulated under the assumption that the num- 
ber of elementary particles involved is conserved, however clearly this is not that case for 
field theory. An extension of the Feshbach formalism has been developed for the pion- 
deuteron system [69] and general pion-nucleus reactions [70]. This technique bears some 
resemblance to Okubo's methods [71], which is consistent with meson field theory, but is 
developed in terms of an effective Schrodinger equation so as to remain in close contact 
with conventional nuclear physics. This approach was already pursued by Schiitte and 
Providencia [72] in the framework of the CCM for the Lee model. However, in the Lee 
model, because of an inherent Tamm-Dancoff approximation (which limits the number of 
mesons present at any instant), and exact solvability of the model, the issue of renormal- 
izibility of the nuclear matter properties is unclear. Here, we follow their approach in an 
extended version of Lee model which is not exactly solvable and does not display these 
trivialities. The renormalization of the extended Lee model in the few-body sector was 
already investigated in section 2.10. 

Notice that the CCM as introduced by Schiitte and Providencia based on Rayleigh- 
Ritz- Variational principle (e.g., see [72]) is different from the CCM introduced by Arpo- 
nen and Bishop from the standpoint of variational principle and the Hellman-Feynman 
theorem [3, 4]. As we have already illustrated in the first two chapters, we believe that the 
latter has more advantages and is more suitable for a field theoretical application since it is 
naturally embedded in the modem effective field theoretical framework. Having said that 
our main goal in this chapter is to introduce other quantum many-body theory techniques 
and challenge if they are adoptable for a field theoretical application. 



3.2 The extended Lee Model 

The extended Lee model (ELM) is a simple model connecting elementary and compos- 
ite particles [63]. Although this model is not a chiral model but it exhibits many field- 
theoretical features. We define four kinds of particles. The V- and N- particles, two 
different types of fermions, and the 9 scalar boson and the 9 anti-boson. We take for the 
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Hamiltonian the expression: 



a p k 

H'l = T.^U^lN^^^k + Y.^U'^l^-^k' + h.c. . (3.1) 

a/3k aPk' 

Here afc,a|, and hk,h\ are the annihilation and creation operators which satisfy boson 
commutation rules. The V^, Va and A^^, Nj^ define the fermion sector and obey the usual 
anticommutator rules. The a, k and k' are abbreviations for all quantum numbers (such 
as momentum, spin and isospin, etc). Within our formal investigation, we leave open the 
specification of E^^, E'^ and ujk- It can be taken as either a relativistic or non-relativistic 
expression. The bare kinetic energies E^ and E^ are renormalized to E^ and Ep by the 
interaction. The matrices W^^^ and are the bare interaction strength renormaliz- 

ing to Wapk and Wapk', respectively. The interaction strength is defined by the kind of 
bosons exchanged (the scalar, pseudoscalar or vector bosons). The exchange of higher 
spin (J ^ 2) bosons, such as the f(1260), A2(1310), f (1514) and g(1680) seems to have 
little influence on the low-energy NN data. The reason is that their poles are located 
far away from the physical region. In other words, they give rise to the contributions of 
the very short range which are masked by the very strong repulsion coming from vector 
meson-exchange. These contributions, however, are essentially masked or parametrized 
by form factors necessary to regularize the one-boson-exchange diagrams. The interac- 
tion term in H describes the process 

V^N + e, (3.2) 
N^V + e. (3.3) 

The Lee model can be recovered if we decouple the anti-boson 9, W^^^,, — > (i. e. to 
remove crossing symmetry) [73, 64, 65]. In the Lee model the virtual process Eq. (3.3) 
is not included and thus the A/^-particle state remains unrenormalized and the model be- 
comes exactly solvable [73]. We have already shown that the ELM can not exhibit asymp- 
totic freedom at space-time dimension D=A. It is well known that such a model in four 
dimension exhibit a ghost state as the cutoff is removed. It has been shown that the 
composite-particle theory, the meson pair theory of Wentzel [74], can be obtained as 
a strong-coupling limit of the ELM, in which limit the wave function renormalization 
constant of the F-particle vanish [63]. One can find two operators commuting with the 
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Hamiltonian, 

B = 5] V;Va + iV^iV^, (3.4) 



a j3 

P k' k 



a\ak. 



Clearly i? is a baryon number operator and Q is a charge operator. In setting up the 
charge operator, we have assigned the charges 1,0,-1 and 1 to the particles A^, V, 9 and 
9, respectively. The sectors of the model are labeled by the eigenvalues (6, q) of the 
operators {B,Q). The most trivial sectors are: (0,0) = \o), the physical vacuum is 
thus the same as the free particle vacuum for this model, (0, 1) = {6^,10} = \k')} and 
(0, — 1) = {al|0) = \k)} which define respectively anti-meson and meson states, these 
states stay unrenormalized. The sectors (1,0) and (1,1) contain: 

(1,0) = {\/t|0) = |«),iVX|0) = \(3,k),\akk'),...}, 

{1,1) = {\l3),\ak')M3kk'),...}. (3.5) 

In the Lee model stays unrenormalized. This makes the model exactly solvable, but 
this is not the case in ELM. The renormalization of such a model was already studied in 
section 2.5. The one-loop renormalization of the interaction strength is defined as 



In one-loop order, the physical neutron state is a bound state in the (1,0) sector 
with H\il>a) = Ea\i'a), where Ea is renormalized ^-particle energy. We note that in 
the limit that the coupling constant vanishes, the ^-particle state \'4>a) goes over into jo;). 
Therefore, the simplest form of a bound state in the (1,0) sector at one-loop order can be 
written as 

\i^.)^\a) + Y.^{P,k)m, (3.7) 

I3,k 

where the unknown coefficient 0(/5, k) is determined by requiring that \ipa) be an eigen- 
function of Hamiltonian with eigenvalue E^. After straightforward calculations, one ob- 
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Figure 3.1: Diagrams in effective mass operators, h{z). Solid lines and dashed lines 
denote V- and A^-particles, respectively. Dotted and dash-doted lines denote 9- and 9- 
particles. 



tains 

M^) = E , ^T', - (3.8) 

The same argument can be applied for the sector (1,1) and one can obtain in lowest order, 
the physical proton state \il}fi) as a bound state in this sector which can be shown to obey 

1^/3) = 1/3) + E p ;r ^ 



Ef3^ Ep + hi3{Ep), 



haiz) 



The and hi3{z) are the mass operators, and show the off-shell contribution to the 

self-energy, see Fig. 3.1. The mass renormalization in the model is now performed by 
adding corresponding terms as counter terms to the Hamiltonian (these terms will not 
change the conservation of B and Q operators). We assume that the parameters of the 
model are chosen in such a way that there exists one bound state for each a and (3 . 
The form factors contained in W^pk^Wpak') are assumed to make Zq, > 0, > and 
consequently the coupling constant renormalization finite. 



3.3 Nuclear matter equations 

We now wish to consider the binding energy problem of H for the sector (2n, n) with 
n — > oo so-called N — V matter. We assume that the non-interacting ground-state wave 
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function of this system be the Slater determinant built up by an equal number of V- and 
N- particle up to a Fermi momentum Pp, 

10)= n ^aNl\0). (3.10) 

a,l3<PF 

We denote by a{b) occupied ^-particle (A^-particle) states, by A{B) the unoccupied states 
and by a{f3) either states. We write the exact correlated ket ground state, in a coupled 
cluster formulation, as 

|^) = e^+^|0). (3.11) 

Here the cluster operator is separated into two parts, one, S, without mesons which is 
obviously the same as in the CCM with a phenomenological potential and R with mesons 
contributions which contains the additional parts originating from quantum field theory. 
Here we define bra ground state as a hermitian conjugate of ket ground state in every level 
of truncation. Using the symmetry properties Eq. (3.4) and momentum conservation the 
general form of R and S are 

n/2 n/2 
i=l j=0 



Sij = Tj^-j^ {Ai..AiBi..Bj\Sij\ai..aibi..bj)A 



ABah 



X Vl..VlNi^..Ni^N,^..N,,Va,..Va„ 

n/2 n/2 k k' 



k=l 1=0 m=l n=l 



Rklmn — , |2J|2 {^l--^k-m+nBl--Bi_^^_nKi..KjnK'^..K[\Rklmn 



k\2l\ 

ABabKK' 

X Va,..Va,al..alA'rbl'„- (3-12) 

The subscript A of the R- and yS-amplitudes stands for antisymmetrization of the states. 
For obtaining the binding energy of N-V matter one should solve the Schrodinger equa- 
tion by making use of the definitions in Eqs. (3.10)-(3.12), therefore 

e-^e^He^e^\(i)) = E\4)). (3.13) 
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Projecting this equation onto the complete orthonormal set of states 



\<P),Nj,Vaal\<P),VlN,bl,\<P),.. , 



(3.14) 



one obtains a system of coupled integral equations of which the first determines the bind- 
ing energy and others fix the corresponding amplitudes (R) and (S). These coupled 
integral equations from effective interaction view point are the decoupling conditions and 
leads to energy-independent prescription. In order to derive the standard Brueckner the- 
ory we apply the Rayleigh-Ritz- Variational principle. This type of the CCM formulation 
was originally proposed by Providencia and Shakin [75]. This is of course different from 
coupled-cluster theory introduced by Arponen and Bishop from the standpoint of varia- 
tional principle and the Hellman-Feynman theorem [3, 4] (the so-called NCCM scheme, 
see chapter 1.2). 

One of our goals here is to show that in the consistent truncation scheme the many- 
body problem does not require further renormalization and all the bare parameters intro- 
duced in definition of the wave function can be fixed in the nucleon-nucleon scattering. It 
is clear that Eq. (3.12) needs to be truncated. Aiming at a two-hole line and one-meson- 
exchange expansion of the ground-state energy, we choose the following terms for S and 
i?fromEq. (3.12). 



R ^ J2 RaBkNlVaal + J2 RAbk'VlNbbl + I /^RsB'abkNlNl NbVaal 



All other term vanish up to this order because of the symmetries of the ELM Hamiltonian. 
Now we define cluster expansions of expectation value of the Hamiltonian with respect 
to lip), which is a straightforward generalization of the standard expansion [75, 76]. We 
neglect three-body and higher-order cluster. The one-body, two-body, etc., correlated 
wavefunction are defined as 



+ J2 ^/^RAA'abk'VlvlVaNbbl = Ri + R2 + R3 + Ra- 



(3.15) 



|V'a) = e^-^Xt|o) = (i + i?^)V;t|o), 
|V^,) = e^+M|0) = (l + i?2)iV,^|0), 
|« = e^+XVi|0) = (l + 5 + i?KVi|0). 



(3.16) 



The expectation value of the Hamiltonian is then obtained by and shown in Fig. 3.2 
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9 = u = 



+ 




+ 




+ 




R = ■ = 



+ 



+ 




+ 



+ 



= N-particle 



= V-particle 



= 6 -particle 



Figure 3.2: The symbols x and o stand for kinetic energy and different parameters of wave 
function(C, C", D, F, F'), respectively. The symbol • denotes the interaction coupling W. 
The out-going arrows form o denote V, N, ^-particle and in-coming arrows denote V, A^- 
hole and ^-particle. The abbreviation C.C impUes the complex conjugate terms. 
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E ^ ^ti^lf'^ = ^^P^ + ^^P^ + habPaPb, (3.17) 

a b ab 

where the Hamiltonian cluster integrals have been introduced as 

K = ('0a|-H'|'0a), 
hb = (ipblHll/jb): 

hab = {lpab\H\lpab) " ^a^^ft " hUa, (3.18) 

here Ua and rib are the norm of [ijja) and IV^^) respectively. We disregard the other cluster 
integrals like haa' because they would only get contributions from two-mesons intermedi- 
ate states. This makes sense in consistent one-meson exchange approximation. According 
to the selective summation made in the expression Eq. (3.17), the occupation numbers pa 
and pb should satisfy the following algebraic equations: 

riaPa + Y '^abPaPb = 1, 
b 

ribpb + Y ^<^bPbPa = 1, (3.19) 

a 

with notation, 

riab = ii^abli^ab) " ^aHfe. (3.20) 

The basic idea for the treatment of {ilj\H\^) / (^|^) is to expand it term by term using the 
Wick-rule to keep track of all possible contribution. The expansion of {^\H\ip) and (^|^) 
can therefore be characterized by a set of suitable diagrams and one can check that the 
division by (^|^) cancels out all contribution from disconnected diagrams (the linked- 
cluster theorem). Now we minimize the quantity E with respect to the constrains given 
by Eq. (3.19) and denote it E, thus we multiply Eq. (3.19) by the Lagrange multipliers 
and €b and subtract them from Eq. (3.17), then by using the definition Eq. (3.18) we find 

E = Y,PaiE'a-ea) + Ypb(^b-^b) + Yp-(^\^l(^0-^a)Rl 

aba 

+R\h'I + H'lR^\a) + Y Pb{b\Rl{K - ^b)R2 + RIh^i + H^iR2\h) 

b 

+ Y PbPa{ba\S\H^ - e„ - eb)S + i?S(//° - e„ - eb)R3 

ab 

+rI{H^ - e„ - eb)R4 + S^H'^R + R^H'^S\ab). (3.21) 
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It is convenient to treat the pa and pb as independent variables, therefore we minimize E 
respect to i?*^^, R*^^^^,, S\b^^, R*BB'abk^ R*AA'abk^ Pa,Pb which respectively yields: 



RaBkiE^I, +CUk- e„) + W^Bk + Yl PbSABabW'/,, = 0, (3.22) 

Ab 

RAbk' (E'a + OOk' - eb) + W%,, + J2 PaSABabW^*^,, = 0, (3.23) 

Ba 

{AB\{H^ - e„ - et)S + H^R\ab) = 0, (3.24) 

{BB'kliH^ -Ea- eb)Rs + H'^Slab) = 0, (3.25) 

{AA'k'\{H^ - e„ - eb)i?4 + H^S\ab) = 0, (3.26) 



e„ = + {a\Rl{H^ - e„)i?i + rIh'} + H''jR^\a) + ^p6(a6|(i?l + Rl)H^S\ab), 

b 

eb = + {b\Ri{H^ - eb)R2 + RIh'^ + H'^R2\b) + J2Pa{ab\{Rl + Rl)H'}S\ab). 

a 

(3.27) 

Equations (3.22)-(3.27) represent a set of unrenormalized Variational coupled integral 
equations for N — V matter, consistent with two-hole-line truncation. 



3.4 The renormalized nuclear matter equations 

The method that we will use to reformulate Eqs. (3.22)-(3.27) in order to obtain renormal- 
ized equations for iV — \^ matter is based on the Feshbach projection operator formalism 
[56]. Projection operator techniques have also been used to analyze several sectors of the 
Lee model [77]. 

The total energy of N — V matter can be obtained by making use of Eq. (3.21) and 
Eqs. (3.24)-(3.27) and exploiting the constrains Eq. (3.19); 

^ = + S]^" - YPaPb{ab\{Rl + Rl)H'}S\ab). (3.28) 

a b ab 

We can now directly add mass counter terms to the bare Hamiltonian in Eqs. (3.22) and 
(3.23) to dress masses. We expand the last term in the binding energy Eq. (3.28) in terms 
of the coupling constants and S. This form will be used later to achieve full renormaliza- 
tion for the binding energy since there is a subtle relation between S and the renormalized 
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coupling constant, 

N f I. — > n.rti, 



E - E^« + E^^- E PaP^^%^^SABa, (3.29) 



a b abk'AB 

PaPh ^ '^ABab— 7^ PaPbPa' ^ >^ AB' a'b>^ ABab 

abkAB ^a-EB-UJk aa'bt^BB' ^b ' Ea - i^k' 

PaPbPb' ^ JA'BaV'^ABab- 

abVkAA'B ea-Es-UJk 

The last two terms in Eq. (3.29) are a three-hole-line contribution which we ignore, con- 
sistent with our aim for a two-hole-line expansion. We introduce the projection operators 

Q = ^\AB){AB\, Q' = l/2^\BB'k){BB'k\, (3.30) 

Q" ^l/2j2\AA'k'){AA'k'\. 

Inspired by the projection technique of feshbach, one can use the above definition to 
combine Eqs. (3.24)-(3.26), 

H^-z + QH^ (q'^—^Q' + Q»^—^Q")h'}q] S\ab) = -QH'i{R, + R2)\ab), 

(3.31) 

where z — €a + We can show that the right hand side of Eq. (3.31) can be reduced by 
using Eqs. (3.22) and (3.23) 

{AB\H^R,\ab) = J2KbkR aBk 

k 

_y W^AbkKk y SA'BabW'^,,W%, ^^^^^ 

^ea-EB-UJk ea-EB-UJk ' 

{AB\H^R,\ab) = ^W',^,,RAbk' 

k' 

^ ^aBk'^Abk' y p SAB'abW^Bk'^a^B'k' q 33^ 

The second term in Eqs. (3.32) and (3.33) are a three-hole-line contribution which are 
again ignoreable in our approximation. Let us introduce the unrenormalized effective 
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Figure 3.3: The effective two-body interaction, the lines are defined in Fig. 3.1. 



interaction U{z) 

{a(3\U{z)\a[3') = («/5|i/°^=^i/°|«'/?')iinked 

z- Ho 

a/5 k a pk afi k' a pk' 

^ z- ep- €0 - ujk ^ z-e^-e^' - Uk' ' 
Hg = Y. eaVlV^ + 5^ epVlVp + ^ u:^{a\ak + h%). (3.34) 

a p k 

The mass renormalization terms has been taken into account by using renormalized masses 
in the propagators. (Notice that = and = Eb)- This definition correspondences 
to the quasi-potential in a Lippmann-Schwinger type equation if one replaces H^ Hq. 
Having made use of defined effective interaction U (z) one can show 

{AB\H%Ri + R2)\ab) ~ {AB\U {z)\ab) . (3.35) 

This expression is consistent with the two-hole-line approximation. Now we decompose 
the "p-p" and "n-n" interaction involved in Eq. (3.31) after adding mass counter terms in 
bare Hamiltonian, 

Qi/?(Q'^— Q' + Q"—^Q")H'lQ = QU{z)Q + h{z) +q{z), (3.36) 

z — Hq z — Hq 

where the operator h(z) and q{z) are defined via Eqs. (3.8,3-9) as 

h{z)\AB) = [hAiz - Eb) + hB{z - Ea)]\AB), (3.37) 

q{z)\AB) = -{y + y ) \AB). 

^ ' \^ Z-E,-EB-COk ^ z-Ea-EA-UkJ^ ' 

In diagrammatic language Eq. (3.36) contains self-energy diagrams which contain un- 
occupied intermediate states because of the operators Q' , Q" , This is achieved by adding 
q{z) to h{z). We now introduce the operator B{z) to accomplish coupling constant renor- 
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malization. Its structure emerges from renormalization of finite sectors of ELM, 

B{z)\al3) = Zl{z - E(,)Zl{z - (3.38) 



with notations. 



ITI/O 12 

here z = + Ep. The benefit of B[z) can be manifested by following factorization 
property 

Q{z - - h{z))Q = Q{z - Ho)B{z)Q. (3.40) 
One may use B{z) to show that U{z) is related to renormalized effective N — V potential 

U{z) = B{zf'''U{z)B{zf/''. (3.41) 

Therefore, the renormalized effective two-body N — V potential can be readily found ( 
see Fig. 3.2) 

{a^\U[z)\a /?) = -> ^ 

V z-e,3-ef^ -ujk 

where we have introduced twice-subtracted "dressing factor" defining as 

Dressing factor appears independent of truncation made and describes off shell correction 
of a self-energy of V and N particle. In Lee model Vfjiz) = 1, since iV-particle stay un- 
renormalized. Twice-subtracted dressing factors remain finite even without form factor, 
due to the occurrence of a cubic energy denominator, this can be verified by expanding 
Eq. (3.43) using relativistic expression for Ea,Ej3,ujk, ujk' and interaction strength. Now 
one can rewrite Eq. (3.31), by making use of Eqs. (3.35, 3.36) and exploiting the factor- 
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ization property Eq. (3.40) 

{Ho -z + q{z) + U{z))B{zy^^S\ab) = -QU{z)B{zf'^\ah) , (3.44) 
where q{z) is renormalized two-body operators obtained from Eq. (3.37) 

\WAbk?r\{z - EbYUz - Ea) 



q{z)\AB) --[Y. 



z — Eh — Eb — UJk 
y\W.B.K{z-EByBiz-EAh 

^ Z-Ea-EA-UJk' /' ' 

ak' 

It is observed the operator B{z), transfers unrenormaUzed correlation function to renor- 
maUzed one. This is due to the simple structure of vertices renormalization in this model. 
We now turn to Eq. (3.27) for and e;,, one may use Eqs. (3.22,3-23) to obtain 

= + RgBk^aBk + PbRAbk'^aBk'^ABab 

Bk bk'AB 

Bk ^ bkAB a B K 

I ^ ^^aBk'^^Abk' c I \^ ^ ^ ^^aBk'^^ a'B'k' c* q 

2-^ PK PO —'='ABab^ PbPa'- —'^AB'a'b'^ABab, 

bk'AB OAK a'bk'ABB' " A « 

^b = E^ + ^RAbk'^Abk' + XI Pa^*aBk^bAk^ABab 
Ak' akAB 

|TT-n 12 TI-O TT-0* 

— ??0 I \^ I ' ' Abk' I I \^ ^ Abk' ' ' aBk' Q 

- ^b+ 2^ p — + Pc^- P —^ABab 

^e^-E^-Uk' e^-E^-Uk' 

, ^ ^blk^aBk c I \^ ^ ^ ^Mk^A'b'k Q* Q 

^ 1^ P'^Z PO —^ABab+ 2^ PaPb'- —^A'Bab'^ABab- 

bT^B ^a-E^-Uk ea-Ej,-Uk 

(3.46) 

The last term in ea{eb) in Eq. (3.46) are a three-hole-line contribution. From Eqs. (3.8) 
and (3.9) it is obvious that the first two terms in ea{eb) in Eq. (3.46) are the renormalized 
energy of iV(y) particles if we add the contribution of intermediate occupied a (6) states 
in self-energy. The corresponded terms should be subtracted which contribute to the 
definition of €a{€b). Having introduced mass counter terms, one can eliminate 5" by using 
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Eq. (3.44) and obtain the renormalized equation for E, and ea, 

E = ^€a + ^€b -^pbPa{ab\G{€a + €b)\ab) + PaPb{ab\U + eb)\ab), 



ab ab 



(3.47) 

ea^Ea + e- + Pb{ab\G{ea + eb)\ab), (3.48) 

b 

eb^Eb + Sb + Pa{ab\G{ea + eb)\ab), (3.49) 

a 

G{z) = U{z) + U{z) ^ G{z), (3.50) 

z - Hq- q[z) 

}^pb{ab\U{ea + eb)\ab), (3.51) 

bk ^a-Eb-Uk Y 

> Pa(aoF(ea + e6)|a6), (3.52) 

66 — £/„ — ^fc' 

afc' a 

where and e?, are given by self-consistent equation of the Brueckner type and G{z) is 
the solution of a Brueckner-Bethe-Goldstone equation [68, 78]. 



3.5 Conclusion 

The correction to the standard many -body scheme emerges in this model in different steps 
which are connected with the following physical effects. We have to replace the V'^^^ 
by U{z), which accounts for the renormalization effects of the boson propagating in the 
N — V matter system. This can be understood since fermions feel the average potential 
given by ea and during the time when a boson is exchanged vaN — V matter [78]. In 
scattering formalism (T-matrix) we have z — Ea + Ep, the total energy of N — V system, 
whereas there is a medium effect in N — V matter bringing the 2;- value to 2; = €a + €b. 
The energy denominator of the Brueckner Eq. (3.50) contains the correction q{z) which 
is not present for standard N — V scattering. This correction is due to the effect of the 
Pauli principle on the self-energy of the fermions. The comparison between many-body 
solution in the Lee model [72] and the ELM reveals following distinctions: 
The simplest correction in the ELM is due to the self-energy diagrams of iV-particle 
and Pauli principle on this diagrams which are manifested through the occurrence of the 
corresponding dressing factor r^(2;) and q{z) respectively, whereas these corrections in 
the Lee model do not appear because A^-particle remains unrenormalized . In the Lee 
model the lowest order of would appear in a three-hole-line expansion while in the ELM 
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Figure 3.4: The symbol x stands for the renormaUzation counter term which is shown on 
the left land side. 

this term arises in two-hole-line contribution on an equal footing with lowest order of ea- 
The pb renormalization correction of the Lee model (see Fig. 3.4) goes to contribution of 
mass renormalization of A^-particle in the ELM. Therefore the renormalization correction 
of occupation numbers gets contribution from higher order which are not generated in 
our approximation. It is notable that the ELM due to possessing the crossing symmetry 
preserves the symmetry between V and N part of the renormalization correction, whereas 
this is not the case for many-body problem within Lee model. 

This presentation is conclusive that coupled-cluster theory without the decoupling 
property in systematic truncation scheme leads to derivation of generalized Brueckner (in- 
dependent) theory which includes renormalization correction originating from medium 
effect. We showed as well that a combination of the coupled-cluster theory and Feshbach 
projection technique provides a powerful method to renormalize quantum many-body 
problem in a truncated Fock space. 

3.6 Some final remarks 

We have employed the coupled-cluster theory in various versions and investigated the 
merits and short-comings of such techniques for field theoretical applications. We showed 
that the CCM version introduced by Arponen and Bishop can be easily adopted with 
Wilsonian renormalization group and provides very strong tools for describing (non)- 
perturbative phenomena (please see the conclusion of the last chapter 2. 1 1). On the other 
hand, the CCM version introduced by Schiitte and Providencia can be equipped with 
the Feshbach projection formalism and produces a renormalized generalized Brueckner 
theory. One should bear in mind that in the latter, the famous "small-energy denominator" 
problem is not avoidable and is not at all clear how systematically high-energy modes can 
be integrated out for a very complicated system. 

In the rest of this work, we employ an effective QCD model as can be constructed 
from the techniques of the first part of this thesis and address with greater details, non- 
perturbative phenomenological phenomena such as nucleon and diquark solutions, con- 
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Chapter 4 



QCD properties and effective quark 
chiral models 

It is believed that the strong interactions are described by a quantum field theory known 
as quantum chromodynamics (QCD) [50, 79]. QCD from many aspects is an unique 
theory. Quantum electrodynamics (QED), and its expansion to the electroweak standard 
model of particle physics, is also a quantum field theory. However, QED breaks down 
at short distances and is not well-defined. QED is renormalizable theory but it loses all 
his credibility as we approach to Landau pole. On the other hand if the cutoff goes to 
infinity, QED becomes trivial. QED is not the only theory with a Landau pole problem, 
every theory which is not asymptotically free suffers from this problem. The quantum 
field theory of gravity obtained from general relativity suffers from nonrenormalizability. 
QCD is the only known theory which is free from such problems. QCD only needs few 
parameters to be defined completely, one universal coupling strength and one mass for 
each kind of quark. 

In following chapter we introduce part of QCD which is relevant to nuclear physics. 
We refrain from discussing all details since it can be found in many quantum-field theory 
textbooks, see e.g., Refs. [50, 79]. 

4.1 QCD Symmetries 

In this section we will introduce the underlying symmetries of QCD. The Lagrangian of 
QCD is given by [50, 79], 

L = q{i^>^d^ - im')q - -^{F^^f + gq^A^q, (4.1) 

83 
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where q is the quark field which is defined in the fundamental representation of the color 
and flavor group, and the conjugate Dirac field is defined as g = ^^7"^. The gluon field 
matrix = A^A"/2 is defined in the fundamental SU {Nc = 3) representation, A" being 
the generators of the gauge group which satisfies [A"/2, A''/2] = if'^^'^X'/2. We define g 
as strong coupling constant. The field strength is given by 

F;, = d^Al - d^Al + gr^^A^Al. (4.2) 

The non-Abelian nature of QCD is manifested by the quadratic term in the gauge field 
strength. The color and flavour indices of the quark field are suppressed. m° is the current 
quark mass which is not directly observable if QCD confines quarks. The current quark 
mass is color-independent and can be brought diagonal in flavour space. There are six 
flavours of quarks, each of which has a different mass. The three light quarks are called 
up (u), down (d) and strange (s), while the three heavy quarks are called charm (c), bottom 
(b) and top (t). The following values for the light current quark masses are found from 
the Particle Data group [80], 

m° = 2 to 8 MeV, = 5 to 15 MeV, m° = 100 to 300 MeV. (4.3) 

Notice that the quark masses are renormalization-scheme dependent. The above values 
are obtained in a subtraction scheme at a renormalization scale (9(lGeV). In addition to 
flavour, quarks carry another quantum number known as colour. Each quark comes in 
three colours, red, green and blue. 

The Lagrangian Eq. (4.1) has a large classical symmetry: we have the local gauge 
symmetry SU{Nc) hy construction, 

q ^ U,q, q ^ qU^ U,{x) = exp(ir(x)(y ),), 

A^ ^ U^A^Ul - -U.id^Ul 

We have also global flavour symmetry which does not affect the gluon fields, 

q^Uvq, q^qUl, t/y = exp(^^» (y)^). (4.4) 

where (^)f denotes the generators of the flavor group U{Nf) and Nf denotes the number 
of flavors. The above symmetry is referred to as vector flavor symmetry Uv{Nf). When 
the generator matrix is taken unit matrix we have C/y(l) symmetry associated with con- 
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servation of baryon number. There is another global symmetry which is exact at mP — 0, 
namely chiral symmetry. This symmetry is very similar to vector flavor symmetry, apart 
from an extra factor of 75 in the generator of the transformation. 

q UaQ, q qUA, Ua = exp (^i75^A(y)F) • (4-5) 

Notice that due to the factor 75 the quark field and its conjugate partner are transformed 
by the same matrix in contrast to vector transformation Eq. (4.4). This transformation is 
called the axial-vector transformation and can be combined with the vector transforma- 
tion to define a bigger symmetry at chiral mP — which is then called chiral symmetry 
Uv{Nf) X UA{Nf). One may alternatively define right- and left-handed quark fields by 
following transformation 

1-75 1 + 75 

Ql = — q, qR = ^ g, (4.6) 

The right- and left-handed massless fermions are eigenvalues of the helicity or chirality 
(with eigenvalue ±1) and are not mixed together. The chiral symmetry can be equivalently 
written as x C/r(A^/). 

It is believed that intermediate-energy hadronic physics, say over range of energy 
MeV-GeV is adequately described by the dynamics of the lowest-mass quarks, u and d. 
The overall classical symmetry of the Lagrangian with Nf — 2 becomes 

^t/(iVe)iocal X {SU{2)l X SU{2)r X U{l)v X C/(l)A)global, (4.7) 

Not all above-mentioned symmetries survive quantization. Particles with opposite he- 
licity are related by a parity transformation, therefore in a chirally symmetric world the 
hadrons should come in parity doublets. However, in a real life we do not observe such de- 
generacy. Therefore, one can conclude that chiral symmetry is not realized in the ground 
state and chiral symmetry is spontaneously broken. The Goldstone theorem tell us that 
the spontaneous breaking of a continuous global symmetry implies the existence of as- 
sociated massless spinless particles. This indeed confirmed due to the existence of the 
light pseudoscalar mesons in nature (pions, kaons and etas) as the corresponding Gold- 
stone bosons [81]. Moreover, the existence of quark condensate (qq) implies that the 
SU{Nf)L X SU{Nf)ji symmetry is spontaneously broken down to SU{Nf)v. There- 
fore one may conceive QCD quark condensate as an order parameter for chiral symmetry 
breaking. The concept of spontaneous broken chiral symmetry is the cornerstone in the 
understanding of the low-energy hadronic spectrum. 



86 CHAPTER 4. QCD PROPERTIES AND EFFECTTVE QUARK CHIRAL MODELS 



The U{1)a symmetry implies that all hadrons should come with opposite parity part- 
ners. However, this is not the case, therefore this symmetry must be broken somehow. 
If the spontaneous symmetry breaking mechanism works here, then one should observe 
Goldstone boson associated with U{1)a, namely an / = pseudoscalar meson having 
roughly the same mass as the pion. Surprisingly there is no such Goldstone boson. This 
problem is sometime called U{1)a puzzle. It turned out that the symmetry is ex- 

plicitly broken by quantum effects. This effect is known as the axial anomaly [82]. It was 
shown by 't Hoof that due to instanton effects, the U{1)a symmetry is not manifested in 
nature [82]. 

Finally, at mP — 0, the QCD Lagrangian is invariant under a scale transformation 
which is called dilatational symmetry: 

q{x) ^ e^/^q{e-^x), A^x) ^ eA^e-^x), ^ V (4-8) 

This symmetry is again broken at quantum level due to the trace anomaly [83]. 



4.2 Non-perturbative features of QCD 

In this section we shall recapitulate the most important features of QCD which are not 
accessible perturbatively. Let us firstly elaborate why perturbation theory in terms of 
coupling g can not be used in the low-energy regime of the theory. Having introduced the 

gauge fixing term and an associated ghost term by means of the Faddeev-Popv procedure 
[50, 79, 84], one can carry out perturbation theory in terms of coupling. Due to the 
renormalization process, a renormalization scale fi enters the algebra [85]. Therefore the 
running coupling is described by the RG equation 

^ = (4.9) 

din II 

Where the coupling is small, the P function can be computed perturbatively, 
with 

/3o = ll-^iV/, A = 102-yAr;. (4.11) 
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Therefore, one can readily calculate the effective running coupling 

^ - 9^^^) _ l^TT (918 - lUN,) ln(ln(^VA^)) 

^^^^ 47r (33 - 2Nf) ln(/i7A2) ^ (33 - 2NfY ln(/iVA2) ^' ^ ' ^ 

where A is a scale parameter of QCD and depends on the subtraction scheme and the 
number of active flavours, 

A^ = (208+i)MeV, (4.13) 



where the symbol MS stands for minimal subtraction scheme [79] and the superscripts 
indicate the number of active flavours. This value is taken from an analysis of the various 
high energy processes, see Ref. [86]. The most striking feature of the running coupling is 
that it decreases logarithmically as ji increases. Therefore perturbation theory works very 
well for large /x. This phenomenon is called asymptotic freedom [87]. However, if ji is 
near Aj^g, perturbation theory does not work anymore and non-perturbative phenomena 
enter the stage. Admittedly, there is no unambiguous method available to connect small 
and large distances in QCD. 

One of the most important non-perturbative features of QCD is dynamical chiral sym- 
metry breaking which is responsible for generation of a quark mass from nothing^ In 
order to show that this phenomenon is purely non-perturbative, we employ the QCD gap 
equation [90], 

S{p)-' = {il.p + m') + J ^yD,,{p - q)^^^S{q)r:{p, q), (4.14) 

where and g are the current-quark bare mass and the coupling constant, respectively. 
Dfj.u{p—q) is the dressed-gluon propagator and r°(p, q) is the dressed-quark-gluon vertex. 
The general solution of the gap equation is a dressed-quark propagator of the form 

One may now use the gap equation to work out the fermion self -energy perturbatively 
[91]. One obtains, 

S(p2) = m° (^1 - ^ ln(pV^') + -) ■ (4.16) 
It is observed that at all orders of loop expansion terms are proportional to the current- 



'There is another very different way to generate mass from vacuum, the so-called Casimir effect [88] 
which originates from the response of vacuum in the presence of non-perturbative boundary condition. The 
existence of boundary conditions in quantum field theory is not always free of problems (see e. g., [89]). 
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quark mass and consequently vanish as — > 0. Therefore, no mass (the quark mass 
is defined as a pole of the dressed-quark propagator) is generated at current-quark mass 
equal to zero, i.e. the dynamical chiral symmetry breaking is impossible in perturbation 
theory and there is no mixing between left- and right-handed quarks "perturbatively". 
Notice that apart from the trivial solution -B(p^) = at m = 0, a non-trivial solution 
B{p'^) 7^ can indeed be found at the chiral point, albeit accessible non-perturbatively. 
The renormalization effect is not included in Eq. (4.14), but it does not change the above 
argument [91]. As we already mentioned, there is a close relationship between the gener- 
ation of the quark mass, -B(p^) 7^ 0, and the fact that (qq) 7^ 0. The quark condensate in 
QCD is given by the trace of the full quark propagator Eq. (4.15), 

(qq) = —ilimTiS{x,y). (4.17) 

Notice that since qq is a gauge invariant object, one may take any gauge to obtain the 
dressed quark propagator which has a general form as equation (4.15). It is obvious when 
we have -B(p^) = 0, never does the quark condensate take place, simply because of the 
identity Tr7^ = 0. It has been shown in Landau gauge that the dynamical quark mass, 
M{p^) is large in infrared, M(0) ~ 0.5 GeV, but is power-law suppressed in the ultraviolet 
[92], 

where 7^ = 12/(33 — 2A^/) is the mass anomalous dimension and {qqY is the renor- 
malization group invariant vacuum quark condensate. The dressed-quark mass-function 
Eq. (4.18) is a longstanding prediction of the Dyson-Schwinger equation [90] and has 
been recently confirmed by quenched lattice QCD [93]. It has been shown in many non- 
perturbative approaches that the emergence of a dynamical quark mass leads to the non- 
vanishing of quark condensate and vice versa, e.g., see chapter 6. 

Another important non-perturbative feature of QCD is color confinement [94]. Loosely 
speaking, confinement is defined as the absence of any colored object in nature. But it 
is possible that there exists a composite colored particle which can form colorless bound 
states with another colored particle like quarks. The color confinement is still not prop- 
erly understood, and a clear and indisputable mechanism responsible for this effect yet 
remains to be discovered ^. Confinement originates non-perturbatively, since it is asso- 

r dg 

ciated with a linear potential with a string tension a oc A^e ^ 1^(3) which is obviously 



^The Clay Foundation is offering $1 million prize to anyone who can provide a mathematical proof of 
confinement. 
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non-perturbative in the coupling^. 

One may wonder if there is a non-trivial solution for gap equation B{p'^) =^ which 
gives rise to a pole of the quark propagator, which would contradict QCD confinement 
since the quark is colored Indeed this is one of the subtle point in every QCD model 
and can not be easily resolved. In principle, there will be a long-range force between 
massive quarks to confine them and also a short range spin-spin interaction between mas- 
sive dressed quarks. The former will modify the low momentum part of the propagator 
to remove the quark from being on-shell. Actually, this describes a phenomenologically 
motivated picture of a constituent quark model based on the dynamical symmetry break- 
ing. Having said that it is very hard to incorporate the dynamical symmetry breaking and 
the confinement into a QCD model. In fact, many models constructed to describe the 
low-energy properties of hadrons are assumed to be only dominated by the quark flavor 
dynamics and dynamical symmetry breaking and are indeed reliable only at intermediate 
scales, between confinement scale few hundred MeV up to a scale about 1 GeV. 



4.3 Effective low-energy quark interaction 

Physical hadrons are colorless objects and their properties seem to be determined by the 
flavor dynamics which is induced by an effective QCD interaction. The first step toward a 
construction of such effective theory is to integrate out gluonic degrees of freedom, then 
by standard bozonization and hadronization methods [95] derive the desired effective low- 
energy theories based on relevant degrees of freedom. There have been many attempts to 
attack this difficult problem. Two very well established methods are the global color 
model approach introduced by Cahill and Roberts [96], and the so-called field strength 
approach introduced by Reinhardt and collaborators [97]. Our main goal in this section 
is to give a taste of such approaches and focus only on the main themes rather than de- 
tails. First we rewrite the quark- gluon interaction term in QCD Lagrangian Eq (4.1) by 
rewriting 

qrA,q = A';^J^, J^^q^^q. (4.19) 

^Note that the string picture of quark confinement is not free of flaws, since string breaking will occur 
once the potential energy approaches the quark pair creation threshold. 

"'it is well-known that for confinement it is sufficient that no colored Schwinger function possesses 
a spectral representation. It is equivalent to say that all colored Schwinger functions violate reflection 
positivity [20, 90]. This is one way of reaUzation of QCD confinement. There are in fact many different 
ways that the confinement can be realized such as monopole condensation effect, infrared enhancement of 
the ghost propagator etc. For a review of this subject see Ref.[94]. 
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The full quantum theory of QCD is solved by computing the functional integral describing 
the vacuum-to-vacuum transition amplitude, 

Zqcd = J VqVq J p^«e^ ^ '^'^^'^^^ . (4.20) 

In order to handle the gluonic part of the QCD functional integral, one first has to define 
gauge inequivalent orbits by using a standard gauge fixing procedure and the Faddeev- 
Popov method [50, 79] (in what follows we assume that this procedure has already been 
carried out). We now split the above generating functional integral as 

Zqcd = J VqVqeicp (^i J d^xq{i-t''d^ - m% + r[J]^ , (4.21) 

where the gluon part of action is defined in 

r[J] = log j PA» exp (^-^jF' + gj A^J^^ . (4.22) 

If we could evaluate exactly the gluonic part of the functional integral then we would be 
done. But unfortunately this integration can not be handled unless we resort to some sort 
of systematic approximation (this is due to the presence of cubic and quartic terms of A^ 
in the Lagrangian). One possibility to proceed is to expand the effective action in powers 
of the quark current as suggested by Cahill and Roberts [96], 

r[j] = r[j = o] + ^|r«v„^dxi + ^|r2(xi,x2);:r4J„^^^ 

+ i_ y r(^){x,,x2,..x^)z:XJi:^---J^:dxi...dxr,, (4.23) 



where the coefficients 



rH(.„ X,, xn)-:;; = (Jlf^) , (4-24) 

V^-'MI ■■^'^Hn J J=0 

are defined as one-particle irreducible gluon correlation functions in the absence of quarks. 
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In the above expression, the brackets denotes the functional average over the gluon field 

The zeroth order term T[J = 0] does not depend on the quark field and is therefore an 
irrelevant constant. The first order term T'^^^ gives the expectation value of the gluon field 
and is zero in the absence of external fields. The leading non-trivial term is r^(,Ti, .^2)^1% = 
D'^'l{x — y) which is the exact gluon propagator and includes all gluon self-interactions 
and gluon-ghost interactions. Notice that the quark loops are incorporated through the 
quark current attached as legs to the exact gluon propagator. The main approximation 
in this approach is to ignore all terms n > 3. Note that none of the gluon correlation 
functions is gauge or Lorentz invariant. While each term in the expansion is separately 
invariant under Lorentz and "global color symmetry", nevertheless the whole expansion 
in Eq. (4.23) is invariant under local gauge symmetry. We truncate Eq. (4.23) up to n = 2, 
and with this simplification the QCD generating functional is approximated by 

Zqcd ~ j 'Dq'Dqexjp{iSQFD), (4.26) 

where the Sqfd the induced non-local QCD action which describes the quark flavour 
dynamics, 

Sqfd^ J q{ird^.-m^)q+^ J J q{x)^^.^q{x)D{x - y)q{y)^^^q{y). (4.21) 

The exact form of D{x) is not available at the moment. Hence, the main phenomenol- 
ogy task is to simulate the simplest form of gluon propagator in order to reproduce the 
confinement and asymptotic freedom of quarks. Although there is no a priori reason to be- 
lieve that the effective interaction of quarks propagating in the QCD vacuum should retain 
the form of a one gluon exchange interaction, it has been proved that already this simple 
approximation reproduces most phenomenological models such as the instanton liquid 
model [99], Nambu-Jona-Lasinio (NJL) model [98], various chiral bag or topological- 
soliton models [96] and the Dyson-Schwinger equation approximation [90], etc. As an 
example, a useful starting point for low energy effective action is to employ a gluon prop- 
agator which is reduced to its crudest form, 

X 
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where is a constant with dimension of (energy)^. This choice leads to a local NJL type 
model with interaction: 



This interaction describes a system of quarks interacting via a two-body-force. The lo- 
cal form of the interaction of course causes ultraviolet divergences, which introduces an 
energy scale A^^, breaking the scale invariance of the classical Yang-Mills Lagrangian (at 
zero current quark mass), in an anomalous fashion. We have sketched how an effective 
quark theory can be approximately obtained from QCD by eliminating gluonic degrees of 
freedom. This is slightly different from the RG approach discussed in the first part of the 
thesis in which our main goal was to eliminate the high-energy degrees of freedom. Hav- 
ing said that, both have the same foundation, namely eliminating the irrelevant degrees of 
freedom and as we already discussed can be implemented at the same time. 

Unfortunately, there is no economic way to derive effective theories from QCD for 
a given system. Therefore, one may write down the most general Lagrangian based on 
relevant degrees of freedom, having imposed some general constraints such as symmetry 
properties. This approach was introduced by Weinberg [100] and later by Gasser and 
Leutwyler [101]. For an example, an effective chiral quark theory can be presented as 



where C„ are the local chiral-invariant operators consisting of quark fields and c„ are 
dimensionless coupling constants. The theory is only valid below the scale A, and the 
momenta of the loop integrals are cut off at A. One may now truncate the above series 
and by first obtaining the coupling constant through experimental input, calculate other 
quantities. This EFT approach has been discussed in detail in Ref. [101]. 



4.4 Fierz-transformation and the effective quark interac- 
tion 




(4.29) 




n 



(4.30) 



The NJL model Lagrangian Eq. (4.29) contains the color octet flavor singlet currents of 
quark. Since hadrons are color singlets, it is desirable to recast the Lagrangian in such 
to act in color singlet channels. This can be accomplished by a Fierz transformation, see 
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e.g., Ref. [98, 102], using the relation 

If we take the interaction in qq channel, then the first part of the above expression creates 
color singlet mesons, while the second part is color octet. However, the signs behind the 
color singlet and the color octet are opposite, therefore if we choose the negative sign for 
color singlet (in the Lagrangian level), i.e. an attractive interaction, then the interaction 
for color octet will be repulsive and consequently no bound state exist for this channel 
consistent with nature. Moreover, at large Nc the color octet qq can be neglected. In 
the same fashion one can recast the interaction in the qq channel. In order to make a 
color- singlet baryon out of three quarks, we first couple two quarks in the fundamental 
triplet representation 3c which leads to either a sextet 6c or an antitriplet 3c. But only an 
antitriplet can be combined with another quark in 3c to make a color singlet state. Notice 
as well that 3c is antisymmetric and 6c is symmetric in color quantum numbers. We can 
now use the following Fierz transformation, 

(y) . . (y) ,r " ^) ^ w/-'''-^- ^^-^^^ 

It is obvious that the interaction in the qq 3c channel becomes attractive, therefore diquark 
formation in 3c is in principle possible. In the same way, one can "Fierz" the flavor and 
Dirac quantum numbers. For the meson channel we use 




(4.33) 



Notice that in contrast to color group SU{Nc) the flavor group is U (Nf), hence the flavor 
generators includes A°/2 = 1 / y^2Nf. One may immediately read off from the above 
decomposition that for Nj = 3, mesons occur as nonets under flavor transformation. For 
the diquark channel we use 

where A™ denotes the symmetric generators of U{Nf ), i. e. A"^, — A^, and A^. The anti- 
symmetric part A^' stands for A°'^'^'^'^'^ . Finally the Dirac indices can be rearranged for 
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meson channels by making use of 



{lix)ijWki = ^iAj + {ilb)ii{ilb)k3 - 2 {{lij)iiWk3 + {lnlb)ii{Ylb)kj) ■ (4.35) 



For the diquark channel we employ 

{lu)ij{l^)ki = CikCij + {n5C)ik{Ci-f5)ij 



1 



(4.36) 

where C = ^7^7° is the charge conjugate matrix. Using all the above transformations, 
one can now rewrite the gluon exchange Lagrangian Eq. (4.29) (for Nf — Nc — 3) in the 
following compact form 



9^ 



2AI 



Cf + £f , 



3A^ 



(4.37) 



where the charge conjugate spinor spinor is given by g*^ = Cq^ and the other notations 
are defined as follows 



Ma 
E„ 



2 



A = 0,...8, 
a = 1,2, 3; 4e{KAs}, 



(4.38) 



We will study with full detail the mesons, diquarks and the baryons structures within a 
non-local version of the above Lagrangian in the chapter 6. 

Here some remarks are in order. First of all, the Fierz transformation does not spoil 
chiral symmetry. The meson channel Cf and diquark Cf are separately chirally invari- 
ant. The main difference between mesonic and diquark interaction stems from the Pauli 
principle. The diquark qq vertices are antisymmetric in color indices, hence the Pauli prin- 
ciple requires a certain combinations of Dirac and flavor vertices. Notice that for Nc 7^ 3, 
one has to replace the factor 1/3 with [Nc — l)/2Nc and l/Nc for meson and diquark 
channels, respectively. This indicates that in large Nc limit, the colored diquark channels 
are suppressed by a 1/A^c factor and only color singlet mesons and baryons survive. This 
is in accordance with Witten's conjecture that in the large- A^c limit, QCD transforms into 
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a theory of weakly interacting mesons and baryons emerging as solitons of this theory 
[103]. We will elaborate on this conjecture in the next chapter. 
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Chapter 5 

QCD inspired pictures for baryons 



5.1 Introduction 

Despite all efforts to describe hadron physics in terms of its underlying QCD theory, a 
unified and unambiguous description of hadrons is still missing. This is due to the com- 
plexities and the non-perturbative features of low-energy sector of QCD which prohibits a 
straightforward computation of hadronic properties in this regime. One well-established 
method toward understanding the hadronic physics is the use of QCD sum rules [104], 
which aim to interpolate between the calculable high-energy behaviour of the QCD and 
low-energy phenomenology. Although this approach has been reasonably successful phe- 
nomenologically, there are many uncertainties induced by the choice and formulation of 
the phenomenological part, and the result can exhibit a significant dependence on the 
mass scale at which the matching is performed. More importantly, the confinement phe- 
nomenon and instanton effects are not incorporated in this formalism. Another method 
is to simulate QCD on a lattice of spacetime points. Some appreciable progress has been 
made in ab initio calculations of low-lying baryon resonances using lattice QCD simula- 
tion on a computer. This is limited by computational resources. This leads to pion masses 
of usually more than 400 MeV (thus, the chiral point is not accessible in lattice QCD) or 
to simulations within the quenched approximation, where sea quark effects are neglected. 
Moreover, there is yet incomplete understanding of systematic errors, e.g., finite-size ef- 
fects ^ 

'For understanding QCD phase structure, one should describe the important non-perturbative nature 
of QCD and the hot/dense QCD in a unified way. Thus far, lattice regularization have not be able to 
overcome the issue associated with the fact that at non-zero chemical potential the fermionic determinant 
is complex. On the other hand due to computational difficulties, a single lattice of equally spaced points is 
forced to span all distance scale which does not allow a sequence of descriptions intermediate between the 
constituent quark model and QCD. Full consideration of QCD phase structure is not possible. Nevertheless, 
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Therefore, any models which simplify the QCD dynamics and is capable of describing 
part of the hadronic physics (which might not be accessible to ab initio methods, like 
lattice QCD) is greatly appreciated. 

Our main goal in the following is to study baryon pictures in a relativistic framework 
by employing only quark degrees of freedom. There are two distinct possibilities to build 
a model for baryons, one is in the limit of a large number of colors based on the picture of 
baryons as soliton, the second one is to describe baryons in term of bound state of a finite 
number of colors. The key ingredient of the former will be introduced in section 5.2, and 
the latter will be described in sections 5.3,4. 

5.2 Baryons as chiral solitons; Skyrme model 

The idea to describe baryons as solitons was first introduced by Skyrme [105] before the 
advent of the quarks and gluons. However, this was not fully appreciated until the con- 
jecture of Witten [103]. It was 't Hooft's proposal that the inverse of the number of colors 
1/iVc, can be treated as an effective expansion parameter [106]. Later, Witten pursued 
this idea and showed that QCD is reduced to an effective theory of weakly interacting 
mesons (with an effective four-meson vertex scaling like 1/Nc), and that baryons emerge 
as soliton solutions of this meson field without any further reference to their quark content 
[103]. Since then, this subject has been well studied. Here, we intended to explore this 
approach in a selective way, following Ref. [107] closely . 

At low-energy one expects that the effective meson theory is a type of non-linear a 
model with pions as the lighter mesons. We shall focus on the massless two flavor model. 
One can combine the four fields (tt, a) into one unitary 2x2 matrix U (chiral field) with 
only one isovector field (p{x, t). 



where the isovector r contains the pauli matrices. In terms of U {x) the non-linear a 
model is defined by the Lagrangian 



where c is a constant. The elements of chiral SU{2) x SU{2) transform any chiral field 

the effective field theories approach based on simple models embodied essential ingredients of QCD at low- 
energy to give a qualitative (and even sometimes quantitative) prescription. 




(5.1) 




(5.2) 



5.2. BARYONS AS CHIRAL SOLITONS; SKYRME MODEL 



99 



as 

U{x) LU{x)R\ (5.3) 

where R and L are arbitrary SU (2) matrices. The Lagrangian is invariant under such 
transformations. The vacuum configuration = i.e. [/ = 1) is only invariant under the 
coset L = R which reflects the spontaneous breaking of chiral symmetry. One can readily 
construct the Noether current associated with the symmetry transformation Eq. (5.3). The 
vector and axial-vector current correspond to R = L and R^ = L, respectively. The 
matrix element of the axial-vector current between the vacuum and a one pion state relates 
the unknown coefficient c in Eq. (5.2) to the pion decay constant. It turns out that c = 
U = 93 MeV [107, 108]. 



Now, we want to find a solitonic solution for our meson theory. A fundamental re- 
quirement for a solitonic solution of the equation of motion is finiteness of energy. The 
static soliton configuration U{r) represents mappings U : TZ'^ SU{2). A necessary 
condition is to require the boundary condition 

lim U(r) = 1. (5.4) 

r— >cxD 

This means that spatial infinity is mapped to one point in flavor space, i.e. TZ^ is compact- 
ified to the hyper sphere S^. 

U -.S^ ^ S\ (5.5) 

which induces a topological invariant, the winding numbers. The associated topological 
current has been given by Skyrme [105] 

= ^e''''^'Tr[{U^dM){U^d,U){U^d^U)], (5.6) 

with conservation law d^B^ = 0. The integral over the zero component defines the 
topological charge 




(5.7) 



which as we will show, introduces the winding number n. A static configuration of the 
soliton is given by the spherically symmetric hedgehog Ansatz 

C/(r) =exp(iT.fe(r)). (5.8) 



In order to understand the geometrical meaning of the above Ansatz, we use this expres- 
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sion to compute the topological charge from Eq. (5.7). 

B = ^2'''" J Tr[{U^diU){U^d,U){U^d,U)]d'x, 

^^drO £x = - / sin2(e)de = n, (5.9) 
Je(o) 

where we imposed the boundary conditions 0(0) = —mi and Q{oo) = 0. Solitons with 
different winding numbers are topologically distinct, and there thus exists no continuous 
deformation connecting solitons of different winding number. It has been conjectured by 
Skyrme [105] that the topological current can be related to the baryon current and the 
winding number with the baryon number. 

One may scale the spatial coordinate r in U{r) by Ar i.e. U — > U{Xr), this leads to 
scaling the potential part of Lagrangian Eq. (5.2): C^'^^ — > jC^^^. Therefore the minimal 
energy is only obtained for A — > oo, e. i., no stable solitons can be found and the soliton 
collapse to zero size. However, if the Lagrangian contains a term containing products of 
four spatial derivatives (but only quadratic in time derivative for sake of quantization), it 
will scale as C^'^^ — > jC^^^ + Ai2^^) which stabilizes at A = C^'^yC^'^\ A possible form of 

£(4) 

is given by 

C^'^ = ^ {Tr[{U^d,U), {U^dMMU^d^U), (U^d'^U)]) , (5.10) 

where e is an extra parameter which determines the size of the particle. 

A natural question is how baryons get their half integer spin and isospin within soliton 
picture since the pion field possesses spin zero and isospin one. An immediate answer is 
quantization. We refrain to go through details here and concentrate only on main points. 
The time dependent soliton solution should be obtained as a first step toward quantization. 
However, such solutions are hardly available and one needs to resort to an approximation. 
A reasonable Ansatz for such a time-dependent solution is given by [107] 

U{r,t)^A{t)Uo{r)A\t), (5.11) 

where A{t) E SU{2) is referred to as the collective rotation and contains collective co- 
ordinates (it can be parametrized in terms of Euler angles). This Ansatz does not change 
the potential energy of the hedgehog. The time derivative in the Lagrangian, produces 
terms which represent the rotational energy of a rotating skyrmions. Having used hedge- 
hog properties and a suitable definition of angular velocities as canonical variables, one 
can show that the absolute values of spin and isospin are equal, = J^, and the energy 
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eigenvalues of the system forms a rotational spectrum Ej — ^J{J + 1) + M, where M 
is the static energy of the soliton and 9 denotes the moment of inertia. One can imme- 
diately observe that the quantum numbers of the low energy baryons (e.g. the nucleon 
with J = T = 1/2 and delta J = T = 3/2)are consistent with this picture. One of 
the obvious shortcomings of this presentation is that it does not give any reason in favor 
or against half-integer or integer values of J for the quantized skyrmions. We will argue 
later that the fermionic character of the quantized baryon can be revealed in a model with 
SU (3) X SU (3) symmetry by inclusion of the Wess-Zumino-Witten term. 

An extension of the theory to three flavors with chiral SU (3) x SU (3) symmetry is 
essential if one wants to consider the baryon octet and decuplet, see for details Ref. [107, 
108]. We assume the generic form of the underlying Lagrangian remains unchanged 
C = jC'--'^^ + C'^'^K However, the chiral field U needs to be increased to a SU (3) field with 
the mesons fields 0", a = 1, .., 8, which in addition to the pions, contains the kaons and 
the octet component of the rj 



where A" denotes the Gell-Mann matrices. The decay constants /„ are defined through 
the gradient expansion of the axial-vector current, analogous to the case with only two 
flavors. In a similar fashion to the two flavor case, one can obtain solitonic solutions of 
the theory and quantization can be done by introducing the collective rotations. 

In order to link the effective meson theory to QCD, one should firstly consider if all 
symmetries of the Lagrangian C are in accordance with QCD. Witten [109] observed that 
the Lagrangian Eq. (5.2) possesses an extra discrete symmetry that is not a symmetry of 
QCD. Under parity transformation P, the pseudoscalar meson fields as described by QCD 
should obey Pn{x, t) = — 7r(— x, t). In our meson theory this means. 



But it is obvious that the Lagrangian Eq. (5.2) is invariant under x —x and U U\ 
separately. In order to break this unwanted symmetry, one needs to add some extra term 
to the meson action. Unfortunately, there exists no local term in four spacetime dimension 
which can be added to the Lagrangian, so as to get rid of this separate symmetry. However, 
it is very easy to look for such extra term by considering the equation of motion. Witten 
[109] suggested that the simplest term (with lowest possible number of derivatives) which 




(5.12) 



P : X 



-X, t^t, U ^ UK 



(5.13) 
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needs to be added are as follows 

S-di^C^ + Xe^^'P'^C^CCpC, = 0, (5.14) 

o 

where = U^d^U. The new term is odd under x — > —x while the first and sim- 
ilar higher-order terms are even. However, the new term is even under transformation 
U ^ W while the first term is odd. Therefore, Eq. (5.14) is invariant only under the 
combined action of P Eq. (5.13). The problem now is that the four-dimensional action 
corresponding to the new extra term can not be written in a chirally invariant form. How- 
ever, this action can be rewritten in such a form in five dimensions. Therefore, we extend 
the coordinate space to five-dimensional manifold M5 in such a way that our conven- 
tional four-dimensional spacetime M4 is the boundary of a five-dimensional volume i.e. 
dM^ = M4. Therefore, one can write 

r^xf e'^''^'^TT{CiCjCkCiCm)d^x. (5.15) 

This action indeed leads to the equations of motion Eq. (5.14) where written in four- 
dimensional spacetime by means of Stokes 'theorem. Witten in his remarkable paper 
[109], argued that the coefficient A in above equation must be integer multiple of a nor- 
malization factor, 

A = »^. (5.16) 

This is comprehensible, since the path-integral formulation requires the action to be 
changed by a multiple of 27r when going from M5 to its complement, which has the 
identical boundary with opposite orientation. This is indeed in the same spirit of Dirac's 
quantization of a magnetic monopole. The physical meaning of the integer n is fixed 
through a connection to the Wess-Zumino action ^ [110]. Witten included the photon 
fields in a gauge invariant way which generates a vertex for the decay 7r° — > 77 with 

-7^°6„..,F^'^F^^ (5.17) 



where F^*^ is the field strength tensor of the photon. One can immediately compare this 
result with the well-known triangle anomaly in QCD and find that n — Nc- In this way, 
the effects of anomalies in QCD are correctly reproduced by the action Eq. (5.15) which 
is called Wess-Zumino-Witten (WZW) term. Notice that if one considers an adiabatic 27r 



^The Wess-Zumino action was introduced to account for anomalies which occur through the renormal- 
ization of fermion loops in quantum field theories where pseudoscalar mesons are coupled to fermions 
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rotation of the soliton, the WZW term produces a contribution NcTt to the action while 
other terms do not contribute. Therefore, the soliton acquires a phase (—1)^° for such a 
rotation as required for fermions with Nc odd or bosons with Nc even. It is interesting that 
for two flavours, this conclusion can not be made since the WZW term is zero. There- 
fore, for three flavors, the WZW term provides some hint about the statistical difference 
between baryons and mesons (fermions or bosons). The fact that the Skyrme Lagrangian 
needs to be augmented by the WZW term indicates that the underlying physics of the 
Skyrme is a model of quarks and gluons which possess QCD properties. 

As we discused in the last chapter, one of the simplest but viable quark model at the 
present is the NJL model. The solitonic solutions of NJL models have been extensively 
investigated [108]. Notice that, in contrast to Skyrme-like models with infinite energy 
barriers separating sectors with different winding numbers, chiral quark models, such as 
the NJL models have finite energy barriers separating the different sectors, and they give 
rise to so-called non-topological solitons. 

Despite all appealing features of soliton models, there are some shortcomings: It is 
well-known that solitonic models are not very accurate, e.g., in leading order of Nc, the 
quasi-classical soliton configuration with quantized collective variables can produce bary- 
onic observables with errors of about 20% — 30%, and corrections due to mesonic fluctu- 
ations seem to be very important. Another point of weakness is that exotic states are very 
controversial in these models^. 

5.3 Bag models 

In 1974 MIT group [86] developed a new picture of hadrons based on the simple as- 
sumption that the physical vacuum prohibits free quarks and gluons, but instead creates 
bubbles of hadronic size in which quarks and gluons may propagate ordinarily. This idea 
has been employed in various models: a hybrid bag model where the nucleon consists of 
a quark bag surrounded by a meson cloud, the little bag model where in contrast to the 
hybrid bag, pions are not allowed to propagate inside the bag, a cloudy bag model where 
the mesons are constrained to the chiral circle and are allowed inside bag, and finally the 
chiral bag model where the constrained mesons outside bag are described by the Skyrme 
Lagrangian. In all these models, there is one extra parameter in the model, the bag radius, 

■'As Cohen [1 1 1] argued the main reason is that the rigid-rotor quantization is not vahd for such states. 
In other words, the assumption that the collective motion is orthogonal to vibrational motion is only true 
for non-exotic motion, but the Wess-Zumino term induces mixing at leading order between collective and 
vibrational motion with exotic quantum numbers. Recent discovery of Pentaquark ^+ which was already 
predicted based on soliton model have brought a lot of activity on this subject. 
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which provides some hint about the quark and pion distributions. A review of various bag 
models can be found in Ref. [113]. 

Chiral bag models seem to interpolate between two different aspects of QCD, the long 
range-low energy (non-perturbative) and the small distance (perturbative) behaviours. 
This idea was proposed by Rho et al [114]. The very small volume V represents the 
perturbative domain of QCD containing quarks and gluons only, as opposed to its com- 
plementary piece containing the confined phase of QCD with color- singlet objects such 
as mesons, 

^ — ^Ug^V + '^meson(l " ^v) + '^boundary (^F; 
boundary = -qexp {--f5X.TQ{r)) q, 

where inside the bag we have massless quark fields q and outside we have chiral meson 
fields U which obeys the Skyrme Lagrangian (for simplicity the hedgehog configurations 
U = exp {iT.x<d{r)) with spherical bags of radius R are assumed ). The boundary term 
cause the full Lagrangian to be invariant under a combined chiral symmetry; 

q' ^ exp {i^^a.T) q, C/' — > exp (— ict.r) C/ exp (— ia.r) . (5.18) 

We already identified the topological charge carried by the meson field as baryon number. 
On the other hand, the quarks inside the bag each carry one third amount of baryon charge, 
therefore, it seems puzzling that the total baryon number is not integral. Goldstone and 
Jaffe [1 15] proved that a conjecture of Rho et al [1 14] that baryon number in the hybrid 
model remains one is indeed right. The crucial observation they made is that the charge of 
the vacuum baryon number, inside the bag is changed due to boundary effects. This shows 
that the baryon number remains one regardless of the profile of the Skyrme fields and the 
size of bag. It was later proved that the total energy from bag, chiral fields and vacuum 
(Casimir energy) is insensitive to the bag size as well [116]. Therefore, it is tantalising 
to assume that observables should not depend on the details of the bag (e.g., its size) 
also. This statement is called "The Cheshire Cat Principle (CCP)" [117]. Topological 
quantities, like the baryon number satisfy an exact CCP in (3+1) dimensions while non- 
topological observables such as masses, static properties and also non-static properties 
satisfy it approximately well [117]. 

The main problem with various bag models is that they are not fully covariant and 
possible large modification of observables due to quantum fluctuations make model pre- 
diction less reliable. 
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5.4 Diquark-quark picture; Relativistic Faddeev approach 

In the previous sections, we reviewed the basic foundations of two well established pic- 
ture of baryons, motivated from QCD properties. Two main shortcomings within these 
approaches, namely the uncertainty of computed quantities and lack of covariance, make 
these not viable for phenomenological usage in the intermediate-energy regime where a 
fully covariant formulation is required. A new generation of continuous beam facilities 
such as CEBAF at TJNAF, ELSA, COSY, MAMI, etc, which are designed to explore the 
intermediate energy between non-perturbative and perturbative regime of QCD, needs ac- 
curate covariant formulation to describe the forthcoming data. The diquark-quark picture 
of baryons based on a relativistic Faddeev approach, is a framework for such a fully co- 
variant approach. This approach has been extensively employed during the last decade, 
and is phenomenologically very successful [118, 119, 120, 121, 122, 123, 124, 125]. In 
this picture, a bound state of a baryon ^ is obtained as a pole of the three-quark correla- 
tion by summing over infinitely many interaction graphs. This process is very similar to 
obtaining a two-body bound state which leads to the Bethe-Salpeter equation. 

In the context of local quantum field theory, the few-body problem seems to be ill- 
defined, since any restriction on degrees of freedom (e.g., particle numbers) may spoil 
the Lorentz invariance (e.g., on the equal-time quantization, the boosts generators involve 
interactions and change the number of particles, therefore limiting the number of particles 
is against the Lorentz invariance). Nevertheless, we know from a phenomenological point 
of view that a fixed number of constituent quarks might be enough to describe baryons. 
Therefore, we introduce the notion of baryon wave functions as matrix elements of three 
quark operators between the physical vacuum |Q) and a (nucleon) bound state \Pn)', ^ ~ 
(Q \T{qqq) \ Pjv) . Having said that, QCD vacuum is non-perturbative and indeed possesses 
non-trivial condensates, and thus the wave function contains sea quark and gluonic parts. 
However, it may be reasonable to ignore all other operator matrix elements which involve 
an arbitrary number of quarks and gluons as irrelevant in comparison to the dominant 
amplitude ^. One can now solve the three-body problem by means of the Green's function 
formulation of quantum field theory. As we will show, the non-perturbative feature of the 
vacuum can then be effectively incorporated into the formalism in a systematic fashion. 

In the following we introduce an approximation scheme based on the relativistic Fad- 
deev approach to simplify the three-quark problem in form of a diquark and quark inter- 

^There is another very old-fashioned approach to obtain baryon bound state which was invented in 
the early sixties, the constituent quark models or quark potential models. In these models one starts with 
simple potential e.g, the hyperfine type interaction and employs 3-particle Schrodinger (or Dirac ) equation 
to obtain the spectrum. This approach is not covariant and does not incorporate the minimal field theoretical 
aspect of QCD, hke quarks degree of freedom. We refrain from discussing this approach here. 
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acting via quark exchange. For simplicity, we use a formal presentation; all color, flavor 
and Dirac indices are implicit in the single particle labels. We use a Euclidean metric in 
momentum space. We denote a dressed single quark propagator by Si, with 

{27Ty6\ki-pi)S,{k,;pi) = J dXd^Xp,e^('=-^'=>-^'^*>)(0|Tg,(a;,)g,(a;,)|0). (5.19) 

In our definition of Green's functions and bound state matrix elements, we always take 
out one 5-function corresponding to conservation of energy-momentum. We define the 
full quark six-point function (or the three-quark correlation function) as 



1=1 



The three-quark correlation function satisfies the Dyson equation, 

G ^Go + Go®K®G, (5.21) 

where Go denotes the disconnected three dressed quark propagator and K stands for the 
three-quark scattering kernel containing all two and three-body irreducible diagrams. The 
symbol "®" denotes summation and integration over all internal and dummy indices. 
A bound state of mass M with wave function ^ emerges as a pole of the three-quark 
correlation function, 

^{ki,k2,k3)-^{pi,P2,P3) 

G{h:Pi) p2^^2 ' (5.22) 

where P = P1+P2+P3 and we defined ^ as a three-body wave function which represents 
the transition matrix element between the vacuum and a bound state with mass M, 

(5.23) 

We now substitute the bound state parametrization Eq. (5.22) into Eq. (5.21) and compare 
the residues. This leads to the homogeneous bound state equation, 

^ = Go ® /s: ® < — > G"^ ® ^ = 0. (5.24) 



5.4. DIQUARK-QUARK PICTURE; RELATIVISTIC FADDEEV APPROACH 107 



Solving this equation exactly is almost impossible since neither the detail of all two- and 
three-particle irreducible graphs appearing in K, nor the full dressed quark propagator 
contained in Gq are known. It is well known that the problem becomes more tractable if 
one employs the Faddeev approximation, by discarding all three-body irreducible graphs 
from the interaction kernel K. In this way one can write the kernel as a sum of three 
two-body interaction kernels, 

K^Ki + K2 + K3, (5.25) 

where Ki with i = 1, 2, 3 refers to the interactions of quark pairs (jk) with a spectator 
quark (i). The two-quark propagators gi satisfy their own Dyson equation with kernel K^, 

gi^Go + Go0Ki0 Qi, (5.26) 

where gi and Ki are defined in three-body space, since Go is defined in three-body space. 
Hence, the former contains a factor Si (the propagator of the spectator quark), and the 
latter contains a factor S^^^ i.e., Ki = k^q S[^. One may associate a disconnected 
scattering amplitude to every spectator quark (i), i.e. Ti = ti® S^^, where tj describes 
the scattering between the quarks (j) and (k) in two-quark subspace. The matrix Ti is 
obtained by amputating all incoming and outgoing quark legs from the connected part of 
9i, 

gi^Go + Go®Ti® Go. (5.27) 
By combining the two previous equations, the Dyson equation for Tj can be found as 

Ti^Ki + Ki®Go® Ti. (5.28) 

Now we define Faddeev components via Eqs. (5.24,5.25) 

^i^Go^Ki® (5.29) 

where we have ^ = X] ^i- We rewrite the Eqs. (5.27) as gi (g) Gq^ = I + Gq^Ti and plug 
this expression into Eq. (5.29) and make use of Eq. (5.26). We then find the well-known 
Faddeev bound state equations, 

"^i^Go^Ti® (jjjj + *fe) = SjSk ® + *fe) . (5.30) 

We have shown that the complicated three-quark problem can be systematically simpli- 
fied by employing the full two-quark correlation function ti, instead of the kernel K. In 
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this way the eight dimensional Eq. (5.24) is reduced to a set of coupled four-dimensional 
equations (5.30). A further simplification of this equation can be achieved by approxi- 
mating the full two-quark correlation ti as a sum of separable correlations, 

ti{ki, h;pi,P2) = xKki, h)Dt{h + k2)x'}{pi,P2), (5.31) 

a 

where the function Xi is the vertex function of two-quark with a diquark and Xj denotes 
its complex conjugate. The index "a" denotes the different channels of the diquarks, and 
Df is the corresponding diquark propagator. Note that separability implies that U does not 
depend on any of the scalar products ki.pj. The diquarks parametrization to some extend 
contains the unknown non-perturbative physics within the baryon structure. A natural 
Ansatz for the Faddeev component is given by 

'^f'{p^,PJ,Pk) = J2^r'SfS^XtMP^^Pk)Dt{p,+Pkma,{P^,PJ+Pk), (5.32) 
a 

where summation over repeated indices is assumed. The Greek multi-indices a, (3,.. de- 
note color, flavor and Dirac indices, and i,j,k indicate the type of quarks. The quark 
(i) and diquark labels (jk) are fixed. The quantity $ is the baryon-quark-diquark ver- 
tex function and depends only on the relative momentum between the momentum of the 
spectator quark, pi and the momentum of the diquark quasi-particle, Pj + Pk. In a rela- 
tivistic formulation of a few-body system there is no unique definition of the momentum 
transfered between individual particles. Therefore, we introduce a new parameter r] which 
parametrises this ambiguity and shows the distribution of the total momentum within the 
system (diquark and quark). We define a relative momentum between the quark (i) and 
the diquark consisting of the quarks (jk) by 

p = (1 - ri)pi - ri(pj + pk) ^Pi- vP, (5.33) 

where P = pi + P2 + Ps- The physical properties of the baryons will indeed not depend 
on T]. One can employ the definition of (5.32) to rewrite the Faddeev equation (5.30) in 
terms of a vertex function $, 

- T.inP.Sr'xi,'a][D'sf^%,] + U ^ k). (5.34) 

b 

In the above derivation we assumed that the quark-diquark vertex function is antisymmet- 
ric under the exchange of the quark labels (x"^^ = ~Xi,-yi3) a consequence of the Pauli 
exclusion principle. Eq. (5.34) resembles a Bethe-Salpeter equation; the first bracket can 
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be conceived as a quark-diquark interaction kernel since it contains the exchange of a sin- 
gle quark between a quark and a diquark, the second term couples the interaction kernel 
to baryon-quark-diquark vertex via diquarks and quarks propagator. For identical quarks 




Figure 5.1: The coupled Bethe-Salpeter equations for the vertex function $. 



the antisymmetrization of the vertex functions is essential, we can now sum over the type 
of particle and drop the index like (i). Therefore, equation (5.34) can be rewritten as, 

Go,yfs{k,P) = 2S^p{k,)D\kd), 
Ka^{k,p,P) = x'p'-y{k,q)Sp'Yxi-y'{Q,Pq), (5.35) 

where we define the spectator quark and diquark momenta as kq{pq) = rjP + k{p), 
kd{Pd) — — v)P ~ k{p), respectively. Momentum conservation fixes the momen- 
tum of exchanged quark q — —k — p + {1 — 2r])P, see fig. 5.1. The factor of two in 
the above presentation originates from the summation over the type of quarks and can be 
absorbed into the definition of diquark propagator. 

In conclusion, we managed to recast the three-body Faddeev equation in the form 
of an effective two-body Bethe-Salpeter equation between the diquarks and the quarks, 
having summed over the ladder-type quark exchange diagrams between the quarks and 
the diquarks. The only assumptions we have made are; 1) we neglect all three-particle 
irreducible contributions 2) we model the connected two-body correlation as a sum of 
separable terms which are identified by diquark channels, see Eq. (5.31). In order to 
obtain an equation for physical baryons, we have to project Eq. (5.34) onto the baryon 
quantum numbers. This will be carried out for a model in the next section. 
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Chapter 6 



Baryons structure in a non-local quark 
confining N JL model 

6.1 Introduction 

The NJL model is a successful (low-energy) phenomenological field theory inspired by 
QCD [98]. The model is constructed to obey the basic symmetries of QCD in the quark 
sector, but unlike the case of low-energy QCD, quarks are not confined. The basic ingre- 
dient of the model, apart from the standard bilinear Lagrangian in the quark fields, is a 
zero-range interaction containing four fermion fields. This means that the model is not 
renormalizable. If we make the standard one over the number of colours {\/Nc) loop 
expansion, already at one-loop level an ultraviolet subtraction (usually implemented by 
a cut-off) supplemented with a regularisation method is required. The value of the cut- 
off can be related to the scale of physical processes not included in the model, and thus 
determines its range of validity. Consequently, processes involving a large momentum 
transfer, such as anomalous decay, can not be described by the model. At higher orders in 
the loop expansion, which are necessary for calculating mesonic (baryonic) fluctuations 
[126, 127], one needs extra cut-off parameters. It is hard to determine these parameters 
from independent physics, and thus to build a viable phenomenology. A similar problem 
appears in the diquark-quark picture of baryons where an additional cut-off parameter is 
required to regularise the diquark-quark loops [127]. It has been shown that a renormal- 
izable extension of the NJL model (at least to one-loop level) [128] can be constructed 
by matching the NJL-contact interaction at low energy with a one-gluon exchange type 
interaction above the Landau pole. Another way to cope with the non-renormalisability of 
the model is to embed this model into a renormalizable theory such as the linear cr-model 
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and apply a renormalisation group approach [129, 130]. However, all these approaches 
add extra complexity which make them difficult to apply for anything but very simple 
problems. 

Another drawback of the model is the absence of confinement, which makes it ques- 
tionable for the description of few-quark states and for quark matter. If energetically 
allowed, the mesons of the model can decay into free quark-antiquark pairs, and the pres- 
ence unphysical channel is another limitation of the applicability of NJL model. At the 
same time, it is also known that the NJL model exhibits a zero-temperature phase transi- 
tion at unrealistically low baryon density [102]. This problem is caused by the formation 
of unphysical coloured diquark states. These may be explicitly excluded at zero density 
by a projection onto the physical channels, but dominate the behaviour at finite density. 
The model is not able to describe nuclear matter, even in the low-density regime [131]. 

We do not know how to implement colour confinement in the model and, anyway, 
the exact confining mechanism of QCD is still unknown. In the context of an effective 
quark theory, a slightly different mechanism of "quark confinement" can be described by 
a quark propagator which vanishes due to infra-red singularities [132] or when it does not 
produce any poles corresponding to asymptotic quark states [133, 90]. Another realisation 
of quark confinement can be found in Ref. [134]. It has been shown that a non-local 
covariant extension of the NJL model inspired by the instanton liquid model [99] can 
lead to quark confinement for acceptable values of the parameters [135]. Here the quark 
propagator has no real pole, and consequently quarks do not appear as asymptotic states. 
The quark propagator has infinitely many pairs of complex poles corresponding to quarks 
which have a finite lifetime. This phenomenon was also noticed in Schwinger-Dyson 
equation studies in QED and QCD [136, 137, 138]. 

We can simply accept the appearance of these poles as an artifact of the naive trunca- 
tion scheme involved. However, it has been recently suggested that it might be a genuine 
feature of the full theory, and be connected with the underlying confinement mechanism 
[137, 138]. For example, it has been shown by Maris that if we remove the confining 
potential in QED in 2-i-lD the mass singularities are located almost on the time axis, and 
if there is a confining potential, the mass-like singularities move from the time axis to 
complex momenta [138]. In this chapter, we study this kind of confinement from another 
viewpoint. We show that when we have quark confinement in the non-local NJL model, 
the baryons become more compact, compared to a situation where we have only real poles 
for quark propagator. 

There are several other advantages of the non-local version of the model over the 
local NJL model: the dynamical quark mass is momentum-dependent and also found 
in lattice simulations of QCD [139]. More importantly, the non-locality regularises the 
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model preserving anomalies [135], and the regulator makes the theory finite to all orders 
in the 1/Nc expansion, and leads to small next-to-leading order corrections [140]. As a 
result, the non-local version of the NJL model may have more predictive power. 

The instanton-liquid model is only one way to motivate such a model [99]. Many 
effective field theories constructed by the Wilsonian renormalisation group approach lead 
to a non-locality, at least as irrelevant terms in the renormalisation group sense. Non- 
locality also emerges naturally in the Schwinger-Dyson resummation [90]. Considerable 
work has been done on these nonlocal NJL models including applications to the mesonic 
sector [135, 141], phase transitions at finite temperature and densities [142], and the study 
of chiral solitons [143]. 

In this chapter we present our first results from a calculation of the relativistic Faddeev 
equation for a non-local NJL model [144], based on the co variant diquark-quark picture of 
baryons [118, 119, 120, 121, 122, 123, 124, 125]. Such an approach has been extensively 
employed to study baryons in the local NJL model, see, e.g., Refs. [118, 119, 120, 121]. 
We include both scalar and the axial-vector diquark correlations. We do not assume a 
special form for the interaction Lagrangian, but we rather treat the coupling in the diquark 
channels as free parameters and consider the range of coupling strengths which lead to a 
reasonable description of the nucleon. We construct diquark and nucleon solutions and 
study the possible implications of the quark confinement in the solutions. Due to the 
separability of the non-local interaction, the Faddeev equations can be reduced to a set of 
effective Bethe-Salpeter equations. This makes it possible to adopt the numerical method 
developed for such problems in Refs. [122, 123, 124, 125]. 

6.2 A non-local NJL model 

We consider a non-local NJL model Lagrangian with SU (2)/ x SU (3)c symmetry. 

C — '4'{i^ — m^tl) -\- Ci, (6.1) 

where rric is the current quark mass of the u and d quarks and £/ is a chirally invariant 
non-local interaction Lagrangian. Here we restrict the interaction terms to four-quark 
interaction vertices. 

There exist several versions of such non-local NJL models. Regardless of what version 
is chosen, by a Fierz transformation one can rewrite the interaction in either the qq or qq 
channels, and we therefore use the interaction strengths in those channels as independent 
parameters. For simplicity we truncate the mesonic channels to the scalar (0+, T — G) and 
pseudoscalar (0~, T = 1) ones. The qq interaction is truncated to the scalar (0+, T = 0) 
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and axial vector (1+, T = 1) colour 3 qq channels (the colour 6 channels do not contribute 
to the colourless three-quark state considered here). We parametrise the relevant part of 
interaction Lagrangian as 

>C7 = ^g7,ja{x)ja{x) + 9sJs{x)Js{x) + gaJa{x)Ja{x), 

X)'llj{xi)r^'llj{x3), 

x)^{Xi)[-f5CT2p^]i^ {X^), 

X)1P^{X2)[C-'J,T2(3^]1P{X,). 

x)'iIj{xi) [7;^CTjT2/3^]V'^(x3), 

a;)V'^(x2)[C-S''T2T,/3^]V(x4), (6.2) 

where r„ = (1, i75r). The matrices (3^ = ^/3J2X^{A = 2, 5, 7) project onto the colour 
3 channel with normalisation tr = 35^^' and the t, 's are flavour SU (2) matrices 

with tr (TjTj) = 26ij. The object C = ^7275 is the charge conjugation matrix. 

Since we do not restrict ourselves to specific choice of interaction, we shall treat the 
couplings Qs, Qa and g^, as independent parameters. We assume gT,^s,a > 0, which leads 
to attraction in the given channels (and repulsion in the qq colour octet and qq colour 
antisextet chaimels). The coupling parameter gr^ is responsible for the pions and their 
isoscalar partner a. The coupling strengths Qs and Qa specify the behaviour in the scalar 
and axial-vector diquark chaimel, respectively. 

For simplicity, we assume the form factor f{x — Xi) to be local in momentum space, 
since it leads to a separable interaction 

f{x-x,) = J ^e-(---')-V(p)- (6.3) 

It is exactly this separability that is also present in the instanton liquid model [145]. The 
dressed quark propagator S{k) is now constructed by means of a Schwinger-Dyson equa- 
tion (SDE) in the rainbow-ladder approximation. Thus the dynamical constituent quark 
mass, arising from spontaneously broken chiral symmetry, is obtained in Hartree approx- 
imation^ (the symbol Tr denotes a trace over flavour, colour and Dirac indices and tr d 

'Notice that as we demonstrated in section 4.4, the exchange diagrams (for four-fermion interactions) 
can always be cast in form of direct diagram via a Fierz transformation. This means, the Hartree-Fock ap- 





J d'^Xid^Xsf{x 


- X3)f{xi 


Js{x) = 


J d'^Xid'^Xsflx 


- X3)f{Xi 


Js{x) = 


J d'^X2d'^X4f{x 


- Xi)f{x2 


Ja{x) = 


j d'^xid^X3f{x 


- X3)f{xi 




J d'^X2d^X4f{x 


- Xi)f{x2 
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denotes a trace over Dirac indices only) 



M{p) ^mc + ig. 



'TT, 



fip) J ^Jr[S{k)]f{k), 



(6.4) 



where 



S-\k) = I^-M{k), 
one can simplify this equation by writing M{p) in the form 



(6.5) 



M{p) = m, + {M{0) - m,)f{p). 



(6.6) 



The non-linear equation can then be solved iteratively for M(0). 

Following Ref. [135], we choose the form factor to be Gaussian in Euclidean space, 
/(pe) = exp(— p|,/A^), where A is a cutoff of the theory. If one assumes that A is re- 
lated to the average inverse size of instantons 1/p, then its value parametrises aspects of 
the non-perturbative properties of the QCD vacuum [99] . This choice respects Poincare 
invariance and for certain values of the parameters it leads to quark, but not colour, con- 
finement. For values of M(0) satisfying 



the dressed quark propagator has no poles at real in Minkowski space (p^ + M^(p^) ^ 
0). The propagator has infinitely many pairs of complex poles, both for confining and non- 
confining parameter sets. This is a feature of these models and due care should be taken in 
handling such poles, which can not be associated with asymptotic states if the theory is to 
satisfy unitarity. One should note that the positions of these poles depend on the details of 
the chosen form factor and the cut-off, hence one may regard them as a pathology of the 
regularisation scheme. Since the choice of the cut-off is closely related to the truncation 
of the mesonic channels, (for example, if one allows mixing of channels, the cut-off and 
the positions of poles will change. In Fig. 6.5 we have shown the positions of the first 
poles of the quark propagator for various cutoff. We have examined that in the presence of 
nai mixing, these positions will change, but it follows very similar trend). Even though 
the confinement in this model has no direct connection to the special properties of the 
pion, there is an indirect connection through the determination of the parameters from the 

proximation is equivalent to the Hartree approximation with properly redefining coupling constants. There- 
fore, Hartree approximation is as good as Hartree-Fock one as long as the interaction terms in Lagrangian 
are not fixed by some underlying theory. 




(6.7) 
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Figure 6. 1 : Quark condensate for various sets of parameters as a function of current quark 
mass (right). The dynamical quark mass at zero momentum as a function of coupling (left) 
for parameter set: Afo(O) = 300 MeV, A = 860.58 MeV = 10.98, the mass generation 
occur at critical coupling g'^ = 27.0 [GeV~^]. 




Figure 6.2: A graphical representation of the Bethe-Salpeter equation for the qq T-matrix 
in RPA approximation. The solid lines denote the dressed quark propagators Eq. (6.5) 
and shaded boxes denote meson propagators. 



pionic properties. From the gap equation Eq. (6.4) one can show that dynamical symmetry 
breaking occurs for l/g^^ — < 0, where l/g'^^ = -^^^2—, Nc and Nj are the number 
of colours and flavours, respectively. For g,, > g^' fermions become massive and the 
vacuum accommodates a non- vanishing condensate (V^V^), and consequently there exists 
a massless Nambu-Goldstone boson, see Fig 6.1. 



6.3 Meson channel 



The quark-antiquark T-matrix in the pseudoscalar channel can be solved by using the 
Bethe-Salpeter equation in the random phase approximation (RPA), as shown in Fig. 6.2, 
seeRef. [135]. 
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+ 




Figure 6.3: One-pion-to- vacuum matrix element in RPA, contributing to the weak pion 
decay. The lines are as defined in Fig. 6.2. The wavy line denotes a weak decay. 



x5{pi +P2 -Ps-Pa), 



[«75Ti]/(p3)/(P4) 



(6.8) 



where 



J,(g') = iTr 
= 6i 



f{k)^5TiS{k)^5riS{q + k)f{q + k), 



tr Db5S{k)j,S{k + q)]f{k)f{q + k), 



{2n) 



(6.9) 



where q denotes the total momentum of the qq pair. The pion mass corresponds to 
the pole of T-matrix. One immediately finds that m-,^ = if the current quark mass rric 
is zero, in accordance with Goldstone's theorem. The residue of the T-matrix at this pole 
has the form 

^"(pi,P2) = ig.,,[tc (8) r" (8) 75]/(pi)/(p2), (6.10) 

where g^^qq is the pion-quark-antiquark coupling constant and is related to the correspond- 
ing loop integral J^^ by 



-2 _ 
9nqq 



dq"^ 



(6.11) 



Notice that Z — g^^^ can be regarded as a pion wavefunction renormalisation constant. 
For space-time dimension D = A, one can show a (for the local NJL model we 
have Z^'^ oc LnA), therefore in the continuum limit A — > cxd we have Z = which is 
precisely the compositeness condition [146]. In this extreme limit pions become pointlike. 
The cutoff for spacetime D = A can be removed only at the expense of making the 
theory trivial in the continuum limit. It has been shown that for four-fermion theories the 
renormalisability, nontriviality and compositeness are intimately related [130, 147]. 

The pion decay constant /tt is obtained from the coupling of the pion to the axial- 
vector current. Notice that due to the non-locality the axial-vector current is modified 
[135, 148] and consequently the one-pion-to-vacuum matrix element gets the additional 
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Table 6.1: The parameters for the sets A and B, fitted to = 92.4 MeV and ttItt = 139.6 
MeV. Resulting values of the dynamical quark mass M(0) are also shown. 

Parameter set A set B 
M(0) (MeV) 297.9 351.6 
Mo(0) (MeV) 250 300 
mc(MeV) 7.9 11.13 
A (MeV) 1046.8 847.8 
g^(GeV~^) 31.6 55.80 



Table 6.2: The first two sets of poles of the quark propagator (in magnitude) in the 
Minkowski frame. 

set A set B 
±391 MeV ±408 ± 238i MeV 
±675 MeV ±1575 ± 307i MeV 



contribution shown in Fig. 6.3. This new term is essential in order to maintain Gell- 
Mann-Oakes-Renner relation [135] and makes a significant contribution. The pion decay 
constant is given by 

^ ^/ (0i^[^^5|(5(P-))75r„(5(p+))]/(p-)/(p+) 

+ l|/ / (0iTr[T/'^(p_,p^)5(p_)75r„5(p^)] 

x[f{k) {f{p^) + f{p.)) - f{p+)f{p.)f{k) {f{k + q) + f{k - q)X^.l2) 

where T4(P-,P+) is defined in Eq. (6.10), with notation p± — p±^q. 

Our model contains five parameters: the current quark mass rric, the cutoff (A), the 
coupling constants Qj, , gg and We first fix g^, and the current quark mass rric for 
arbitrary values of A by fitting and to their empirical values. In this way, we 
can consider the entire parameter space of the model. The corresponding solution of gap 
equation are shown in Fig. 6.4. In the left panel the constituent mass M at zero momentum 
is shown as a function of the cutoff. It is obvious that for very small cutoff, there is no 
solution for the gap equation. On the right panel of Fig. 6.4, we show the corresponding 
values of the quark condensate. The quark condensate (^'V') — iT^S{0) is closely related 
to the gap equation Eq (6.4). In the latter there appears an extra form factor inside the 
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Figure 6.4: The dynamical quark mass at zero momentum (a) and quark condensate (b) 
with respect to the cutoff for fixed = 93 MeV and = 140 MeV. The solid line 
denotes the result when is not zero and the dashed line denotes the chiral limit rric = 0. 



loop integral. The quark condensate with non-zero current-quark mass is quadratically 
divergent and is regulated by a subtraction of its perturbative value. These values can 
fall within the limits extracted from QCD sum rules 190MeV < -{qq)^/^ < 260MeV 
at a renormalisation scale of 1 GeV [149] and lattice calculation [150], having in mind 
that QCD condensate is a renormalised and scale-dependent quantity. In right-hand side 
of Fig. 6.1 we show the quark condensate with respect to current quark mass for various 
parameter sets, it is seen that for large coupling (or large M(0)), the magnitude of quark 
condensate decreases as the current quark mass increases. This feature is consistent with 
the behaviours of lattice and sum rule results. 

We analyse two sets of parameters, see Table 6.1. The set A is a non-confining param- 
eter set, while set B leads to quark confinement (i.e., it satisfies the condition Eq. (6.7). 
The quark condensate in the chiral limit is — (207MeV)^ and — (186 MeV)'^ for sets 
A and B, respectively. At non-zero current quark mass one obtains — (215MeV)'^ and 
— (IQlMeV)^ for sets A and B respectively. The position of the quark poles are given in 
Table 6.2 for two sets of parameters. In Fig. 6.5 we show the position of the first pole of 
quark propagator with respect to various cutoff, which indicates that for large cutoff, we 
have only real poles and we have indeed complex poles for reasonable range of cutoff. 
The real part of the first pole of dressed quark propagator can be considered in much the 
same as the quark mass in the ordinary NJL model. Since we do not believe in on-shell 
quarks or quark resonances, this is also a measure for a limit on the validity of the theory. 
The real part of the first quark propagator pole is larger than the constituent quark 
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Figure 6.5: the position of the first poles of quark propagator in complex plane for various 
cutoff A (MeV) for fixed = 93 MeV and = 140 MeV For every point, there 
is a complex conjugate partner. The imaginary part and real part denote and m^, 
respectively. Some values of cutoff are given on the plot. 



mass at zero momentum M(0), as can be seen in Table 6.1. 

As we will see the mass appears as an important parameter in diquark and nucleon 
solution rather than the constituent quark mass. The same feature has been seen in the 
studies of the soliton in this model, where determines the stability of the soliton [143]. 
In contrast to the local NIL model, here the dynamical quark mass Eq. (6.6) is momentum 
dependent and follows a trend similar to that estimated from lattice simulations [139]. 
Although this is less fundamental since one is free to choose the form factor, nevertheless 
the quark mass is a gauge dependent object and is not directly observable. 

The parameters gs and ga are yet to be determined, we shall treat them as free parame- 
ters, which allows us to analyse baryon solutions in terms of a complete set of couplings. 
The coupling-constant dependence appears through the ratios = gs/g-K and = ga/g-K- 
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6.4 Diquark channels 



In the rainbow-ladder approximation the scalar qq T-matrix can be calculated from a very 
similar diagram to that shown in Fig. 6.2 (the only change is that the anti-quark must be 
replaced by a quark with opposite momentum). It can be written as 

T(P1,P2,P3,P4) = f{Pl)f{P2)b,CT2(3^]T{q)[C-'^,T2(3^]f{p,)f{p,) 

X5(pi +P2 -P3 -P4), (6.13) 

with 

r(g) = T^^^T^y (6.14) 

where ? = pi + P2 = Pa + P4 is the total momentum of the qq pair, and 

Uq') = ^Tr|^f(-A;)[75^^r2/?^]5(-A;)^[C-ST2/?1% + ^)f(g + A;), 
= Q^J-0-,trnb,S{k)^,S{k + q)]f{k)f{q + k). (6.15) 

In the above equation the quark propagators S{k) are the solution of the rainbow SDE 
Eq. (6.5). The denominator of Eq. (6.14) is the same as in the expression for the pion 
channel, Eq. (6.8), if gg — Qn- One may thus conclude that at = 1 the diquark and pion 
are degenerate. This puts an upper limit to the choice of r^, since diquarks should not 
condense in vacuum. 

One can approximate r(g) by an effective diquark exchange between the external 
quarks, and parametrise T{q) around the pole as 

T{q) = 2igl^^V\q)D{q), D-\q) ^ q' - MJ„ (6.16) 

where Mds is the scalar diquark mass, defined as the position of the pole of r(g). The 
strength of the on-shell coupling of scalar diquark to quarks, gdsqq is related to the polari- 
sation operator Jg by 

9dsqq = ^\q^=Ml, (6-17) 

and V^{q) is the ratio between the exact T-matrix and on-shell approximation. It is obvi- 
ous that we should have (q) \q2=Ml = 1- 

Here, there is no mixing between the axial-vector diquark channel with others, there- 
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Figure 6.6: The denominator of the diquark T matrix for the longitudinal and transverse 
axial vector channel, for parameter set A at = 0.44. Note that there is no longitudinal 
pole. 



fore in the same way one can write the axial- vector diquark T-matrix in a similar form 

T(pi,P2,P3,P4) = f{pi)f{P2) [i,CnT2f3''y{q) [C-'^.T2ni3^]f{ps)f{p,) 

X5(pi +P2 -P3 -P4), (6.18) 



with 



9 



+ 



(6.19) 



Ll + ^„JJ(g2) l + gjj^{q^) 

Here we prefer to decompose the axial polarisation tensor into longitudinal and transverse 
channels as well, 



zTr 



6i 



f{~k) [rCr.r^P^] S{-kf [C-^Yr^nP^] S{q + k)f{q + k), 



tr D[YS{k)YS{k + q)]f\k)f\q + k) 



Jl{q'W - q'q'/q') + Jaiq')q'q''/q'- 



(6.20) 



We find that the longitudinal channel does not produce a pole (see Fig. 6.6), and thus the 
bound axial- vector diquark solution corresponds to a pole of the transverse T-matrix. The 
transverse component of r'^'^(g) matrix is now approximated by as. 



T^''{q)^2zgl^^V^{q)D^^^{q), 



D^^{q) 



9 



q'^q" /q^ 



q' - Ml 



(6.21) 
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Figure 6.7: The scalar and axial- vector diquark mass as a function of and r^, respec- 
tively, for both parameter sets. The dotted and broken-dotted lines lines denote the quark- 
quark pseudo-threshold for set A and B, respectively. 



where includes the off-shell contribution to the T-matrix. The coupling constant 

9daqq IS related to the residue at the pole of the T-matrix, 

9daqq = ^U^=Mj^- (6.22) 



6.4.1 Diquark Solution 

The loop integrations in Eq. (6.15, 6.20) are evaluated in Euclidean space^. For the current 
model, the usual analytic continuation of amplitudes from Euclidean to Minkowski space 
can not be used. This is due to the fact that quark propagators of the model contain 
many poles at complex energies leading to opening of a threshold for decay of a diquark 
(or meson) into other unphysical states. Any theory of this type need to be equipped 
with an alternative continuation prescription consistent with unitarity and macrocausality. 
Let us define a fictitious two-body threshold as twice m^. For a confining parameter 
set, each quark propagator has a pair of complex-conjugate poles. Above the two-body 
pseudo-threshold < — 4(m^)^, where q is meson (diquark) momentum, the first pair of 
complex poles of the quark propagator has a chance to cross the real axis. According to 
the Cutkosky prescription [151], if one is to preserve the unitarity and the microcausality, 
the integration contour should be pinched at that point. In this way, one can ensure that 

^We work in Euclidean space with metric g^" — S^" and a hermitian basis of Dirac matrices {7^, ji^} = 
2(5^iy, with a standard transcription rules from Minkowski to Euclidean momentum space: — > ik^, 

k^' ^ -fc^ 
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Figure 6.8: The scalar diquark-quark-quark coupling as a function of r^. The dotted and 
broken-dotted lines indicate the quark-quark pseudo-threshold for set A and B, respec- 
tively. 

there is no spurious qq (or qq) production threshold, for energies below the next pseudo- 
threshold, i.e. twice the real part of the second pole of the quark propagator. Note that 
it has been shown [152] that the removal of the qq pseudo-threshold is closely related to 
the existence of complex poles in the form of complex-conjugate pairs. Since there is no 
unique analytical continuation method available for such problems, any method must be 
regarded as a part of the model assumptions [135, 141, 152]. Here, we follow the method 
used inRef. [135]. 

We use the parameter sets determined in the mesonic sector shown in table 6.1. Our 
numerical computation is valid below the first qq pseudo-threshold. Note that the longi- 
tudinal polaris ability Jf'(g) defined in Eq. (6.20) does not vanish here. However this term 
can probably be ignored since it does produce any poles in the T-matrix, and moreover 
there is no conserved current associated to this channel. We find that for a wide range 
of and Ta, for all parameter sets, a bound scalar and axial-vector diquark exist (the 
results for additional sets can be found in [32]). This is in contrast to the normal NJL 
model where a bound axial- vector diquark exists only for very strong interaction [120]. 
The diquark masses for various values of and are plotted in Fig. 6.7. As already 
pointed out, the scalar diquark mass is equal to the pion mass at = 1. It is obvious from 
Fig. 6.7 that for = the axial-vector diquark is heavier than the scalar diquark, and 
consequently is rather loosely bound. For very small and one finds no bound state 
of either diquark. This kind of diquark confinement is due to the screening effect of the 
ultraviolet cutoff and can not be associated with confinement in QCD which originates 
from the infrared divergence of the gluon and ghost propagators. Having said that, it is 
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Figure 6.9: The ratio of the on- shell approximation compared to the exact diquark T- 
matrix for the various scalar and axial diquark masses. 



possible that real diquark confinement may arise beyond the ladder approximation [154]. 
There, in order to preserve Goldstone's theorem at every order, we must include addi- 
tional terms in the interaction. Although these new terms should have minimal impact on 
the solutions for the colour-singlet meson channels, they can provide a repulsive contri- 
bution to the colour-antitriplet diquark channels which removes the asymptotic-diquark 
solutions from the spectrum. This would indicate that diquark confinement is an inde- 
pendent phenomenon and is not related to the particular realisation of quark confinement. 

One should note that the nucleon bound state in the diquark-quark picture does not 
require asymptotic-diquark states since the diquark state is merely an intermediate device 
which simplifies the three-body problem. Nevertheless, evidence for correlated diquark 
states in baryons is found in deep-inelastic lepton scatterings and in hyperon weak decays 
[157]. At the same time, diquarks appear as bound states in many phenomenological mod- 
els. It is puzzling that diquarks are even seen in lattice calculations [155, 156]. In contrast 
to our perception of QCD colour confinement, the corresponding spectral functions for 
these supposedly confined objects in the colour anti-triplet channel are very similar to 
mesonic ones [156]. 

In Fig. 6.8, we show the scalar diquark-quark-quark coupling defined in Eq. (6.17) 
with respect to various scalar diquark couplings. A pronounced change in behaviour 
around the quark-quark pseudo-threshold is observed in the confining set B, and this 
seems to justify our emphasis on the pseudo-threshold defined by twice the real part of 
the quark pole. 

Next we study the off-shell behaviour of the diquark T-matrix. In Fig. 6.9 we show 
the discrepancy between the exact T-matrix and the on-shell approximation V^'°'{q). At 
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the pole we have by definition that V^'"'{q)\q2=M^^ = 1- We see that elsewhere the off- 
shell contribution is very important due to the non-locality of our model. We find that the 
bigger the diquark mass is, the bigger the off-shell contribution. The off-shell behaviour 
of the scalar and the axial-vector channel for both parameter sets A and B are rather 
similar. 



6.5 Three-body sector 

In order to make three-body problem tractable, we discard any three-particle irreducible 
graphs (this is sometimes called the Faddeev approximation). The relativistic Faddeev 
equation can be then written as an effective two-body BS equation for a quark and a 
diquark due to the locality of the form factor in momentum space (see Eq. (6.3)) and 
accordingly the separability of the two-body interaction in momentum-space. We adopt 
the formulation developed by the Ttibingen group [122, 123, 124] to solve the resulting 
BS equation. In the following we work in momentum space with Euclidean metric. The 
BS wave function for the octet baryons can be presented in terms of scalar and axial vector 
diquarks correlations, 

^(p,P)^(P,.)=(^^|^'^;J^j)«(P,.), (6.23) 

where u{P, s) is a basis of positive-energy Dirac spinors of spin s in the rest frame. The 
parameters p = {1 — ifjpi ~ r]{pj + pk) and P = Pi + pj + pk are the relative and total 
momenta in the quark-diquark pair, respectively. The Mandelstam parameter rj describes 
how the total momentum of the nucleon P is distributed between quark and diquark. 

One may alternatively define the vertex function associated with ip{p, P) by amputat- 
ing the external quark and diquark propagators (the legs) from the wave function; 

0(p, P) = S-\p,)b-\p,) ( ^'1^1 ) , (6.24) 

with 

where D{p), D^''{p) and S{p) are Euclidean versions of the diquark and quark prop- 
agators which are obtained by the standard transcription rules from the expressions in 
Minkowski space, Eqs. (6.16,6.21) and Eq. (6.5), respectively. The spectator quark mo- 
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mentum pq and the diquark momentum pd are given by 



Pq = 7]P+ p, 



(6.26) 
(6.27) 



with similar expressions for kq^a, where we replace p by /c on the right-hand side. In 
the ladder approximation, the coupled system of BS equations for octet baryon wave 
functions and their vertex functions takes the compact form. 



where K^^{p, k; P) denotes the kernel of the nucleon BS equation representing the ex- 
change quark within the diquark with the spectator quark (see Fig. 5.1), and in the colour 
singlet and isospin ^ channel we find (see Ref. [120]) 



where x and (and their adjoint x and x^) stand for the Dirac structures of the scalar 
and the axial- vector diquark-quark-quark vertices and can be read off immediately from 
Eqs. (6.13, 6.16) andEqs. (6.18, 6.21), respectively. Therefore we have 



x'{Pi,kd) = gdsgg{l'C)y/ 2V^f ip, + (1 - a)kd)f{-pi + akg), 
x'^iPuka) = 9da,,{l''C)^j2V-{kd)f{p^ + (1 - a)kd)f{-pi + aka). (6.30) 



We have used an improved on- shell approximation for the contribution of diquark T- 
matrix occurring in the Faddeev equations. Instead of the exact diquark T-matrices we 
use the on-shell approximation with a correction of their off-shell contribution through 
V^'"-{p). What is neglected is then the contribution to the T-matrix beyond the pseudo- 
threshold. As we will see this approximation is sufficient to obtain a three-body bound 
state. In order to evaluate the structure of the diquark T-matrix completely, one nor- 
mally employs the dispersion relation, however, this is not applicable here, due to non- 
analyticity of the diquark T— matrix. Notice, as we already pointed out for the confining 
set B, we do not have qq continuum, however, there exists many complex poles beyond the 
pseudo-threshold which might be ignored, provided that they lie well above the energies 
of interest. 




(6.28) 



K'''ip,k;P) 



-3 



xHPi,kd)S^{q)x^iP2,Pd) 
VSx'(pi,kd)S^(q)r(P2,Pd) 



V3x''{pi,kd)S^{q)x'{p2,Pd) 

-X''(pi,kd)S^(q)r(p2,Pd) 

(6.29) 



The relative momentum of quarks in the diquark vertices x and x^ are defined as 
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^ p + k/2 - {1 - 3ri)P/2 and = -k - p/2 + (1 - 2>r])P/2, respectively. The 
momentum kd of the incoming diquark and the momentum pd of the outgoing diquark are 
defined in Eq. (6.27) (see Fig. 5.1). The momentum of the exchanged quark is fixed by 
momentum conservation at g = —p — /c + (1 — 2r])P. 

It is interesting to note the non-locality of the diquark-quark-quark vertices naturally 
provides a sufficient regularisation of the ultraviolet divergence in the diquark-quark loop. 
In the expressions for the momenta we have introduced two independent Mandelstam pa- 
rameters T], a, which can take any value in [0, 1]. They parametrise different definitions of 
the relative momentum within the quark-diquark (77) or the quark-quark system (a). Ob- 
servables should not depend on these parameters if the formulation is Lorentz covariant. 
This means that for every BS solution ip{p, P; rji, ai) there exists a equivalent family of 
solutions. This provides a stringent check on calculations, see the next section for details. 

We now constrain the Faddeev amplitude to describe a state of positive energy, positive 
parity and spin s = 1/2. The parity condition can be immediately reduced to a condition 
for the BS wave function: 

/^5(p,P)\^/ 7V(P,^)7' ^\P,P) \ ..... 

V r{p. p) ) V I'^'vriP. p)i' ) V -r{p. p) ) ' ^ ^ 

where we define p = A-pp and P = A-pP, with Ai^ = diag(— 1, —1,-1,1). In order to 
ensure the positive energy condition, we project the BS wave function with the positive- 
energy projector A+ = {1 + f), where the hat denotes a unit four vector (in rest frame 
we have P = P/iM). Now we expand the BS wave function P) in Dirac space 
r e {1, 75, 7^^, 757^^, a^"}. The above-mentioned conditions reduce the number of inde- 
pendent component from sixteen to eight, two for the scalar diquark channel. Si, {i — 1,2) 
and six for the axial-diquark channel, Ai,{i — 1, ...6). The most general form of the BS 
wave function is given by 

^P\p,P) = {s,-ii^S2)A-^, 

V{p. P) = (iP^i>TA^ + PM2 - p'^T^i^A^ + ip^A4 + (p^^T - It) ^5 

-(^7t^t + «Pt)^6)75A+. (6.32) 

Here we write 7^ = 7^" — fP^. The subscript T denotes the component of a four- vector 
transverse to the nucleon momentum, pt — p — P{p.P). In the same way, one can 
expand the vertex function in Dirac space, and since the same constraints apply to the 
vertex function, we obtain an expansion quite similar to Eq. (6.32), with new unknown 
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coefficients Sj and Aj which are substituted the coefficients Si and Ai, respectively. The 
unknown scalar function Si{E>i) and Ai{Ai) depend on the two scalars which can be built 
from the nucleon momentum P and relative momentum p, z = P.p = cos a; (the cosine 
of the four-dimensional azimuthal angle of p^) and p"^. 

In the nucleon rest frame, one can rewrite the Faddeev amplitude in terms of tri- 
spinors each possessing definite orbital angular momentum and spin [122]. It turns out 
that these tri-spinors can be written as linear combinations of the eight components de- 
fined in Eq. (6.32). Thus from knowledge of Si and A^, a full partial wave decomposition 
can be immediately obtained [122]. Notice that although the diquarks are not pointlike 
objects here, they do not carry orbital angular momentum i. e. L'^x^'^{q) — 0- This is 
due to the fact that the off-shell contribution V^'°'{q) is a function of scalar q'^. Moreover, 
the form factor in our model Lagrangian is also scalar, hence the total momentum de- 
pendent part of the diquark-quark-quark vertices are scalar functions and carry no orbital 
angular momentum. Therefore, the partial wave decomposition obtained in Ref. [122] for 
pointlike diquarks can be used here. Note that no such partial wave decomposition can 
be found if one uses the BS vertex function 0^'^ since the axial- vector diquark propagator 
mixes the space component of the vertex function and time component of the axial-vector 
diquark. 



6.5.1 Numerical method for the coupled BS equations 

For solving the BS equations we use the algorithm introduced by Oettel et al [125]. The 
efficiency of this algorithm has already been reported in several publications, see for ex- 
ample Refs. [122, 123, 124]. We will focus here only on the key ingredients of this 
method. The momentum dependence of quark mass in our model increases the complex- 
ity of the computation significantly. 

As usual, we work in the rest frame of the nucleon P — (0, iM^). In this frame we 
are free to chose the spatial part of the relative momentum p parallel to the third axis. 
Thus the momenta p and k are given by 

p^ = b|(o,o,vT^^)_ 

A;^ = |/c|(sin^'sin0Vl - ^'^^sin^'cos^Vl - z'^cos^Vl - z''^,z'), (6.33) 
where we write z — cos a; and z' — cos a;'. The wave function Eq. (6.32) consists of 
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2 X 2 -blocks in Dirac space can be simplified to 



iasMP^z) 0\ j,2rp^^( ia2Mf,z) 

The great advantage of this representation is that the scalar and the axial-vector compo- 
nents are decoupled. Therefore the BS equation decomposes into two sets of coupled 
equations, two for the scalar diquark channel and six for the axial diquark channel. We 
expand the vertex (wave) functions in terms of Chebyshev polynomials of the first kind, 
which are closely related to the expansion into hyperspherical harmonics. This decompo- 
sition turns out to be very efficient for such problems [122, 123, 124, 125]. Explicitly, 



Ffip^z) = Y.'"Fr\p')Tn{zl 

n=0 

^max 

Ft{p\z) = ^ i'^i^^('")(/)T^(^), (6.35) 

m=0 

where T„(^) is the Chebyshev polynomial of the first kind. We use a generic notation, the 
functions Ff{ and Ff) substituting the function S^, A^ (and Sj, Aj), 

§1,2 ^ Fi^2, A1...6 ^ Ft^. (6.36) 

We truncate the Chebyshev expansions involved in F'^ and Ff at different orders rimax and 
rrimax, respectively. We also expand the quark and diquark propagators into Chebyshev 
polynomials. In this way one can separate the P • p and P ■ k dependence in Eqs. (6.28- 
6.24). Using the orthogonality relation between the Chebyshev polynomials, one can 
reduce the four dimensional integral equation into a system of coupled one-dimensional 
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Figure 6.10: The dependence of the nucleon mass on the Mandelstam parameter rj for a 
few values of the cut-off on the expansion. Here we use set B, with Mds = 725 MeV and 
Mda = 630 MeV. The variables rji and ri2 denote the position of the singularities. 



equations. Therefore one can rewrite Eqs. (6.24, 6.28) in the matrix form 

8 fT^max 

j=l m=0 




(6.37) 



Here g^-^ and H^"^ are the matrix elements of the propagator and the quark exchange 
matrices, respectively. The indices n, m give the Chebyshev moments and i, j denote the 
individual channels. To solve Eq. (6.37), we first rewrite it in the form of linear eigenvalue 
problem. Schematically 

X{P^)ip = K{P')ip, (6.38) 

with the constraint that A(P^) = 1 at = — M^. This can be used to determine the 
nucleon mass iteratively. 

As already pointed out the BS solution should be independent of the Mandelstam 
parameters r],a. As can be seen in Fig. 6.10, there is indeed a large plateau for the 77 
dependence if we use a high cut-off on the Chebyshev moments. The limitations on 
the size of this area of stability can be understood by considering where the calculation 
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contains singularities due to quark and diquark poles. 



Tj e 



1 

1 - 



Mat' Mat. 
Mn' Mn 



if 
if 



Mds < Mda, 
Mda < Mds. 



(6.39) 



A similar plateau has been found in other applications [122, 123, 124]. [Other complex 
poles lies out side the minimal region given in Eq. (6.39).] The singularities in the quark- 
exchange propagator put another constraint on the acceptable range of ri;r] > |(1 — 
No such constraint exists for a, which relates to the relative momentum between two 
quarks. To simplify the algebra we take a — 1/2. 

In what follows we use a momentum mesh of 60 x 60 for p, k, mapped in a non-linear 
way to a finite interval. In the non-singular regime of Mandelstam parameter r] Eq. (6.39), 
the Faddeev solution is almost independent of the upper limit on the Chebyshev expan- 



sion, and for rrir, 



10, Tljymj, 



12, see the Fig. 6.10, this seems to be satisfied. This 



limit is some higher than the reported values for simple models [122, 123, 124, 125]. 



6.5.2 Nucleoli Solution 

In order to understand the role of the axial diquark in nucleon solution, we first consider 
the choice ra = 0. For this case we find that the non-confining set A can not generate a 
three-body bound state in this model without the inclusion of the off-shell contribution. 
For the confining set B one also has to enhance the diquark-quark-quark coupling gdsqq 
by a factor of about 1.73 over the value defined in Eq. (6.17) (as we will show, this extra 
factor is not necessary when the axial-vector is included). The situation is even more 
severe in the on-shell treatment of the local NIL model, since one needs to include the 
full qq continuum contribution in order to find a three-body bound state when the axial- 
vector diquark channel is not taken into account [118]. 

As can be seen from Fig. 6.7 a decrease in leads to a larger diquark mass, and an 
increase in the off-shell contribution to the qq T-matrix (see Fig. 6.9). It is this off-shell 
correction is need for a bound nucleon. 

The nucleon result is shown in Fig. 6.11. We also show a fictitious diquark-quark 
threshold defined as Mds + ^'r- The nucleon mass can be seen to depend roughly linearly 
on the scalar diquark mass. A similar behaviour is also seen in the local NIL model [119]. 
Increasing the diquark mass (or decreasing r^) increases the nucleon mass, i.e. the scalar 
diquark channel is attractive. In order to obtain a nucleon mass of 940 MeV, we need 
diquark mases of 608 MeV and 623 MeV for set A and B, respectively. The corresponding 
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Figure 6.11: The nucleon mass without inclusion of the axial diquark channel. The dotted 
lines indicate the diquark-quark threshold. All values are given in GeV. 



nucleon binding energy measured from the diquark-quark threshold are 56 MeV and 91 
MeV for set A and B, respectively, compared to the binding of the diquarks (relative to 
the qq pseudo threshold) of about 174 and 193 MeV for set A and B, respectively. Such 
diquark clustering within the nucleon is also observed in the local NJL model [119], and is 
qualitatively in agreement with a instanton model [158] and lattice simulations [159]. The 
nucleon solutions for sets A and B behave rather differently with respect to the diquark- 
quark threshold, see Fig. 6.12. This indicates that for non-confining set A, the nucleon 
solution is rather sensitive to the diquark-quark threshold and tends not avoid it. However, 
for confining set, since there is no well-defined threshold, this tendency is absent and as 
we approach to the diquark-quark threshold, the nucleon binding energy decreases and 
can approach to zero. 

Next we investigate the effect of the axial- vector diquark channel on nucleon solution. 
We find that the axial-vector diquark channel contributes considerably to the nucleon 
mass^ and takes away the need for the artificial enhancement of the coupling strength 
for set B. In Figs. 6.13, 6.14 we show the nucleon mass as a function of the scalar and 
axial-vector diquark mass. Similar to the scalar diquark channel, we define the axial- 
vector diquark-quark threshold as Mda + We see that as one increase the axial- vector 
diquark (and scalar diquark) masses, the qq interaction is weakened and consequently the 
nucleon mass is increases. Therefore the contribution of the axial-vector channel to the 
nucleon mass is also attractive. 



^Note however that neglecting TrA^-loops may lead to a quantitative overestimate of the axial-vector 
diquark role in the nucleon [160]. 
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Figure 6.12: The nucleon binding energy measured from the diquark-quark threshold with 
respect to scalar diquark mass. Only scalar diquarks are included in the calculation The 
smooth behaviour of solid line indicates the absence of threshold effects for the confining 
case. 



In Fig. 6.15 we plot the parameter space of the interaction Lagrangian with variable 
Ts and Ta which leads to the nucleon mass = 0.940 GeV. The trend of this plot for 
the non-confining set A is very similar to the one obtained in the local NJL model [120] 
(although we use a different parameter set) and roughly depends linearly on the ratios 
and ra, 

M^[0.940GeV] = -r, - 0.94r„ + 1.2. (6.40) 

Therefore, any interaction Lagrangian with and which satisfies the relation + 
0.94ra = 0.30 gives a nucleon mass at about the experimental value. This relation shows 
how interaction is shared between scalar and axial- vector channels. If the scalar diquark 
interaction is less than 0.14, we need the axial- vector interaction to be stronger than the 
scalar diquark channel > in order to get the experimental value of nucleon mass. For 
set B, as we approach to = 0, the curve bends upward, reflecting the fact that we have 
no bound state with only scalar diquark channel. In Fig. 6.15 we see for the confining 
set B that the interaction is again shared between the scalar and the axial-vector diquark 
and for small < 0.19 one needs a dominant axial- vector diquark channel ra > r^. It 
is obvious that the axial- vector diquark channel is much more important in the confining 
than the non-confining phase of model. 
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Figure 6.13: The nucleon mass as a function of the scalar diquark mass. The dotted lines 
are the scalar diquark-quark threshold, the broken lines are the axial vector diquark-quark 
threshold. All values are given in GeV. 
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Figure 6.14: The nucleon mass as a function of the scalar diquark mass for various axial 
vector diquark masses for both parameter sets. The scalar diquark-quark threshold are 
shown by the doted lines. 
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Figure 6.15: Range of parameters (r^, ra) where we find a nucleon mass of 940 MeV. 
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Table 6.3: Diquark masses and coupling of diquarks to quarks obtained for — 0.940. 
All masses are given in MeV. E^siEda) denote the binding energy of diquarks in the 
nucleon, Eff{Eff) denote the binding energy of the nucleon measured from scalar (axial) 
diquark mass^ 





Set Al 


Set A 
Set A2 


Set A3 


SetBl 


SetB 
SetB2 


SetB3 


as 


115 


748 


698 


802 


705 


609 


Qdsqq 


0.74 


0.83 


1.04 


0.73 


1.31 


1.79 




0.09 


0.12 


0.17 


0.06 


0.14 


0.24 


Eds 


7 


34 


84 


14 


111 


207 




226 


199 


149 


270 


173 


77 


Mda 


705 


725 


775 


604 


660 


725 


Qdaqq 


1.08 


0.98 


0.79 


1.99 


1.67 


1.28 


Ta 


0.20 


0.17 


0.11 


0.32 


0.23 


0.15 


Eda 


11 


57 


7 


212 


156 


91 




156 


176 


226 


72 


123 


193 


„RMS 

Pi. 


194.88 


181.51 


163.70 


283.99 


232.86 


209.95 



In order to study the implications of the quark confinement for the description of the 
nucleon, we compare in Table 6.3 three representative cases for both the non-confining 
and confining parameter sets, which all give nucleon mass about 940 MeV. The first three 
columns contain results for set A, and the last three columns for the confining set B. 

Given the definition of diquark-quark thresholds, in the presence of both scalar and 
axial- vector diquark channels, the diquarks in the nucleon can be found much more 
loosely bound, although one obtains a very strongly bound nucleon solution near its ex- 
perimental value, see table 6.3. Next we study the nucleon BS wave function for the 
various sets given in Table 6.3. The nucleon wave and vertex function are not physi- 
cal observables, but rather they suggest how observables in this model will behave. In 
Figs. 6.16-21 we show the leading Chebyshev moments of the scalar functions of the nu- 
cleon BS wave function for various sets (Al-3 and Bl-3) which describes the strengths 
of the quark-diquark partial waves with 5* as a total quark-diquark spin and L as a total 
orbital angular momentum. They are normalised to Ff^^\pi) = 1, where pi is the first 
point of the momentum mesh. It is seen that the contribution of higher moments are con- 
siderably small, indicating a rapid convergence of the wave function amplitudes in terms 
of Chebyshev polynomials. In the confining case Figs. 6.19-21 there is a clear interfer- 
ence which is not present in the non-confining Figs. 6.16-19. Therefore in the confining 
case, all wave function amplitudes are somehow shifted to higher relative four-momenta 



13SCHAPTER 6. BARYONS STRUCTURE IN A NON-LOCAL QUARK CONHNING NIL MODEL 

between diquark and quark. In order to understand the role of this interference we obtain 
a mass density for the various channels. This density is defined as 



where p± stands for space component of relative momentum p and D^^{pd) defined in 
Eq. (6.25). This definition corresponds to a diagram occurring for the calculation of 
the isoscalar quark condensate in the impulse approximation [121]. In the above defini- 
tion of the density function, we have integrated over the time component of the relative 
momentum. In this way the density function becomes very similar to its counterpart in 
Minkowski space. Although the above definition of density is not unique, it does provide 
a useful measure of the spatial extent of the wave function (we have examined the possi- 
bility of taking matrix elements of other operators between the BS wave function, since 
this does not lead to any significant effect, the results are not presented here). The results 
are plotted in Figs. 6.22 and 6.23. 

It is noticeable that in various sets, the s-wave is the dominant contribution to the 
ground state. The relative importance of the scalar and the axial diquark amplitude in 
the nucleon changes with the strength of the diquark-quark couplings gdsqqiddaqq) and 
accordingly with rs{ra). We see that in the confining sets, the nucleon density extends 
to higher relative momentum between the diquark and the quark. This indicates a more 
compact nucleon in the confining case. In order to find a qualitative estimation of the 
confinement effect in our model, we calculate = {{p]} — {p±)'^y^'^, the results can 
be found in table 6.3. This can be related to the mean-square radius of the nucleon, if 
we assume minimal uncertainty. We see in the both confining and non-confining cases 
a decrease in p^^ with weakening axial-vector diquark interaction (and consequently 
increasing the scalar diquark interaction strength). If we compare p^^ for the two sets 
A2, B2, which have very similar interaction parameters rs(ra), an increase about 25% 
is found. This effect can not only be associated with the non-locality of our interaction, 
since that is present in both confining and non-confining cases. 

6.6 Summary and Outlook 

In this chapter we investigated the two- and three-quark problems in a non-local NJL 
model. We have truncated the diquark sector to the bound scalar and the axial- vector chan- 
nels. We have solved the relativistic Faddeev equation for this model and have studied the 
behaviour of the nucleon solutions with respect to various scalar and the axial-vector in- 




(6.41) 
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teractions. We have also investigated a possible implication of the quark confinement of 
our model in the diquark and the baryon sector. 

Although the model is quark confining, it is not diquark confining (at least in the 
rainbow-ladder approximation). A bound diquark can be found in both scalar and the 
axial-vector channel for a wide range of couplings. We have found that the off-shell 
contribution to the diquark T-matrix is crucial for the calculation of the structure of the 
nucleon: without its inclusion the attraction in the diquark channels is too weak to form 
a three-body bound state. We have also found that both the scalar and the axial-vector 
contribute attractively to the nucleon mass. The role of axial-vector channel is much more 
important in the confining phase of model. The nucleon in this model is strongly bound 
although the diquarks within nucleon are loosely bound. The confining aspects of the 
model are more obvious in three-body, rather than the two-body sector. By investigating 
the nucleon wave function we showed that quark confinement leads to a more compact 
nucleon. The size of nucleon is reduced by about 25% in confining phase. 

For both confining and non-confining phases, an increase in the scalar diquark channel 
interaction leads to a lower nucleon mass, see Figs. 6.13 and 6. 14, but the mass of the A 
remains unchanged since it does not contain scalar diquarks. In the standard NIL model 
where the axial-vector diquark does not contribute significantly to the nucleon binding 
[127, 161, 162], the difference between the A and nucleon mass is directly related to the 
scalar diquark interaction. In the current model where the axial- vector diquark makes a 
larger contribution to the nucleon mass, therefore a detailed calculation for the delta states 
is needed to understand the mechanism behind the A-A^ mass difference. In the standard 
NIL model this leads to a contradiction, since for an axial- vector interaction which gives 
a reasonable description of nucleon properties, both the axial- vector diquark and more 
importantly the A are unbound [162]. The crucial role of the axial- vector diquark corre- 
lation in the non-local NIL model, especially in the confining phase of the model, indicate 
that this model might do better. 

In order to understand the implications of this model in baryonic sector fully one 
should investigate properties of the nucleon such as the axial vector coupling constant, 
the magnetic moment, etc. On the other hand, the role of quark confinement in this model 
can be better clarified by investigating quark and nuclear matter in this model. One of 
the long standing problem in four-Fermi chiral quark models is the fact that quark/nuclear 
matter does not saturate [163], mainly due to the strongly attractive quark interactions 
responsible for the spontaneous chiral symmetry breaking. Recently, Bentz and Thomas 
[164] have shown that a sufficient strong repulsive contribution can arise if confinement 
effects are incorporated, albeit in the cost of introducing a new parameter into model. 
It was shown that such repulsive contribution can lead to saturation of nuclear matter 
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equation of state. It is indeed of interest to investigate the stability of nuclear matter 
within this quark confining model. Such problems can be studied based on the Faddeev 
approach, see e. g., Ref. [131]. 
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Figure 6.16: Different Chebyshev moments (labeled by n*^) of scalar and axialvector 
(AV) diquark amplitudes of the nucleon BS wave function given by Set Al. 
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Figure 6.17: Different Chebyshev moments (labeled by n*^) of scalar and axialvector 
diquark amplitudes of the nucleon wave function given by Set A2. 
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Figure 6.18: Different Chebyshev moments (labeled by n*^) of scalar and axialvector 
diquark amplitudes of the nucleon wave function given by Set A3. 
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Figure 6.19: Different Chebyshev moments (labeled by n ) of scalar and axialvector 
(AV) diquark amplitudes of the nucleon BS wave function given by Set Bl. 
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Figure 6.20: Different Chebyshev moments (labeled by n*^) of scalar and axialvector 
diquark amplitudes of the nucleon wave function given by Set B2. 
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Figure 6.21: Different Chebyshev moments (labeled by n*^) of scalar and axialvector 
diquark amplitudes of the nucleon wave function given by Set B3. 
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Figure 6.22: Shows the nucleon density (M = 0.940 GeV, Set A) with respect to relative 
momentum between diquark and quark for different set of Al, A2 and A3. 
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Figure 6.23: Shows the nucleon density (M = 0.940 GeV, Set B) with respect to relative 
momentum between diquark and quark for different set of 51, B2 and B3. 
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